bzoj 1453 双面棋盘

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<cstdlib>
 5 #include<cmath>
 6 #include<algorithm>
 7 #include<queue>
 8 #include<vector>
 9 #include<map>
10 #define ll long long
11 #define inf 2139062143
12 #define MAXN 220
13 using namespace std;
14 inline int read()
15 {
16     int x=0,f=1;char ch=getchar();
17     while(!isdigit(ch)) {if(ch==‘-‘) f=-1;ch=getchar();}
18     while(isdigit(ch)) x=x*10+ch-‘0‘,ch=getchar();
19     return x*f;
20 }
21 int n,g[MAXN][MAXN],f[MAXN<<2],t[MAXN<<2];
22 int find(int x) {return f[x]==x?x:f[x]=find(f[x]);}
23 struct data {int ans[2],m[MAXN<<1];}tr[MAXN<<2];
24 void solve(int k,int p)
25 {
26     tr[k].ans[g[p][1]]=1,tr[k].ans[g[p][1]^1]=0,tr[k].m[1]=tr[k].m[1+n]=1;
27     for(int i=2,j=1;i<=n;i++)
28     {
29         if(g[p][i]!=g[p][j]) tr[k].ans[g[p][i]]++,j=i;
30         tr[k].m[i]=tr[k].m[i+n]=j;
31     }
32 }
33 void upd(int k,int p)
34 {
35     tr[k].ans[0]=tr[k<<1].ans[0]+tr[k<<1|1].ans[0];
36     tr[k].ans[1]=tr[k<<1].ans[1]+tr[k<<1|1].ans[1];
37     //cout<<k<<" "<<p<<endl;
38     for(int i=1;i<=n*2;i++) f[i]=tr[k<<1].m[i];
39     //for(int i=1;i<=n*2;i++) cout<<i<<" "<<f[i]<<" "<<tr[k<<1].m[i]<<endl;
40     for(int i=1;i<=n*2;i++) f[i+n*2]=tr[k<<1|1].m[i]+n*2;
41     //cout<<f[1]<<" "<<f[2]<<endl;
42     for(int i=1;i<=n;i++) if(g[p][i]==g[p+1][i]&&find(i+n)!=find(i+n*2))tr[k].ans[g[p][i]]--,f[f[i+n]]=f[i+n*2];
43     for(int i=1;i<=n*4;i++)
44     {
45         f[i]=find(i);
46         if(i>n*3) t[f[i]]=i-n*2;
47     }
48     for(int i=1;i<=n;i++) tr[k].m[i]=i;
49     for(int i=1;i<=n;i++) tr[k].m[i+n]=t[f[i+n*3]];
50 }
51 void build(int k,int l,int r)
52 {
53     if(l==r) {solve(k,l);return;}
54     int mid=(l+r)>>1;
55     build(k<<1,l,mid),build(k<<1|1,mid+1,r);
56     upd(k,mid);
57 }
58 void mdf(int k,int l,int r,int x)
59 {
60     if(l==r) {solve(k,l);return;}
61     int mid=(l+r)>>1;
62     if(x<=mid) mdf(k<<1,l,mid,x);
63     else mdf(k<<1|1,mid+1,r,x);
64     upd(k,mid);
65 }
66 int main()
67 {
68     n=read();int a,b;
69     for(int i=1;i<=n;i++)
70         for(int j=1;j<=n;j++) g[i][j]=read();
71     build(1,1,n);
72     int T=read();
73     while(T--) a=read(),b=read(),g[a][b]^=1,mdf(1,1,n,a),printf("%d %d\n",tr[1].ans[1],tr[1].ans[0]);
74 }