hdu3549Flow Problem poj1273Drainage Ditches poj3281Dining sgu185 Network flow

最新推荐文章于 2017-07-14 22:34:53 发布
最新推荐文章于 2017-07-14 22:34:53 发布
阅读量379 收藏
点赞数
分类专栏： Network flow 文章标签： Network flow
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/coder_coder_coder/article/details/38886171
版权

hdu3549Flow_Problem

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 5005;
const int M = 30005;
struct Enode{
    int y, c, next;
} e[M * 2];
//struct Point
//{
//	int son,cur,pre,lim,d;
//} a[maxn];
int n, m, tot, head[N], now[N], h[N], vh[N], augc, found, flow;
//h为距离,vh为距离标号的数目,augc当前边的流量,flow总流量
//int tot, n, m, st, ed, cnt[maxn];
void Addedge(int x, int y, int c){
    e[++tot].y = y; e[tot].c = c; e[tot].next = head[x]; head[x] = tot;
    e[++tot].y = x; e[tot].c = 0; e[tot].next = head[y]; head[y] = tot;
}
void Init(){
    int x, y, c, i;
    scanf("%d%d", &n, &m);
    tot = -1; memset(head, -1, sizeof(head));
    for (i = 0; i < m; i++){
        scanf("%d%d%d", &x, &y, &c);
        Addedge(x, y, c);
    }
    memcpy(now, head, sizeof(head));
}

void Aug(int x, int st, int ed, int n){//x当前点,st源点,ed汇点,n,点最大数目
    int p = now[x], minh = n - 1, augco = augc;
    if (x == ed){
        found = 1;
        flow += augc;
        return;
    }
    while (p != -1){//每条边
        if (e[p].c > 0 && h[e[p].y] + 1 == h[x]){//e[p].y下一点,可行弧
            augc = min(augc, e[p].c);
            Aug(e[p].y, st, ed, n);//下一点
            if (h[st] >= n) return;
            if (found) break;
            augc = augco;
        }
        p = e[p].next;//下一边
    }
    if (found){
        e[p].c -= augc;
        e[p ^ 1].c += augc;//边标号从0开始,反向边,p^1
    }else{
        p = head[x];
        while (p != -1){
            if (e[p].c > 0 && h[e[p].y] < minh){
                minh = h[e[p].y];//最小标号
                now[x] = p;
            }
            p = e[p].next;
        }
        vh[h[x]] --;//标号数目-1
        if (!vh[h[x]]) h[st] = n;//无符合标号
        h[x] = minh + 1;//修改当前节点标号
        vh[h[x]] ++;
    }
}
void Maxflow(int st, int ed, int n){
    flow = 0;
    memset(h, 0, sizeof(h));
    memset(vh, 0, sizeof(vh));
    vh[0] = n;//距离标号的数目
    while (h[st] < n){
        found = 0;
        augc = 1 << 30;//假设流量无穷大
        Aug(st, st, ed, n);
    }
}
int main(){
    int T,cas=0;
    scanf("%d",&T);
    while(T--){
        Init();
        Maxflow(1, n, n);
        printf("Case %d: %d\n",++cas, flow);
    }
    return 0;
}
poj1273Drainage Ditches 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define MAXN 2222
#define MAXM 444444 //邻接表要开边数的2倍
#define inf 1<<30

struct Edge {
    int v,cap,next;
} edge[MAXM];

int n,m,vs,vt,NE,NV;
int head[MAXN];

void Insert(int u,int v,int cap)
{
    edge[NE].v=v;
    edge[NE].cap=cap;
    edge[NE].next=head[u];
    head[u]=NE++;

    edge[NE].v=u;
    edge[NE].cap=0;
    edge[NE].next=head[v];
    head[v]=NE++;
}

int level[MAXN];//标记层次（距离标号）
int gap[MAXN];
//间隙优化，定义gap[i]为标号是i的点的个数
//在重标记i时，检查gap[level[i]],若减为0,这算法结束。

void bfs(int vt)
{
    memset(level,-1,sizeof(level));
    memset(gap,0,sizeof(gap));
    level[vt]=0;
    gap[level[vt]]++;
    queue<int>que;
    que.push(vt);
    while(!que.empty()) {
        int u=que.front();
        que.pop();
        for(int i=head[u]; i!=-1; i=edge[i].next) {
            int v=edge[i].v;
            if(level[v]!=-1)continue;
            level[v]=level[u]+1;
            gap[level[v]]++;
            que.push(v);

        }
    }
}


int pre[MAXN];//前驱
int cur[MAXN];

int SAP(int vs,int vt)//源点,汇点
{
    bfs(vt);//bfs优化,
    memset(pre,-1,sizeof(pre));
    memcpy(cur,head,sizeof(head));//cur保存的是当前弧
    int u=pre[vs]=vs,flow=0,aug=inf;//源点前驱还是它本身,aug表示增广路的可改尽量
    gap[0]=NV;
    while(level[vs]<NV) {
        bool flag=false;
        for(int &i=cur[u]; i!=-1; i=edge[i].next) {
            int v=edge[i].v;//v是u的后继
            //寻找可行弧
            if(edge[i].cap&&level[u]==level[v]+1) {
                flag=true;//找到可行弧
                pre[v]=u;//记录前驱
                u=v;
                //  aug=(aug==-1?edge[i].cap:min(aug,edge[i].cap));
                aug=min(aug,edge[i].cap);
                //如果找到 一条增广路
                if(v==vt) {
                    flow+=aug;//更新最大流
                    for(u=pre[v]; v!=vs; v=u,u=pre[u]) {
                        edge[cur[u]].cap-=aug;//前向弧容量减少
                        edge[cur[u]^1].cap+=aug;//后向弧容量增加
                    }
                    //     aug=-1;
                    aug=inf;
                }
                break;
            }
        }
        if(flag)continue;
        int minlevel=NV;
        //寻找与当前点相连接的点中最小的距离标号(重标号)
        for(int i=head[u]; i!=-1; i=edge[i].next) {
            int v=edge[i].v;
            if(edge[i].cap&&level[v]<minlevel) {
                minlevel=level[v];
                cur[u]=i;
            }
        }
        if(--gap[level[u]]==0)break;//更新gap数组后若出现断层,直接退出
        level[u]=minlevel+1;//重标号
        gap[level[u]]++;//距离标号为level[u]的点的个数+1
        u=pre[u];//转为当前节点的前驱继续寻找可行弧
    }
    return flow;
}
int main()
{
    int u,v,w,w1,w2;
    while(scanf("%d%d",&m,&n)==2){
        vs=0;
        vt=n-1;
        NV=n;//总点数
        NE=0;
        memset(head,-1,sizeof(head));
        for(int i=1; i<=m; i++) {
            scanf("%d%d%d",&u,&v,&w);
            Insert(--u,--v,w);
        //Insert(v,u,w);
        }
        printf("%d\n",SAP(vs,vt));
    }
    return 0;
}
poj3281Dining 
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 5005;
const int M = 30005;
struct Enode{
    int y, c, next;
} e[M * 2];
//struct Point
//{
//	int son,cur,pre,lim,d;
//} a[maxn];
int n, m, tot, head[N], now[N], h[N], vh[N], augc, found, flow;
//h为距离,vh为距离标号的数目,augc当前边的流量,flow总流量
//int tot, n, m, st, ed, cnt[maxn];
void Addedge(int x, int y, int c){
    e[++tot].y = y; e[tot].c = c; e[tot].next = head[x]; head[x] = tot;
    e[++tot].y = x; e[tot].c = 0; e[tot].next = head[y]; head[y] = tot;
}
int f,d;
void Init(){
    int x, y, c, i,tmp;
    scanf("%d%d%d", &n, &f, &d);
    tot = -1; memset(head, -1, sizeof(head));
    for (i = 1; i <= n; i++){
        scanf("%d%d", &x, &y);
        while(x--){
            scanf("%d",&tmp);
            Addedge(0,tmp,1);
            //Addedge(tmp,f+i,1);在这里可能会加重了
        }
        Addedge(i+f,i+f+n,1);
        while(y--){
            scanf("%d",&tmp);
            Addedge(i+n+f,f+2*n+tmp,1);
            //Addedge(f+2*n+tmp,f+2*n+d+1,1);
        }
    }
    for(int j=1;j<=f;j++) Addedge(j,i+f,1);
    for(int j=1;j<=d;j++) Addedge(j+f+n*2,d+f+n*2+1,1);
    memcpy(now, head, sizeof(head));
}

void Aug(int x, int st, int ed, int n){//x当前点,st源点,ed汇点,n,点最大数目
    int p = now[x], minh = n-1, augco = augc;
    if (x == ed){
        found = 1;
        flow += augc;
        return;
    }
    while (p != -1){//每条边
        if (e[p].c > 0 && h[e[p].y] + 1 == h[x]){//e[p].y下一点,可行弧
            augc = min(augc, e[p].c);
            Aug(e[p].y, st, ed, n);//下一点
            if (h[st] >= n) return;
            if (found) break;
            augc = augco;
        }
        p = e[p].next;//下一边
    }
    if (found){
        e[p].c -= augc;
        e[p ^ 1].c += augc;//边标号从0开始,反向边,p^1
    }else{
        p = head[x];
        while (p != -1){
            if (e[p].c > 0 && h[e[p].y] < minh){
                minh = h[e[p].y];//最小标号
                now[x] = p;
            }
            p = e[p].next;
        }
        vh[h[x]] --;//标号数目-1
        if (!vh[h[x]]) h[st] = n;//无符合标号
        h[x] = minh + 1;//修改当前节点标号
        vh[h[x]] ++;
    }
}
void Maxflow(int st, int ed, int n){
    flow = 0;
    memset(h, 0, sizeof(h));
    memset(vh, 0, sizeof(vh));
    vh[0] = n;//距离标号的数目
    while (h[st] < n){
        found = 0;
        augc = 1 << 30;//假设流量无穷大
        Aug(st, st, ed, n);
    }
}
int main(){
    Init();
    Maxflow(0, 2*n+f+d+1, 2*n+f+d+2);
    printf("%d\n", flow);
}
sgu185 
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define MAXN 444
#define MAXM 160002 //邻接表要开边数的2倍
#define inf 1<<30
#define INF 1<<30
#define maxn 444
bool vis[maxn];
int adj[maxn][maxn],dis[maxn];//pre[]记录前驱
int n, m;
void dijkstra(int v)
{
    int i, j, u , min;
    for(i=0;i<=n;i++)  {
        dis[i]=adj[v][i];
        vis[i]=0;
    }
    vis[v]=1;dis[v]=0;
    for(i=1;i<=n;i++)   {
        min = INF;
        for(j=1;j<=n;j++)
            if(!vis[j]&&min > dis[j]){
                min = dis[j];
                u = j;
            }
        if(min == INF)break;
        vis[u]=1;
        for(j=1;j<=n;j++){
            if(!vis[j]&&adj[u][j]!=INF&&dis[u]+adj[u][j]<dis[j])
                dis[j] = adj[u][j] + dis[u];
        }
    }
}
struct Edge {
    int v,cap,next;
} edge[MAXM];

int vs,vt,NE,NV;
int head[MAXN];

void Insert(int u,int v,int cap)
{
    edge[NE].v=v;
    edge[NE].cap=cap;
    edge[NE].next=head[u];
    head[u]=NE++;

    edge[NE].v=u;//反向边
    edge[NE].cap=0;
    edge[NE].next=head[v];
    head[v]=NE++;
}
void make_graph(){
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++){
        if(i==j||adj[i][j]==INF) continue;
        if(dis[i]+adj[i][j]==dis[j]){
            Insert(i-1,j-1,1);
            //Insert(j-1,i-1,1);
            //printf("edge:%d %d\n",i,j);
        }
    }
}
int level[MAXN];//标记层次（距离标号）
int gap[MAXN];
//间隙优化，定义gap[i]为标号是i的点的个数
//在重标记i时，检查gap[level[i]],若减为0,这算法结束。

void bfs(int vt)
{
    memset(level,-1,sizeof(level));
    memset(gap,0,sizeof(gap));
    level[vt]=0;
    gap[level[vt]]++;
    queue<int>que;
    que.push(vt);
    while(!que.empty()) {
        int u=que.front();
        que.pop();
        for(int i=head[u]; i!=-1; i=edge[i].next) {
            int v=edge[i].v;
            if(level[v]!=-1)continue;
            level[v]=level[u]+1;
            gap[level[v]]++;
            que.push(v);
        }
    }
}
int pre[MAXN];//前驱
int cur[MAXN];

int SAP(int vs,int vt)//源点,汇点
{
    bfs(vt);//bfs优化,
    memset(pre,-1,sizeof(pre));
    memcpy(cur,head,sizeof(head));//cur保存的是当前弧
    int u=pre[vs]=vs,flow=0,aug=inf;//源点前驱还是它本身,aug表示增广路的可改尽量
    gap[0]=NV;
    while(level[vs]<NV) {
        bool flag=false;
        for(int &i=cur[u]; i!=-1; i=edge[i].next) {
            int v=edge[i].v;//v是u的后继
            //寻找可行弧
            if(edge[i].cap&&level[u]==level[v]+1) {
                flag=true;//找到可行弧
                pre[v]=u;//记录前驱
                u=v;
                //  aug=(aug==-1?edge[i].cap:min(aug,edge[i].cap));
                aug=min(aug,edge[i].cap);
                //如果找到 一条增广路
                if(v==vt) {
                    flow+=aug;//更新最大流
                    for(u=pre[v]; v!=vs; v=u,u=pre[u]) {
                        edge[cur[u]].cap-=aug;//前向弧容量减少
                        edge[cur[u]^1].cap+=aug;//后向弧容量增加
                    }
                    //     aug=-1;
                    aug=inf;
                }
                break;
            }
        }
        if(flag)continue;
        int minlevel=NV;
        //寻找与当前点相连接的点中最小的距离标号(重标号)
        for(int i=head[u]; i!=-1; i=edge[i].next) {
            int v=edge[i].v;
            if(edge[i].cap&&level[v]<minlevel) {
                minlevel=level[v];
                cur[u]=i;
            }
        }
        if(--gap[level[u]]==0)break;//更新gap数组后若出现断层,直接退出
        level[u]=minlevel+1;//重标号
        gap[level[u]]++;//距离标号为level[u]的点的个数+1
        u=pre[u];//转为当前节点的前驱继续寻找可行弧
    }
    //cout<<'f'<<flow<<endl;
    return flow;
}
bool flag;//flag的作用
void dfs(int u){///输出路径
    if(flag) return ;
    if(u==n-1){
        flag=true;
        cout<<endl;
        return ;
    }
    for(int i=head[u];i!=-1&&!flag;i=edge[i].next){
        int v=edge[i].v;
        if(edge[i].cap==0 && !(i&1)){///
            edge[i].cap=-1;
            cout<<' '<<v+1;
            dfs(v);
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=0;i<=n;i++)for(int j=0;j<=n;j++){
        if(i==j)adj[i][j]=0;
        else adj[i][j]=INF;
    }
    for(int i=1;i<=m;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        adj[x][y]=adj[y][x]=z;
    }
    dijkstra(1);
    //cout<<' '<<dis[n]<<endl;
    vs=0;
    vt=n-1;
    NV=n;//总点数
    NE=0;
    memset(head,-1,sizeof(head));
    make_graph();//位置= =
    if(SAP(vs,vt)<2)printf("No solution\n");
    else{
        cout<<1;
        flag=false;
        dfs(0);
        //cout<<'a'<<endl;
        cout<<1;
        flag=false;
        dfs(0);
    }
    return 0;
}


coder_coder_coder
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
leetcode下载-algorithm-1:力扣、HDU、ZOJ、POJ
06-29
Interview，ZOJ，POJ 等平台。 欢迎Coders对代码加以指正和提议！ 常见问题总结 两整数求平均值 average = min + (max - min) / 2 防止两整数的和越界 整数乘积对比 1.0 * m * m == num 类似乘积对比， 需转为double...
hdu oj poj 题目代码与详解
09-08
【题目代码与详解】是针对在线编程竞赛平台HDU（杭州电子科技大学在线评测系统）和POJ（普林斯顿大学在线判题系统）中的经典题目所编写的解决方案和解析的集合。这些平台是程序员和计算机科学学生提升算法技能、熟悉...
【转载】图论 500题——主要为hdu/poj/zoj
weixin_30448603的博客
06-25 177
转自——http://blog.csdn.net/qwe20060514/article/details/8112550 =============================以下是最小生成树+并查集======================================【HDU】1213 How Many Tables基础并查集★1272 小希的迷宫基础并查集...
POJ 图论【转自夏天的风blog】
Ⅶ的专栏
07-24 1731
POJ 图论【转自夏天的风blog】 目录(?)[+] 杂题 POJ 2449 Remmarguts' Date(中等)  AC http://acm.pku.edu.cn/JudgeOnline/problem?id=2449 题意：经典问题：K短路 解法：dijkstra+A*(rec)，方法很多 相关：http://acm
POJ 图论
07-16 4993
POJ 2449 Remmarguts' Date(中等)  AC http://acm.pku.edu.cn/JudgeOnline/problem?id=2449 题意：经典问题：K短路 解法：dijkstra+A*(rec)，方法很多 相关：http://acm.pku.edu.cn/JudgeOnline/showcontest?contest_id=1144
POJ图论列表
新博客地址:sensirly.github.io
04-09 1326
做题记录：     AC           典型    好题或较难         待做            WA        多次wa后ac 1062* 昂贵的聘礼 枚举+dijkstra加点构图 1087A Plug for UNIX 2分匹配+Floyd  1094 Sorting It All Out floyd +拓扑 1112* Team Them Up! 2分图染
POJ 图论 基本概括
muxi@Achilles的专栏
08-27 5130
转载出处：http://blog.csdn.net/shahdza/article/details/6632538   POJ 2449 Remmarguts' Date(中等)  AC http://acm.pku.edu.cn/JudgeOnline/problem?id=2449 题意：经典问题：K短路 解法：dijkstra+A*(rec)，方法很多 相关
图论500题
Gzbgreat的博客
07-14 407
膜...
POJ图论
梦想起航的地方
01-16 2055
POJ 2449 Remmarguts' Date(中等)  AC http://acm.pku.edu.cn/JudgeOnline/problem?id=2449 题意：经典问题：K短路 解法：dijkstra+A*(rec)，方法很多 相关：http://acm.pku.edu.cn/JudgeOnline/showcontest?contest_id=1144
poj 图论！！！
代码随想录
12-08 2162
POJ 2449 Remmarguts' Date(中等)  AC http://acm.pku.edu.cn/JudgeOnline/problem?id=2449 题意：经典问题：K短路 解法：dijkstra+A*(rec)，方法很多 相关：http://acm.pku.edu.cn/JudgeOnline/showcontest?contest_id=1144
各大Oj经典图论500题
多反思，多回顾，要坚持。
03-10 4750
=============================以下是最小生成树+并查集====================================== 【HDU】 *1213         How Many Tables            基础并查集★ *1272         小希的迷宫            基础并查集★ *1325&&poj1308     Is It
HDU&&POJ图论题集
04-24
图论〔Graph Theory〕是数学的一个分支。它以图为研究对象。图论中的图是由若干给定的点及连接两点的线所构成的图形，这种图形通常用来描述某些事物之间的某种特定关系，用点代表事物，用连接两点的线表示相应两个...
jihe.rar_2289_POJ 3714_poj3714_poj3714 Ra_visual c
09-22
标题中的"jihe.rar_2289_POJ 3714_poj3714_poj3714 Ra_visual c" 提到了一个压缩文件，可能包含有关编程竞赛或算法解决的资源，特别是与POJ（Problem On Judge）平台上的问题3714相关的。"Ra_visual c"可能指的是...
nginx资源useless，don't download
07-09
nginx资源useless，don't download
util-linux-2.37.2.tar.xz
07-09
util-linux-2.37.2.tar.xz
“人力资源+大数据+薪酬报告+涨薪调薪”
最新发布
07-09
人力资源+大数据+薪酬报告+涨薪调薪，在学习、工作生活中，越来越多的事务都会使用到报告，通常情况下，报告的内容含量大、篇幅较长。那么什么样的薪酬报告才是有效的呢？以下是小编精心整理的调薪申请报告，欢迎大家分享。相信老板看到这样的报告，一定会考虑涨薪的哦。
API-pour-les-nuls-WSM14025FRFR-3-of-5.pdf
07-09
api
acm部分模板
01-29
图论 dp 部分模板 内涵几道poj以及hdu上的例题..............................................................

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值