sap模板 别人给的
来源:互联网 发布:万科 股权 知乎 编辑:程序博客网 时间:2024/04/28 19:03
#include <cstdio>
#include <cstring>
#pragma comment(linker, "/STACK:1024000000,1024000000") //栈溢出时,强行开栈
const int MAXN=80000;
const int inf = 0x7fffffff;
const int s = 0;
struct edge{
int v,next,w;
}edge[1000000];
int head[2*MAXN],cnt;//for sap
void addedge(int u, int v, int w)
{
edge[cnt].v = v;
edge[cnt].w = w;
edge[cnt].next = head[u];
head[u] = cnt++;
edge[cnt].v = u;
edge[cnt].w = 0;
edge[cnt].next = head[v];
head[v] = cnt++;
}
int sap(int t)
{
int pre[2*MAXN],cur[2*MAXN],dis[2*MAXN],gap[2*MAXN]; //gap标记断层 这里可能会暴 如果点的范围非常大的话 就放在外面 每次清空一下。
int flow = 0 , aug = inf ,u;
bool flag;
for (int i = 0 ; i <= t ; ++i)
{
cur[i] = head[i];
gap[i] = dis[i] = 0;
}
gap[s] = t+1;
u = pre[s] = s;
while (dis[s] <= t)
{
flag = 0 ;
for (int &j = cur[u] ; ~j ; j = edge[j].next)
{
int v = edge[j].v;
if (edge[j].w > 0 && dis[u] == dis[v]+1)
{
flag = 1;
if( edge[j].w < aug )aug = edge[j].w;
pre[v] = u;
u = v;
if (u == t)
{
flow += aug;
while (u != s)
{
u = pre[u];
edge[cur[u]].w -= aug;
edge[cur[u]^1].w += aug;
}
aug = inf;
}
break;
}
}
if ( flag )continue ;
int mindis = t+1;
for (int j = head[u]; ~j ; j = edge[j].next)
{
int v = edge[j].v;
if (edge[j].w > 0 && dis[v] < mindis)
{
mindis = dis[v];
cur[u] = j;
}
}
if(--gap[dis[u]] == 0)break;
gap[ dis[u] = mindis+1 ]++;
u = pre[u];
}
return flow;
}
void init ()
{
memset (head , -1 , sizeof(head));
cnt=0;
}
- sap模板 别人给的
- 给别人的回帖
- 别人给的
- 别人给的读书计划!
- 别人给我的建议
- 别人给我的概念
- 大数模板(别人的)
- 如何查看别人的模板?
- 给SAP初学者的建议!
- 周杰伦给别人写的歌
- 把自己的热情传递给别人
- 不妨扯下自己的面子给别人
- “给别人加水”的职场智慧
- 开放自己的Apache给别人访问
- 伤口是不用给别人看的,
- 把自己的库共享给别人
- 最好的辞职模板(别人的)
- sap模板
- ABAP 程序完成EXCEL文件的导入导出
- 总是在告诉自己
- 显示文件扩展名
- 科学家发现超光速粒子 或颠覆爱因斯坦狭义相对论
- windows proxy 注册表
- sap模板 别人给的
- PKU ACM 1009 Edge Detection 边缘提取算法
- Java线程:线程的调度-优先级
- CString成员函数用法大全
- android 任务和Activity详解
- Win2003、XP上修改Oracle链接数,并重启数据库
- 消息服务器ActiveMQ-5.2
- C++中对hash_map自定义哈希函数和比较函数的理解
- SSH下直接向主机传送文件