day1
t1的结论貌似在哪见过,自己稍微验证了一下貌似没记错就没有管了。
t2一道很好(keng)的模拟题
t3自己做题好慢啊,想出来dp打上去最后几分钟才过了大样例,我写的是记忆化搜索,判-1很好判,没时间加上去了可惜了,不过还是自己做题太慢了。然后由于没拍,不确定自己dp对不对,就特判了k=0粘了一个暴力,没想到暴力错了,70变成50...基础不扎实gg
100+100+50=250 好菜啊t3没A掉
day2
t1看完题就开始打,没注意算距离时减一下就变成两倍了会爆long long,但好像数据中z都>=0 ??
t2写完搜索后发现好像可以把所有点深度变成一个n进制数,可以记忆化,不过打完发现n=12要跑10s+,果断删掉换成爆搜。并没有想到如何状压...
t3打完前50分后有一个想法,打完后发现自己想的是错的,于是老老实实写x=1的30分,但耗费了好多时间,没时间打出个数据结构了,于是用了vector,还以为vector有多快,CCF的老爷机上1e5都跑不过,于是t3就挂了,暴力分都没拿全。
100+70+50=220 太菜了
自己做较难题不行,另外noip中还犯了很多知识性错误,主要原因是学的不扎实,没想到第一次也是最后一次noip就这样结束了(本菜鸡去年没考noip)...
题解
day1
t1 答案就是x*y-x-y,证明可以见usaco4.1麦香牛块的题解。
t2 细心模拟即可。
t3 设f(i,j)表示到i点从1~i距离为dis[i]+j的路径个数。对于边u->v,有dis[v]+jv=dis[u]+ju+w(u,v),所以
f(v,jv)=∑ f(u,dis[v]-dis[u]-w(u,v)-jv) {u | u->v},具体实现可以建出反图记忆化搜索,答案为-1当且仅当同一状态在dfs栈中出现了两次。
day2
t1 模拟。
t2 f1(i,s)表示点i到集合s中点的最小权值,这个可以O(n^2*2^n)处理出。 f2(s1,s2)表示集合s2与s1相连花费的最小代价,O(n*3^n)处理。
f(i,s)表示已经确定了深度<=i的点,且这些点的集合为s,将来要付出的最小代价,枚举下一层选哪些点,f(i,s)=min{ f(i+1,s|s‘)+f2(s,s‘),虽然s‘可能会连向深度<i的点,但如果这样转移了答案会变大,而且如果真要这样连,之前就会转移了的,所以不会计算出错的答案也不会漏掉最优解,这一部分复杂度是O(n*3^n)的。
t3 x=1可以用线段树完成,只需记当前区间未被删除的点的个数,删除第k个数直接查到删去即可,插入同理。
对于一般情况可以对每一排的前m-1个元素建出线段树,然后对最后一列再建出线段树,操作和x=1是一样的,但发现空间开不下,因为这n+1棵线段树最开始叶子节点是等差数列,所以动态开点即可。
附d2t3代码
#include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> #include<algorithm> #include<vector> #include<cmath> #define P puts("lala") #define cp cerr<<"lala"<<endl #define pb push_back #define fi first #define se second #define mkp make_pair #define ln putchar(‘\n‘) using namespace std; typedef pair<int,int> pii; inline void read(int &re) { char ch=getchar();int g=1; while(ch<‘0‘||ch>‘9‘) {if(ch==‘-‘)g=-1;ch=getchar();} re=0; while(ch<=‘9‘&&ch>=‘0‘) re=(re<<1)+(re<<3)+ch-48,ch=getchar(); re*=g; } typedef long long ll; inline void read(ll &re) { char ch=getchar();ll g=1; while(ch<‘0‘||ch>‘9‘) {if(ch==‘-‘)g=-1;ch=getchar();} re=0; while(ch<=‘9‘&&ch>=‘0‘) re=(re<<1)+(re<<3)+ch-48,ch=getchar(); re*=g; } const int N=300050; int n,m,q; int rt[N],ch[N*80][2],siz[N*80],sz=0; ll val[N*80]; ll kth_las(int o,int l,int r,int k) { if(l==r) {return val[o];} int mid=l+r>>1; if(!ch[o][0]) { ch[o][0]=++sz;ch[o][1]=++sz; siz[ch[o][0]]=(siz[o]+1)>>1; siz[ch[o][1]]=siz[o]>>1; if(l==mid) val[ch[o][0]]=(ll)l*m; if(mid+1==r) val[ch[o][1]]=(ll)r*m; } if(k<=siz[ch[o][0]]) return kth_las(ch[o][0],l,mid,k); else return kth_las(ch[o][1],mid+1,r,k-siz[ch[o][0]]); } ll kth(int o,int l,int r,int k,int wh) { if(l==r) {return val[o];} int mid=l+r>>1; if(!ch[o][0]) { ch[o][0]=++sz;ch[o][1]=++sz; siz[ch[o][0]]=(siz[o]+1)>>1; siz[ch[o][1]]=siz[o]>>1; if(l==mid) val[ch[o][0]]=(ll)(wh-1ll)*m+l; if(mid+1==r) val[ch[o][1]]=(ll)(wh-1ll)*m+r; } if(k<=siz[ch[o][0]]) return kth(ch[o][0],l,mid,k,wh); else return kth(ch[o][1],mid+1,r,k-siz[ch[o][0]],wh); } void del(int o,int l,int r,int k) { if(l==r) {siz[o]=0;val[o]=0;return ;} int mid=l+r>>1; if(k<=siz[ch[o][0]]) del(ch[o][0],l,mid,k); else del(ch[o][1],mid+1,r,k-siz[ch[o][0]]); siz[o]=siz[ch[o][0]]+siz[ch[o][1]]; } void ins_las(int o,int l,int r,int x,ll k) { if(l==r) {siz[o]=1;val[o]=k;return ;} int mid=l+r>>1; if(!ch[o][0]) { ch[o][0]=++sz;ch[o][1]=++sz; siz[ch[o][0]]=(siz[o]+1)>>1; siz[ch[o][1]]=siz[o]>>1; if(l==mid) val[ch[o][0]]=(ll)l*m; if(mid+1==r) val[ch[o][1]]=(ll)r*m; } if(x<=siz[ch[o][0]]) ins_las(ch[o][0],l,mid,x,k); else ins_las(ch[o][1],mid+1,r,x-siz[ch[o][0]],k); siz[o]=siz[ch[o][0]]+siz[ch[o][1]]; } void ins(int o,int l,int r,int x,ll k,int wh) { if(l==r) {siz[o]=1;val[o]=k;return ;} int mid=l+r>>1; if(!ch[o][0]) { ch[o][0]=++sz;ch[o][1]=++sz; siz[ch[o][0]]=(siz[o]+1)>>1; siz[ch[o][1]]=siz[o]>>1; if(l==mid) val[ch[o][0]]=(ll)(wh-1ll)*m+l; if(mid+1==r) val[ch[o][1]]=(ll)(wh-1ll)*m+r; } if(x<=siz[ch[o][0]]) ins(ch[o][0],l,mid,x,k,wh); else ins(ch[o][1],mid+1,r,x-siz[ch[o][0]],k,wh); siz[o]=siz[ch[o][0]]+siz[ch[o][1]]; } void wj() { freopen("phalanx.in","r",stdin); freopen("phalanx.out","w",stdout); } int main() { wj(); int i,j,opt,T; read(n);read(m);read(q); for(i=1;i<=n;++i) { rt[i]=++sz;siz[rt[i]]=m-1+q; ch[rt[i]][0]=++sz; ch[rt[i]][1]=++sz; siz[ch[rt[i]][0]]=(siz[rt[i]]+1)>>1; siz[ch[rt[i]][1]]=siz[rt[i]]>>1; } rt[n+1]=++sz;siz[rt[n+1]]=n+q; ch[rt[n+1]][0]=++sz; ch[rt[n+1]][1]=++sz; siz[ch[rt[n+1]][0]]=(siz[rt[n+1]]+1)>>1; siz[ch[rt[n+1]][1]]=siz[rt[n+1]]>>1; for(int cas=1;cas<=q;++cas) { int x,y;read(x);read(y); if(y==m) { ll o=kth_las(rt[n+1],1,n+q,x); del(rt[n+1],1,n+q,x); printf("%lld\n",o); ins_las(rt[n+1],1,n+q,n,o); } else { ll o=kth(rt[x],1,m-1+q,y,x); del(rt[x],1,m-1+q,y); printf("%lld\n",o);//get (x,y) ll las=kth_las(rt[n+1],1,n+q,x);//get last element and push to xth row ins(rt[x],1,m-1+q,m-1,las,x); del(rt[n+1],1,n+q,x);//push o to last column ins_las(rt[n+1],1,n+q,n,o); } } return 0; }