[BZOJ 1877][SDOI2009]晨跑（费用流）

最新推荐文章于 2018-09-01 23:21:37 发布

Zars19

最新推荐文章于 2018-09-01 23:21:37 发布

阅读量266

点赞数

分类专栏：网络流

本文链接：https://blog.csdn.net/zars19/article/details/60315014

版权

网络流专栏收录该内容

5 篇文章 0 订阅

订阅专栏

Description

Elaxia最近迷恋上了空手道，他为自己设定了一套健身计划，比如俯卧撑、仰卧起坐等等，不过到目前为止，他坚持下来的只有晨跑。现在给出一张学校附近的地图，这张地图中包含N个十字路口和M条街道，Elaxia只能从一个十字路口跑向另外一个十字路口，街道之间只在十字路口处相交。Elaxia每天从寝室出发跑到学校，保证寝室编号为1，学校编号为N。 Elaxia的晨跑计划是按周期（包含若干天）进行的，由于他不喜欢走重复的路线，所以在一个周期内，每天的晨跑路线都不会相交（在十字路口处），寝室和学校不算十字路口。Elaxia耐力不太好，他希望在一个周期内跑的路程尽量短，但是又希望训练周期包含的天数尽量长。除了练空手道，Elaxia其他时间都花在了学习和找MM上面，所有他想请你帮忙为他设计一套满足他要求的晨跑计划。

Input

第一行：两个数N,M。表示十字路口数和街道数。接下来M行，每行3个数a,b,c，表示路口a和路口b之间有条长度为c的街道（单向）。

Output

两个数，第一个数为最长周期的天数，第二个数为满足最长天数的条件下最短的路程长度。

Sample Input

7 10
1 2 1
1 3 1
2 4 1
3 4 1
4 5 1
4 6 1
2 5 5
3 6 6
5 7 1
6 7 1

Sample Output

2 11
HINT

对于30%的数据，N ≤ 20，M ≤ 120。
对于100%的数据，N ≤ 200，M ≤ 20000。

Solution

裸的费用流啦

1A率跌破谷底…交了一遍超时？
于是自己测了一组极限数据，真的超时，于是发现SPFA的!inq[v]判断加错地方了= =
于是又交，还是超时？
我找了找感觉数组开小了（可这和超时有什么关系啊）
然后又交，又超时了？
最后发现freopen忘删了

拆点建流量为1的边，以保证每个路口只能经过一次
具体实现还是看代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
#include<queue>
#define Max(a,b) (a>b?a:b)
#define Min(a,b) (a<b?a:b)
#define INF 0x3f3f3f3f
using namespace std;
int n,m;
int head[500],cnt=0,Maxflow,dis[500],inq[500],a[500],Mindis,pre[500];
struct Node{
    int from,next,to,cap,cost;
}Edges[100005];
int read()
{
    int f=1,x=0;
    char c=getchar();
    while(c>'9'||c<'0'){
        if(c=='-')f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        x=x*10+c-'0';
        c=getchar();
    }
    return x*f;
}
void addedge(int u,int v,int c,int w)
{
    Edges[cnt].from=u;
    Edges[cnt].to=v;
    Edges[cnt].next=head[u];
    head[u]=cnt;
    Edges[cnt].cap=c;
    Edges[cnt++].cost=w;
}
void EK(int s,int t)
{
    Maxflow=0;
    Mindis=0;
    for(;;)
    {
        memset(a,0,sizeof(a));
        memset(dis,0x3f,sizeof(dis));
        queue<int>q;
        q.push(s);
        dis[s]=0;
        a[s]=INF;
        inq[s]=1;
        while(!q.empty())
        {
            int u=q.front();
            for(int i=head[u];~i;i=Edges[i].next)
            {
                int v=Edges[i].to;
                if(dis[v]>dis[u]+Edges[i].cost&&Edges[i].cap>0)
                {
                    pre[v]=i;
                    a[v]=Min(a[u],Edges[i].cap);
                    dis[v]=dis[u]+Edges[i].cost;
                    if(!inq[v])
                    {
                        inq[v]=1;
                        q.push(v);
                    }
                }
            }
            q.pop();
            inq[u]=0;
        }
        if(!a[t]||dis[t]==INF)return;  
        else 
        {
            int i=t;
            while(i!=s)
            {
                Edges[pre[i]].cap-=a[t];
                Edges[pre[i]^1].cap+=a[t];
                i=Edges[pre[i]].from;
            }
            Mindis+=dis[t];
            Maxflow+=a[t];
        }

    }
}
int main()
{
    memset(head,-1,sizeof(head));
    n=read();m=read();
    for(int i=1;i<=m;i++)
    {
        int a,b,c;
        a=read();
        b=read();
        c=read();
        addedge(a+n,b,1,c);
        addedge(b,a+n,0,-c);
    }
    for(int i=2;i<n;i++)
    {
        addedge(i,i+n,1,0);
        addedge(i+n,i,0,0);
    }
    addedge(1,1+n,INF,0);
    addedge(1+n,1,0,0);
    addedge(n,n+n,INF,0);
    addedge(n+n,n,0,0);
    EK(1,2*n);
    printf("%d %d\n",Maxflow,Mindis);
    return 0;
}