「TJOI2015」线性代数 解题报告

「TJOI2015」线性代数

和牛客某题很像

在和里面有\(B_{i,j}\)要求是\(A_i,A_j\)都为\(1\)，和里面减去\(C_i\)要求\(A_i\)为\(1\)，然后先把贡献也就是\(B_{i,j}\)加上，然后建图跑最小割，建图方法和牛客那题一样

另外1e6条边为什么也可以跑过去啊，还很快..

---

Code:

#include <cstdio>

#include <cctype>

#include <cstring>

const int N=510;

const int M=1002000;

int min(int x,int y){return x<y?x:y;}

template <class T>

void read(T &x)

{

	x=0;char c=getchar();

	while(!isdigit(c)) c=getchar();

	while(isdigit(c)) x=x*10+c-'0',c=getchar();

}

int head[N],to[M],Next[M],edge[M],cnt=1;

void add(int u,int v,int w)

{

	to[++cnt]=v,edge[cnt]=w,Next[cnt]=head[u],head[u]=cnt;

	to[++cnt]=u,edge[cnt]=0,Next[cnt]=head[v],head[v]=cnt;

}

int q[N],l,r,dep[N],s,t,A[N];

bool bfs()

{

	memset(dep,0,sizeof dep);

	dep[q[l=r=1]=s]=1;

	while(l<=r)

	{

		int now=q[l++];

		for(int v,i=head[now];i;i=Next[i])

			if(edge[i]&&!dep[v=to[i]])

			{

				dep[v]=dep[now]+1;

				if((q[++r]=v)==t) return true;

			}

	}

	return false;

}

const int inf=0x3f3f3f3f;

int dfs(int now,int flow)

{

	if(now==t) return flow;

	int res=flow,yuu;

	for(int v,i=head[now];i&&res;i=Next[i])

		if(edge[i]&&dep[v=to[i]]==dep[now]+1)

		{

			yuu=dfs(v,min(res,edge[i]));

			if(!yuu){dep[v]=0;continue;}

			edge[i]-=yuu,edge[i^1]+=yuu;

			res-=yuu;

		}

	return flow-res;

}

int main()

{

	int ans=0,n;

	read(n);

	s=n+1,t=s+1;

	for(int a,i=1;i<=n;i++)

		for(int j=1;j<=n;j++)

		{

			read(a);

			A[i]+=a;

			ans+=a;

			if(i!=j) add(j,i,a);

		}

	for(int a,i=1;i<=n;i++)

	{

		read(a);

		add(s,i,A[i]);

		add(i,t,a);

	}

	int flow;

	while(bfs()) while(flow=dfs(s,inf)) ans-=flow;

	printf("%d\n",ans);

	return 0;

}

---

2019.2.24