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

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;
}