码迷,mamicode.com
首页 > 其他好文 > 详细

codeforces #305 D Mike and Fish

时间:2016-04-04 22:25:09      阅读:127      评论:0      收藏:0      [点我收藏+]

标签:

正解貌似是大暴搜?

首先我们考虑这是一个二分图,建立网络流模型后很容易得出一个算法

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

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!