bzoj3175: [Tjoi2013]攻擊裝置&&4808: 馬
阿新 • • 發佈:2018-03-15
ostream 最小 染色 truct ret printf mes stream ans
終於知道為啥網絡流這麽受歡迎了。
其實就是構個圖模板一下的事兒,比較好打是吧。
然後這題網絡流黑白染色(其實感覺上匈牙利更加直接好想啊,但是實際上黑白染色給人感覺就是二分圖)
st連白而ed連黑,流量為1
不能同時出現的就建無限流量的邊
然後sum-最小割
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> using namespacestd; const int dx[8]={-2,-2,-1,1,2,2,1,-1}; const int dy[8]={-1,1,2,2,1,-1,-2,-2}; struct node { int x,y,c,next,other; }a[410000];int len,last[51000]; void ins(int x,int y,int c) { int k1,k2; len++;k1=len; a[len].x=x;a[len].y=y;a[len].c=c; a[len].next=last[x];last[x]=len; len++;k2=len; a[len].x=y;a[len].y=x;a[len].c=0; a[len].next=last[y];last[y]=len; a[k1].other=k2; a[k2].other=k1; } int st,ed,h[51000],list[51000]; bool bt_h() { memset(h,0,sizeof(h));h[st]=1; int head=1,tail=2;list[1]=st; while(head!=tail) { int x=list[head];for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(a[k].c>0&&h[y]==0) { h[y]=h[x]+1; list[tail]=y; tail++; } } head++; } if(h[ed]==0)return false; return true; } int findflow(int x,int f) { if(x==ed)return f; int s=0; for(int k=last[x];k;k=a[k].next) { int y=a[k].y; if(a[k].c>0&&h[y]==h[x]+1&&s<f) { int t=findflow(y,min(a[k].c,f-s)); s+=t;a[k].c-=t;a[a[k].other].c+=t; } } if(s==0)h[x]=0; return s; } int n,color[210][210]; char ss[210]; bool mp[210][210]; int point(int x,int y){return (x-1)*n+y;} int main() { scanf("%d",&n); int br; int sum=0; for(int i=1;i<=n;i++) { scanf("%s",ss+1); for(int j=1;j<=n;j++) { if(ss[j]==‘0‘)mp[i][j]=true, sum++; else mp[i][j]=false; } } //-----sc----- memset(color,0,sizeof(color)); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(j==1)color[i][j]=1-color[i-1][j]; else color[i][j]=1-color[i][j-1]; //---paint_color--- //----init----- st=n*n+1,ed=n*n+2; for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(mp[i][j]==true) { if(color[i][j]==0) ins(st,point(i,j),1); else ins(point(i,j),ed,1); } for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(mp[i][j]==true&&color[i][j]==0) for(int t=0;t<=7;t++) { int ti=i+dx[t],tj=j+dy[t]; if(ti>0&&ti<=n&&tj>0&&tj<=n&&mp[ti][tj]==true) ins(point(i,j),point(ti,tj),999999999); } //----composition------ int ans=0; while(bt_h()==true) { ans+=findflow(st,999999999); } printf("%d\n",sum-ans); return 0; }
bzoj3175: [Tjoi2013]攻擊裝置&&4808: 馬