EK费用流居然写错了……

Description

    想必大家都看过成龙大哥的《80天环游世界》,里面的紧张刺激的打斗场面一定给你留下了深刻的印象。现在就有这么
    一个80人的团伙,也想来一次环游世界。
    他们打算兵分多路,游遍每一个国家。
    因为他们主要分布在东方,所以他们只朝西方进军。设从东方到西方的每一个国家的编号依次为1...N。假若第i个人的游历路线为P1、P2......Pk(0≤k≤N),则P1<P2<......<Pk。
    众所周知,中国相当美丽,这样在环游世界时就有很多人经过中国。我们用一个正整数Vi来描述一个国家的吸引程度,Vi值越大表示该国家越有吸引力,同时也表示有且仅
有Vi个人会经过那一个国家。
    为了节省时间,他们打算通过坐飞机来完成环游世界的任务。同时为了省钱,他们希望总的机票费最小。
    明天就要出发了,可是有些人临阵脱逃,最终只剩下了M个人去环游世界。他们想知道最少的总费用,你能告诉他们吗?

Input

第一行两个正整数N,M。
第二行有N个不大于M正整数,分别表示V1,V2......VN。
接下来有N-1行。第i行有N-i个整数,该行的第j个数表示从第i个国家到第i+j个国家的机票费(如果该值等于-1则表示这两个国家间没有通航)。

Output

在第一行输出最少的总费用。

Sample Input

6 3
2 1 3 1 2 1
2 6 8 5 0
8 2 4 1
6 1 0
4 -1
4

Sample Output

27

HINT

1<= N < =100 1<= M <= 79


题目分析

这就是一个无源汇最小费用最大流的模型,关于M个人的限制则是把(TT,SS)的边流量设为M

掌握得还是不够熟练啊

 #include<bits/stdc++.h>
const int maxn = ;
const int maxm = ;
const int INF = 2e9; struct Edge
{
int u,v,f,c,cst;
Edge(int a=, int b=, int c=, int d=, int e=):u(a),v(b),f(c),c(d),cst(e) {}
}edges[maxm];
int n,m,S,T,SS,TT,ans;
int edgeTot,head[maxn],nxt[maxm],bck[maxm],cst[maxn],flw[maxn];
bool inq[maxn<<]; void addedge(int u, int v, int c, int cst)
{
edges[edgeTot] = Edge(u, v, , c, cst), nxt[edgeTot] = head[u], head[u] = edgeTot, ++edgeTot;
edges[edgeTot] = Edge(v, u, , , -cst), nxt[edgeTot] = head[v], head[v] = edgeTot, ++edgeTot;
}
void maxFlow()
{
for (;;)
{
std::queue<int> q;
memset(flw, , sizeof flw);
memset(bck, , sizeof bck);
memset(cst, 0x3f3f3f3f, sizeof cst);
q.push(S), flw[S] = INF, cst[S] = ;
for (int tmp; q.size(); )
{
tmp = q.front(), q.pop(), inq[tmp] = ;
for (int i=head[tmp]; i!=-; i=nxt[i])
{
int v = edges[i].v;
if (cst[tmp]+edges[i].cst < cst[v]&&edges[i].f < edges[i].c){
cst[v] = cst[tmp]+edges[i].cst, bck[v] = i;
flw[v] = std::min(flw[tmp], edges[i].c-edges[i].f);
if (!inq[v]) q.push(v), inq[v] = ;
}
}
}
if (!flw[T]) break;
for (int i=T; i!=S; i=edges[bck[i]].u)
edges[bck[i]].f += flw[T], edges[bck[i]^].f -= flw[T];
ans += flw[T]*cst[T];
}
}
int main()
{
memset(head, -, sizeof head);
scanf("%d%d",&n,&m);
S = n+n+, T = S+, SS = T+, TT = SS+;
for (int i=,v; i<=n; i++)
{
scanf("%d",&v);
addedge(S, i+n, v, );
addedge(i, T, v, );
addedge(SS, i, INF, );
addedge(i+n, TT, INF, );
}
addedge(TT, SS, m, );
for (int i=; i<=n; i++)
for (int j=,x; j<=n-i; j++)
{
scanf("%d",&x);
if (x!=-) addedge(i+n, i+j, INF, x);
}
maxFlow();
printf("%d\n",ans);
return ;
}

END

最新文章

  1. 【转】使用SQL Tuning Advisor STA优化SQL
  2. CYQ.Data V5 分布式缓存MemCached应用开发介绍
  3. Android源码
  4. IIS URL Rewrite redirect from one Domain to another
  5. 利用Python获取ZOJ所有题目的名字
  6. springmvc学习(三)
  7. javascript代码放置位置对程序的影响
  8. bzoj 1191
  9. python数据类型以及模块的含义
  10. 初学VUE2.0
  11. 【sping揭秘】19、关于spring中jdbctemplate中的DataSource怎么来呢
  12. SAP FI CO模块常用事务代码
  13. Java中ArrayList和LinkedList区别 时间复杂度 与空间复杂度
  14. IOS 获取中英文字符串长度
  15. MySQL修改编码设置及乱码问题
  16. C语言文本处理
  17. 关于easyui表格右侧多出来的那一列。
  18. [OpenCV Qt教程] 在Qt图形界面中显示OpenCV图像的OpenGL Widget (第一部分)
  19. 对fgets的理解
  20. 通过 zxing 生成二维码

热门文章

  1. MyBatis二级缓存配置
  2. Java-GC-标记清除算法
  3. python——字符编码
  4. idea 添加yuicompressor压缩js/css
  5. 从I/O事件到阻塞、非阻塞、poll到epoll的理解过程
  6. openstack安装newton版本Glance部署(二)
  7. vue2.0:(五)、路由vue-router
  8. Python一个有意思的地方:reduce、map、filter
  9. LR脚本示例之常用函数
  10. 通过Jenkins调用自动部署war包及jar包到服务器上的Shell脚本