网络流

关键是建图,思路在代码里

/*
最大流SAP
邻接表
思路:基本源于FF方法,给每个顶点设定层次标号,和允许弧。
优化:
1、当前弧优化(重要)。
1、每找到以条增广路回退到断点(常数优化)。
2、层次出现断层,无法得到新流(重要)。
时间复杂度(m*n^2)
*/
#include <iostream>
#include <cstdio>
#include <cstring>
#define ms(a,b) memset(a,b,sizeof a)
using namespace std;
const int INF = ;
int G[INF][INF];
struct node {
int v, c, next;
} edge[INF*INF*];
int pHead[INF*INF], SS, ST, nCnt;
void addEdge (int u, int v, int c) {
edge[++nCnt].v = v; edge[nCnt].c = c, edge[nCnt].next = pHead[u]; pHead[u] = nCnt;
edge[++nCnt].v = u; edge[nCnt].c = , edge[nCnt].next = pHead[v]; pHead[v] = nCnt;
}
int SAP (int pStart, int pEnd, int N) {
int numh[INF], h[INF], curEdge[INF], pre[INF];
int cur_flow, flow_ans = , u, neck, i, tmp;
ms (h, ); ms (numh, ); ms (pre, -);
for (i = ; i <= N; i++) curEdge[i] = pHead[i];
numh[] = N;
u = pStart;
while (h[pStart] <= N) {
if (u == pEnd) {
cur_flow = 1e9;
for (i = pStart; i != pEnd; i = edge[curEdge[i]].v)
if (cur_flow > edge[curEdge[i]].c) neck = i, cur_flow = edge[curEdge[i]].c;
for (i = pStart; i != pEnd; i = edge[curEdge[i]].v) {
tmp = curEdge[i];
edge[tmp].c -= cur_flow, edge[tmp ^ ].c += cur_flow;
}
flow_ans += cur_flow;
u = neck;
}
for ( i = curEdge[u]; i != ; i = edge[i].next)
if (edge[i].c && h[u] == h[edge[i].v] + ) break;
if (i != ) {
curEdge[u] = i, pre[edge[i].v] = u;
u = edge[i].v;
}
else {
if ( == --numh[h[u]]) continue;
curEdge[u] = pHead[u];
for (tmp = N, i = pHead[u]; i != ; i = edge[i].next)
if (edge[i].c) tmp = min (tmp, h[edge[i].v]);
h[u] = tmp + ;
++numh[h[u]];
if (u != pStart) u = pre[u];
}
}
return flow_ans;
}
/*
poj1149 最大流
建图:
每个顾客为一个节点,从源点到每个猪圈的第一个顾客连接一条容量为猪圈猪数目的边
从第一个顾客到同一个猪圈的其它顾客连一条容量无限的边
每个顾客到汇点的边的容量为他最大买的数量 */
int k, c, m, n,x,y,z;
int vis[INF],sum[INF];
void solve(){
nCnt=;
for(int i=;i<=ST;i++)
for(int j=;j<=ST;j++){
if(G[i][j]) addEdge(i,j,G[i][j]);
}
int ans=SAP(SS,ST,ST);
printf("%d\n",ans);
}
int main() {
/*
先用邻接矩阵统计容量,再用前向星存边,表头在pHead[],初始化nCnt=1
SS,ST分别为源点和汇点
*/
scanf("%d %d",&m,&n);
SS=n+,ST=n+;
for(int i=;i<=m;i++)
scanf("%d",&sum[i]);
for(int i=;i<=n;i++){
scanf("%d",&k);
for(int j=;j<=k;j++){
scanf("%d",&x);
if(!vis[x]){
G[SS][i]+=sum[x];
vis[x]=i;
}
else{
G[vis[x]][i]=0xffff;
}
}
scanf("%d",&k);
G[i][ST]=k;
}
solve();
return ;
}

最新文章

  1. WIFI网络操作
  2. 利用Nginx实现域名转发 不修改主机头
  3. django个人总结
  4. .NET Framework 的 Quirk Version
  5. c++ static变量
  6. River Crossing 简单的动态规划 .
  7. javascript中ajax post实例详解
  8. margin折叠
  9. web页面的生命周期
  10. vim setting
  11. 解决NoSuchMethodError with Spring MutableValues异常问题
  12. C/C++ 基础教程
  13. 定位vc运行时库问题 依赖问题,屡试不爽的一招
  14. sql查询优化整理
  15. memcached配置
  16. python 之生成器的介绍
  17. ESXI6时间源快速同步
  18. PAT1135(红黑书的判定)
  19. C#编程(小结)---------- 小总结
  20. PAT-GPLT L1-033 - 出生年 - [简单模拟]

热门文章

  1. TimePicker,TimePickerDialog以及自定义timepicker(一)
  2. HTML5 Canvas JavaScript库 Fabric.js 使用经验
  3. UVAlive 3263 That Nice Euler Circuit(欧拉定理)
  4. bzoj 1188 [HNOI2007]分裂游戏(SG函数,博弈)
  5. vijosP1195“非常男女”计划
  6. 最受IT公司欢迎的50款开源软件
  7. Div 3 - SGU 105(找规律)
  8. ASP.NET与SOAP协议使用记录
  9. 【设计模式 - 18】之备忘录模式(Memento)
  10. java基础之 IO流