奥妙重重的网络流

原创 2016年08月30日 15:26:58

这两个代码都是BZOJ1532的代码,仅仅是建图的差距就导致了4倍时间差。至今不懂为什么,先留下代码。

#include <bits/stdc++.h>
#define v edge[i].to
#define RepG(i,x) for(int i = head[x];~ i;i = edge[i].next)
#define fl edge[i].f
#define vfl edge[i ^ 1].f
#define Rep(i,n) for(int i = 1;i <= n;++ i)
using namespace std;
const int inf = 1 << 29;
const int N = 20005;
const int M = 400005;
int S = 0,T = N - 1,head[N],cnt,dep[N],q[N],ans,n;
struct Edge{int next,to,f;}edge[M];
void save(int a,int b,int c){edge[cnt]=(Edge){head[a],b,c},head[a] = cnt ++;edge[cnt] = (Edge){head[b],a,0},head[b] = cnt ++;}
bool bfs()
{
    int h,t,x;
    memset(dep,-1,sizeof(dep));
    dep[q[h = t = 0] = S] = 0;
    while(h <= t)
    {
        x = q[h ++];
        RepG(i,x)if(dep[v] == -1 && fl)dep[q[++ t] = v] = dep[x] + 1;
    }
    return dep[T] > 0;
}
int dfs(int x,int f)
{
    if(x == T)return f;
    int s = 0,ff;
    RepG(i,x)if(dep[v] == dep[x] + 1 && f && fl && (ff = dfs(v,min(f,fl))))
        s += ff,fl -= ff,vfl += ff,f -=ff;
    if(!s)dep[x] = -1;
    return s;
}
int mxf(){int s = 0;for(;bfs();s += dfs(S,inf));return s;}
int ax[N],ay[N],m;
void ReBuild(int ans)
{
    memset(head,-1,sizeof(head));cnt = 0;
    Rep(i,n)save(i,T,ans);
    Rep(i,m)
    {
        save(S,n + i, 1);
        save(n + i,ax[i], 1);
        save(n + i,ay[i], 1);
    }
} 
bool check(){return mxf() < m;}
int main()
{
    scanf("%d%d",&n,&m);
    Rep(i,m)scanf("%d%d",&ax[i],&ay[i]);
    int l = 0,r = 10000,mid;
    while(l < r)
    {
        mid = l + r >> 1;
        ReBuild(mid);
        if(check())l = mid + 1;
        else r = mid;
    }
    printf("%d\n",l);
    return 0;
}
//slow
#include <bits/stdc++.h>
#define v edge[i].to
#define RepG(i,x) for(int i = head[x];~ i;i = edge[i].next)
#define fl edge[i].f
#define vfl edge[i ^ 1].f
#define Rep(i,n) for(int i = 1;i <= n;++ i)
using namespace std;
const int inf = 1 << 29;
const int N = 20005;
const int M = 400005;
int S = 0,T = N - 1,head[N],cnt,dep[N],q[N],ans,n;
struct Edge{int next,to,f;}edge[M];
void save(int a,int b,int c){edge[cnt]=(Edge){head[a],b,c},head[a] = cnt ++;edge[cnt] = (Edge){head[b],a,0},head[b] = cnt ++;}
bool bfs()
{
    int h,t,x;
    memset(dep,-1,sizeof(dep));
    dep[q[h = t = 0] = S] = 0;
    while(h <= t)
    {
        x = q[h ++];
        RepG(i,x)if(dep[v] == -1 && fl)dep[q[++ t] = v] = dep[x] + 1;
    }
    return dep[T] > 0;
}
int dfs(int x,int f)
{
    if(x == T)return f;
    int s = 0,ff;
    RepG(i,x)if(dep[v] == dep[x] + 1 && f && fl && (ff = dfs(v,min(f,fl))))
        s += ff,fl -= ff,vfl += ff,f -=ff;
    if(!s)dep[x] = -1;
    return s;
}
int mxf(){int s = 0;for(;bfs();s += dfs(S,inf));return s;}
int ax[N],ay[N],m;
void ReBuild(int ans)
{
    memset(head,-1,sizeof(head));cnt = 0;
    Rep(i,n)save(S, i,ans);
    Rep(i,m)
    {
        save(n + i, T, 1);
        save(ax[i], n + i, 1);
        save(ay[i], n + i, 1);
    }
} 
bool check(){return mxf() < m;}//答案会比它大。 
int main()
{
    scanf("%d%d",&n,&m);
    Rep(i,m)scanf("%d%d",&ax[i],&ay[i]);
    int l = 0,r = 10000,mid;
    while(l < r)
    {
        mid = l + r >> 1;
        ReBuild(mid);
        if(check())l = mid + 1;
        else r = mid;
    }
    printf("%d\n",l);
    return 0;
}
//fast

很妙啊。
在都是重边的时候,第一份代码明显慢于第二份(试试就知道了)。

版权声明:呐,转载请交稿费QAQ(*随意转载啦)

日常训练 腐女的生日 奥妙重重的最短路

题意简述:平面上给定nn个矩形建筑物,问从(0,0)(0,0)走到(x0,y0)(x0,y0)的最短距离,只能平行与坐标轴行走且不能穿过建筑。(n≤105 ,y (n\le 10^5\ ,y\ 坐标∈...

《Linux那些事儿之我是USB》我是U盘(34)迷雾重重的批量传输(三)

在usb_stor_Bulk_transport()中,这个函数中调用的第一个最重要的函数,那就是usb_stor_bulk_transfer_buf()。仍然是来自drivers/usb/stroa...

《Linux那些事儿之我是USB》我是U盘(37)迷雾重重的批量传输(六)

usb_stor_bulk_transfer_sglist()函数有一定的“蛊惑性”,我们前面说过,之所以采用sglist,就是为了提高传输效率。我们更知道,sg的目的就是让一堆不连续的buffers...

《Linux那些事儿之我是USB》我是U盘(35)迷雾重重的批量传输(四)

有时候我也被这个问题所困扰,我不知道是我不明白,还是这世界变化太快。连Linux中都引入了过期这么一个概念。设置一个时间,如果时间到了该做的事情还没有做完,那么某些事情就会发生。 比如需要烤蛋糕,现...

问题重重的WINDOWS SERVER 2008+IIS7.0 站点发布

今天在服务器上发布一个WEB站点,折腾了一个上午。我把过程中碰到的问题记录下来,以备以后参考。 这次发布的服务器环境是WINDOWS SERVER 2008(x64)+SQL SERVER 200...

线性规划与网络流题解

  • 2012年09月26日 10:31
  • 1.92MB
  • 下载

线性规划.网络流.二分图匹配

  • 2012年06月09日 11:50
  • 401KB
  • 下载

网络流算法Dinic的Python实现

在上一篇我们提到了网络流算法Push-relabel,那是90年代提出的算法,算是比较新的,而现在要说的Dinic算法则是由以色列人Dinitz在冷战时期,即60-70年代提出的算法变种而来的,其算法...

网络流 资料整理(精华)

  • 2013年02月18日 13:35
  • 3.35MB
  • 下载

基于道路的网络流模型

  • 2011年09月05日 12:23
  • 571KB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:奥妙重重的网络流
举报原因:
原因补充:

(最多只允许输入30个字)