Cod sursa(job #4118)

Utilizator gigi_becaliGigi Becali gigi_becali Data 30 decembrie 2006 18:46:20
Problema Taramul Nicaieri Scor 100
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.46 kb
#include <cstdio>
#include <string>
#define maxn 256

int main()
{
	int c[maxn][maxn], n, source, sink, t[maxn], Q[maxn], viz[maxn], p, q, u;
	int i, j, nr=0, ok=1;
	// citire
	freopen("harta.in", "r", stdin);
	scanf("%d\n", &n);

	sink=2*n+1;
	source=0;
	
	for(i=1;i<=n;i++)
	{
		scanf("%d %d\n", &p, &q);
		
		c[source][i]=p;
		c[i+n][sink]=q;
	}

	for(i=1;i<=n;i++)
		for(j=n+1;j<=n<<1;j++) c[i][j]=1;
		
	for(i=1;i<=n;i++) c[i][i+n]=0;
	//end citire

	//flow
	 nr=0, ok=1;
	
	while(1)
	{
	
		 p=1, q=1;
		memset(t, -1, sizeof(t));
		memset(viz, 0, sizeof(viz));
		Q[1]=source;
		viz[source]=1;
	
		while(p<=q)
		{
			u=Q[p++];
		
			for(i=source;i<=sink;i++)
				if(c[u][i] && !viz[i])
				{
					viz[i]=1;
					t[i]=u;
					Q[++q]=i;
					if(i==sink) goto exit_loop;
			}
	}
	exit_loop:
	

	if(t[sink]==-1) break;
	
			
		for(i=source;i<sink;i++)
		{
			if(t[i]==-1 || c[i][sink]<=0)continue;
			
			ok=1;
			for(j=i;j!=source;j=t[j]) if(c[t[j]][j]==0){ok=0; break;}
			if(ok==0) continue;
	
			c[i][sink]-=1;
			c[sink][i]+=1;
			for(j=i; j!=source; j=t[j])
			{
				c[t[j]][j]-=1;
				c[j][t[j]]+=1;
			}
		}
	}
	//end flow
	
	// afis
	for(i=1;i<=n;i++)
		for(j=n+1;j<=n<<1;j++) if(!c[i][j] && i!=j-n) nr++;
	
	freopen("harta.out", "w", stdout);
	printf("%d\n", nr);
	
	for(i=1;i<=n;i++)
		for(j=n+1;j<=n<<1;j++) if(!c[i][j] && i!=j-n) printf("%d %d\n", i, j-n);
		
	//end afis

	return 0;
}