wa20分求调/kel
查看原帖
wa20分求调/kel
378346
expnoi楼主2023/8/4 09:55
#include<bits/stdc++.h>
#define int long long
using namespace std;
inline int read()
{
	int s=0,w=1;
	char c=getchar();
	while(c<'0'||c>'9')
	{
		if(c=='-')w=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')s=(s<<3)+(s<<1)+(c^48),c=getchar();
	return s*w;
}
inline void print(int x)
{
	if(x<0)x=-x,putchar('-');
	if(x>=10)print(x/10);
	putchar(x%10+48);
}
int n,g[1010][1010];
struct node
{
	int v,w,c,next;
}e[1000010];
int eid=0,head[1000010],pre[1000010],lst[1000010],dis[1000010],flow[1000010],out[1000010],vis[1000010];
inline void insert(int u,int v,int w,int c)
{
	e[eid].v=v;
	e[eid].w=w;
	e[eid].c=c;
	e[eid].next=head[u];
	head[u]=eid++;
}
int tot=0,s,t,ans;
map<pair<int,int>,int> mp,to;
queue<int> q;
const int inf=1e15;
inline bool spfa()
{
	for(int i=1;i<=tot;i++)dis[i]=flow[i]=inf,pre[i]=lst[i]=vis[i]=0;
	dis[s]=0;
	while(q.size())q.pop();
	q.push(s);
	vis[s]=1;
	while(q.size())
	{
		int u=q.front();
		q.pop();
		vis[u]=0;
		for(int i=head[u];~i;i=e[i].next)
		{
			int v=e[i].v,w=e[i].w;
			if(!w)continue;
			if(dis[u]+e[i].c<dis[v])
			{
				dis[v]=dis[u]+e[i].c;
				pre[v]=u;
				lst[v]=i;
				flow[v]=min(flow[u],e[i].w);
				if(!vis[v])
				{
					vis[v]=1;
					q.push(v);
				}
			}
		}
	}
	return pre[t];
}
signed main()
{
	memset(head,-1,sizeof(head));
	n=read();
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
		{
			g[i][j]=read();
			if(g[i][j]==1)out[i]++;
		}
	tot=n;
	s=++tot;
	t=++tot;
	for(int i=1;i<=n;i++)
	{
		for(int j=i+1;j<=n;j++)//考虑对每条不确定的边建图 
		{
			if(g[i][j]<2)continue;
			++tot;
			insert(s,tot,1,0);
			insert(tot,s,0,0);
			insert(tot,i,1,0);
			to[{i,j}]=eid-1;
			insert(i,tot,0,0);
			insert(tot,j,1,0);
			to[{j,i}]=eid-1;
			insert(j,tot,0,0);
		}
	}
	for(int i=1;i<=n;i++)
	{
		ans+=out[i]*(out[i]-1)/2;
		for(int j=out[i];j<=n;j++)
		{
			insert(i,t,1,j);
			insert(t,i,0,-j);
		}
	}
	int ans=0;
	while(spfa())
	{
		int now=t;
		ans+=dis[t]*flow[t];
		while(now!=s)
		{
			e[lst[now]].w-=flow[t];
			e[lst[now]^1].w+=flow[t];
			now=pre[now];
		}
	}
	print(n*(n-1)*(n-2)/6-ans);
	puts("");
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			if(g[i][j]<2)continue;
			if(e[to[{i,j}]].w==0)g[i][j]=1;
			else g[i][j]=0;
		}
	}
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
			print(g[i][j]);
			putchar(' ');
		}
		puts("");
	}
}
2023/8/4 09:55
加载中...