标签:
正解貌似是大暴搜?
首先我们考虑这是一个二分图,建立网络流模型后很容易得出一个算法
S->行 容量为Num[X]/2;
行->列 容量为1 且要求(x,y)这个点存在
列->T 容量为Num[Y]/2
这样子跑网络流之后我们就得到了一组解
但是我们考虑输出方案
对于每一行,如果Num[X]为偶数,那么显然输出方案是正确的
但是如果Num[x]为奇数,多出的那个显然既有可能是红的也可能是蓝的
但关键是我们不能确定他是红的或者蓝的,因为他的状态也会影响对应的列
同样,列的考虑也是同理
所以我们对于残量网络,如果Num[X]或者Num[Y]是奇数,那么就连对应的容量为1的边
对残量网络在进行一次网络流,这样就可以输出方案辣
(话说网络流跑40w的点居然一点也不虚)
#include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> #include<algorithm> #include<queue> #define h(i) i #define l(i) i+200000 using namespace std; const int maxn=200010; const int oo=0x7fffffff; int S,T; int n,cntx,cnty; int X[maxn],Y[maxn]; int NX[maxn],NY[maxn]; int idx[maxn],idy[maxn]; bool visx[maxn],visy[maxn]; int h[maxn<<1],cnt=1; int cur[maxn<<1]; struct edge{ int to,next,w; }G[4000010]; int vis[maxn<<1]; queue<int>Q; void add(int x,int y,int z){ ++cnt; G[cnt].to=y;G[cnt].next=h[x];G[cnt].w=z;h[x]=cnt; ++cnt; G[cnt].to=x;G[cnt].next=h[y];G[cnt].w=0;h[y]=cnt; } void read(int &num){ num=0;char ch=getchar(); while(ch<‘!‘)ch=getchar(); while(ch>=‘0‘&&ch<=‘9‘)num=num*10+ch-‘0‘,ch=getchar(); } bool BFS(){ memset(vis,-1,sizeof(vis)); Q.push(S);vis[S]=1; while(!Q.empty()){ int u=Q.front();Q.pop(); for(int i=h[u];i;i=G[i].next){ int v=G[i].to; if(vis[v]==-1&&G[i].w>0){ vis[v]=vis[u]+1; Q.push(v); } } }return vis[T]!=-1; } int DFS(int x,int f){ if(x==T||f==0)return f; int w,used=0; for(int i=cur[x];i;i=G[i].next){ if(vis[G[i].to]==vis[x]+1){ w=f-used; w=DFS(G[i].to,min(w,G[i].w)); G[i].w-=w;G[i^1].w+=w; if(G[i].w>0)cur[x]=i; used+=w;if(used==f)return used; } } if(!used)vis[x]=-1; return used; } void dinic(){ while(BFS()){ for(int i=S;i<=T;++i)cur[i]=h[i]; DFS(S,oo); }return; } int main(){ read(n);S=0;T=400001; for(int i=1;i<=n;++i){ read(X[i]);read(Y[i]); NX[X[i]]++;NY[Y[i]]++; } for(int i=1;i<=n;++i)add(h(X[i]),l(Y[i]),1); for(int i=1;i<=n;++i){ if(!visx[X[i]]){ add(S,h(X[i]),NX[X[i]]>>1); visx[X[i]]=true; } if(!visy[Y[i]]){ add(l(Y[i]),T,NY[Y[i]]>>1); visy[Y[i]]=true; } } dinic(); memset(visx,false,sizeof(visx)); memset(visy,false,sizeof(visy)); for(int i=1;i<=n;++i){ if(!visx[X[i]]){ visx[X[i]]=true; if(NX[X[i]]&1)add(S,h(X[i]),1); } if(!visy[Y[i]]){ visy[Y[i]]=true; if(NY[Y[i]]&1)add(l(Y[i]),T,1); } } dinic(); for(int i=1;i<=n;++i){ if(!G[(i<<1)+1].w)printf("b"); else printf("r"); }return 0; }
codeforces #305 D Mike and Fish
标签:
原文地址:http://www.cnblogs.com/joyouth/p/5352888.html