您的位置:首页 > 其它

5331. 【NOIP2017提高A组模拟8.23】壕游戏 费用流 动态连边

2017-08-23 14:51 253 查看
题意:n个点m条边,每条边有一个经过次数限制ci,假设经过了x次,第x次费用为ai*x+bi,有s个终点,要跑k次,问最小总费用多少。

巨鬼畜,连边只按照第一次的连,然后费用流的时候每流一条边然后再连下一条边。。

我好菜啊,ymw随手AK。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define T n+1
using namespace std;
const int N=1e7+5;
const int inf=0x3f3f3f3f;
int n,m,k,s;
int head
,go
,next
,c
,num
,w
,A
,B
,C
;
int pre
,vis
,ans;
int dis[100005],tot=1;
queue<int> q;
inline void add(int x,int y,int z,int W,int a,int b,int cost,int Num)
{
go[++tot]=y;
next[tot]=head[x];
c[tot]=z,w[tot]=W,A[tot]=a,B[tot]=b,C[tot]=cost;
num[tot]=Num;
head[x]=tot;
}
inline void ins(int x,int y,int z,int w,int a,int b,int c,int num)
{
add(x,y,z,w,a,b,c,num);
add(y,x,0,-w,a,b,c,num);
}
inline bool spfa()
{
fo(i,0,T)dis[i]=inf;
q.push(0);
dis[0]=0,vis[0]=1;
while (!q.empty())
{
int x=q.front();
q.pop();
for(int i=head[x];i;i=next[i])
{
int v=go[i];
if (c[i]&&dis[x]+w[i]<dis[v])
{
dis[v]=dis[x]+w[i];
pre[v]=i;
if (!vis[v])
{
vis[v]=1;
q.push(v);
}
}
}
vis[x]=0;
}
if (dis[T]==inf)return 0;
else return 1;
}
inline void mcf()
{
ans+=dis[T];
int x=T;
while (x!=0)
{
int a=A[pre[x]],b=B[pre[x]],cost=C[pre[x]],W=w[pre[x]];
if (cost&&num[pre[x]]<cost)
ins(go[pre[x]^1],go[pre[x]],1,W+a,a,b,cost,num[pre[x]]+1);
C[pre[x]]=C[pre[x]^1]=0;
c[pre[x]]--;
c[pre[x]^1]++;
x=go[pre[x]^1];
}
}
int main()
{
freopen("game.in","r",stdin);
freopen("game.out","w",stdout);
scanf("%d%d%d%d",&n,&m,&k,&s);
ins(0,1,inf,0,0,0,0,0);
fo(i,1,s)
{
int x;
scanf("%d",&x);
ins(x,T,inf,0,0,0,0,0);
}
fo(i,1,m)
{
int x,y,A,B,C;
scanf("%d%d%d%d%d",&x,&y,&A,&B,&C);
ins(x,y,1,A+B,A,B,C,1);
}
fo(i,1,k)
{
if (!spfa())
{
puts("-1\n");
return 0;
}
mcf();
}
printf("%d\n",ans);
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: