首先转化题目： 给你 n 个物品，可以选或不选。 选第i个物品需要 c[i] 的代价。 同时选第 i 和第j个物品获得 b[i][j] 的收益 问最大收益。 网络流建图： 考虑 sum=∑ni=1∑nj=1b[i][j] 剩下就变成了算代价最小（也就是最小割） 建点 (i,j) 和 i ，对于每个点(i,j)，向 i,j 连一条容量为 inf 的边。 S 向每个点(i,j)连一条容量为 b[i][j] 的边 每个点 i 向T连一条容量为 c[i] 的边。 然后跑最大流就好。
#include <bits/stdc++.h> #define gc getchar() #define ll long long #define N 509 #define inf 0x3f3f3f3f using namespace std; int n,b[N][N],c[N],sum; int number=1,pos,cur[N*N],dis[N*N]; bool vis[N*N]; vector<int> G[N*N]; struct edge { int from,to,flow,cap; int rest() { return cap-flow; } void add(int x,int y,int z) { from=x,to=y,cap=z,flow=0; } }e[N*N<<3]; void add(int x,int y,int z) { e[++number].add(x,y,z); G[x].push_back(number); e[++number].add(y,x,0); G[y].push_back(number); } #define E e[G[x][i]] bool bfs(int s,int t) { memset(vis,0,sizeof(vis)); queue<int> Q; Q.push(s); dis[s]=0; vis[s]=1; while (!Q.empty()) { int x=Q.front(); Q.pop(); for (int i=0;i<(int)G[x].size();i++) if (!vis[E.to]&&E.rest()>0) { vis[E.to]=1; dis[E.to]=dis[x]+1; Q.push(E.to); } } return vis[t]; } int dfs(int x,int a,int t) { if (x==t||a==0) return a; int flow=0,f; for (int &i=cur[x];i<(int)G[x].size();i++) if (dis[x]+1==dis[E.to]&&(f=dfs(E.to,min(a,E.rest()),t))>0) { E.flow+=f; e[G[x][i]^1].flow-=f; flow+=f; a-=f; if (!a) break; } return flow; } int Maxflow(int s,int t) { int flow=0; while (bfs(s,t)) { memset(cur,0,sizeof(cur)); flow+=dfs(s,inf,t); } return flow; } #undef E int read() { int x=1; char ch; while (ch=gc,ch<'0'||ch>'9') if (ch=='-') x=-1; int s=ch-'0'; while (ch=gc,ch>='0'&&ch<='9') s=s*10+ch-'0'; return s*x; } int main() { n=read(); for (int i=1;i<=n;i++) for (int j=1;j<=n;j++) { b[i][j]=read(),sum+=b[i][j]; add(0,(i-1)*n+j,b[i][j]); add((i-1)*n+j,n*n+i,inf); add((i-1)*n+j,n*n+j,inf); } for (int i=1;i<=n;i++) { c[i]=read(); add(n*n+i,n*n+n+1,c[i]); } printf("%d\n",sum-Maxflow(0,n*n+n+1)); return 0; }