5331. 【NOIP2017提高A组模拟8.23】壕游戏 费用流 动态连边
2017-08-23 14:51
253 查看
题意:n个点m条边,每条边有一个经过次数限制ci,假设经过了x次,第x次费用为ai*x+bi,有s个终点,要跑k次,问最小总费用多少。
巨鬼畜,连边只按照第一次的连,然后费用流的时候每流一条边然后再连下一条边。。
我好菜啊,ymw随手AK。
巨鬼畜,连边只按照第一次的连,然后费用流的时候每流一条边然后再连下一条边。。
我好菜啊,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); }
相关文章推荐
- jzoj5331 【NOIP2017提高A组模拟8.23】壕游戏
- 【jzoj5347】【NOIP2017提高A组模拟9.5】【遥远的金字塔】【斜率优化动态规划】
- 【jzoj5332】【NOIP2017提高A组模拟8.23】【密码】【ac自动机】【动态规划】
- 5332. 【NOIP2017提高A组模拟8.23】密码 AC自动机+数位DP
- 5333. 【NOIP2017提高A组模拟8.23】大新闻
- 【jzoj5333】【NOIP2017提高A组模拟8.23】【大新闻】【可持久化线段树】
- 【JZOJ5332】【NOIP2017提高A组模拟8.23】密码
- 【JZOJ 5395】【NOIP2017提高A组模拟10.6】Count
- JZOJ 5185. 【NOIP2017提高组模拟6.30】tty's sequence
- JZOJ5398. 【NOIP2017提高A组模拟10.7】Adore
- 【jzoj5360】【NOIP2017提高A组模拟9.12】【Shorten Diameter】
- 【jzoj5289】【NOIP2017提高组A组模拟8.17】【偷笑】【数据结构】
- 【JZOJ 5405】【NOIP2017提高A组模拟10.10】Permutation
- [jzoj5307]【NOIP2017提高A组模拟8.18】偷窃
- 【JZOJ 5377】【NOIP2017提高A组模拟9.19】开拓
- 【JZOJ5287】【NOIP2017提高组模拟】最短路
- [jzoj100047]【NOIP2017提高A组模拟7.14】基因变异
- 【JZOJ 5248】【NOIP2017提高A组模拟8.10】花花的聚会
- jzoj100048 【NOIP2017提高A组模拟7.14】紧急撤离 (网格图,分治,bitSet)
- JZOJ5384. 【NOIP2017提高A组模拟9.23】四维世界 组合数学