NOIP2017总结 & 题解

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中还犯了很多知识性错误，主要原因是学的不扎实，接下来的比赛加油吧。

题解

day1

t1 答案就是x*y-x-y，证明可以见usaco4.1麦香牛块的题解。貌似此题坑到了我校不少选手。

t2 细心模拟即可。

t3 设f(i,j)表示到i点从1~i距离为dis[i]+j的路径个数。对于边u->v，有dis[v]+j_v=dis[u]+j_u+w(u,v)，所以$f(v,j_v)=\sum_{u->v}f(u,dis[v]-dis[u]-w(u,v)-j_v)$，具体实现可以建出反图记忆化搜索，答案为-1当且仅当同一状态在dfs栈中出现了两次。

day2

t1 建图后bfs判联通。

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 the last element and push it to the xth row
            ins(rt[x],1,m-1+q,m-1,las,x);

            del(rt[n+1],1,n+q,x);//push o to the last column
            ins_las(rt[n+1],1,n+q,n,o);
        }
    }
    return 0;
}