题意
这题主要要处理两个限制：
1.最多有\(m\)个人，即最多\(m\)条路径。
2.每个点\(i\)的经过次数是给定的，为\(a_i\)。
其中2.可以用限定上下界的方法解决，1.在建图中讲。
建图：
新建三个点：\(S'\),\(T'\),\(S''\)，将每个点\(i\)拆成入点\(i\)和出点\(i'\)。
1.从\(S''\)向每个入点连容量为\([0,inf]\)费用为\(0\)的边。
2.从每个出点向\(T'\)连容量为\([0,inf]\)费用为\(0\)的边。
3.从\(S'\)向\(S''\)连容量为\([0,m]\)费用为\(0\)的边,这样就解决了限制1。
4.从每个点的入点向出点连容量为\([a_i,a_i]\)费用为\(0\)的边。
5.对于能通行的点对\((i,j)(i<j)\)，从\(i\)的出点向\(j\)的入点连边。
之后我们跑\(S'->T'\)的有源汇上下界最小费用流即为答案。
code:
#include<bits/stdc++.h>
using namespace std;
const int maxn=210;
const int inf=1e9;
int n,m,cnt_edge=1,tot,S,T;
int head[maxn],dis[maxn],a[maxn],in[maxn],out[maxn];
int cost[maxn][maxn];
bool vis[maxn];
struct Edge{int u,v,down,up,cost;}E[maxn*maxn+maxn*2];
struct edge{int to,nxt,flow,cost;}e[(maxn*maxn+maxn*4)<<1];
inline int read()
{
char c=getchar();int res=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9')res=res*10+c-'0',c=getchar();
return res*f;
}
inline void add(int u,int v,int w,int c)
{
e[++cnt_edge].nxt=head[u];
head[u]=cnt_edge;
e[cnt_edge].to=v;
e[cnt_edge].flow=w;
e[cnt_edge].cost=c;
}
inline void addflow(int u,int v,int w,int c){add(u,v,w,c);add(v,u,0,-c);}
inline bool spfa()
{
memset(dis,0x3f,sizeof(dis));
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();vis[x]=0;
for(int i=head[x];i;i=e[i].nxt)
{
int y=e[i].to;
if(dis[y]>dis[x]+e[i].cost&&e[i].flow>0)
{
dis[y]=dis[x]+e[i].cost;
if(!vis[y])q.push(y),vis[y]=1;
}
}
}
return dis[T]!=0x3f3f3f3f;
}
int dfs(int x,int lim)
{
vis[x]=1;
if(lim<=0||x==T)return lim;
int res=lim;
for(int i=head[x];i;i=e[i].nxt)
{
int y=e[i].to;
if(dis[y]!=dis[x]+e[i].cost||e[i].flow<=0||vis[y])continue;
int tmp=dfs(y,min(res,e[i].flow));
res-=tmp;
e[i].flow-=tmp,e[i^1].flow+=tmp;
if(res<=0)break;
}
return lim-res;
}
inline int Dinic(int st,int ed)
{
S=st,T=ed;
int res=0,cost=0;
while(spfa())
{
int flow=dfs(S,inf);
res+=flow;cost+=flow*dis[T];
}
return cost;
}
int main()
{
//freopen("test.in","r",stdin);
//freopen("test.out","w",stdout);
n=read(),m=read();
for(int i=1;i<=n;i++)a[i]=read();
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
cost[i][j]=read();
E[++tot]=(Edge){2*n+1,2*n+2,0,m,0};//S'=2*n+1,S''=2*n+2,T'=2*n+3.
for(int i=1;i<=n;i++)E[++tot]=(Edge){2*n+2,i,0,inf,0},E[++tot]=(Edge){i+n,2*n+3,0,inf,0};
for(int i=1;i<=n;i++)E[++tot]=(Edge){i,i+n,a[i],a[i]};
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
if(~cost[i][j])E[++tot]=(Edge){i+n,j,0,inf,cost[i][j]};
E[++tot]=(Edge){2*n+3,2*n+1,0,inf,0};
for(int i=1;i<=tot;i++)addflow(E[i].u,E[i].v,E[i].up-E[i].down,E[i].cost);
for(int i=1;i<=tot;i++)in[E[i].v]+=E[i].down,out[E[i].u]+=E[i].down;
for(int i=1;i<=2*n+3;i++)//st=2*n+4,ed=2*n+5.
if(in[i]>=out[i])addflow(2*n+4,i,in[i]-out[i],0);
else addflow(i,2*n+5,out[i]-in[i],0);
printf("%d",Dinic(2*n+4,2*n+5));
return 0;
}