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的街道（单向）。
N ≤ 200，M ≤ 20000。
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

Source

Day1

几乎就是一个裸的费用流，将路程设为费用，最大周期就是求最大流。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
const int N=205;
const int M=20005;
const int inf=1e9+7;
int n,m,s,t,cnt=1,ans,sum,hd[2*N],pre[2*M+2*N],dis[2*N];
bool inq[2*N];
queue<int>q;
struct edge
{
    int to,nxt,f,w;
}v[2*M+2*N];
void addedge(int x,int y,int z,int w)
{
    v[++cnt].to=y,v[cnt].f=z,v[cnt].w=w;
    v[cnt].nxt=hd[x],hd[x]=cnt;
}
void addedges(int x,int y,int z,int w)
{
    addedge(x,y,z,w),addedge(y,x,0,-w);
}
bool spfa()
{
    memset(pre,0,sizeof(pre));
    memset(dis,0x3f,sizeof(dis));
    dis[s]=0;
    q.push(s);
    inq[s]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        inq[u]=0;
        for(int i=hd[u];i;i=v[i].nxt)
            if(v[i].f&&dis[v[i].to]>dis[u]+v[i].w)
            {
                dis[v[i].to]=dis[u]+v[i].w;
                pre[v[i].to]=i;
                if(!inq[v[i].to])
                {
                    inq[v[i].to]=1;
                    q.push(v[i].to);
                }
            }
    }
    return pre[t];
}
int main()
{
    scanf("%d%d",&n,&m);
    s=1,t=2*n;
    for(int i=1;i<=n;i++)
        if(i==1||i==n)
            addedges(i,i+n,inf,0);
        else
            addedges(i,i+n,1,0);
    for(int i=1;i<=m;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        addedges(x+n,y,1,z);
    }
    while(spfa())
    {
        int flw=inf;
        for(int i=pre[t];i;i=pre[v[i^1].to])
            flw=min(flw,v[i].f);
        sum+=flw;
        ans+=dis[t]*flw;
        for(int i=pre[t];i;i=pre[v[i^1].to])
            v[i].f-=flw,v[i^1].f+=flw;
    }
    printf("%d %d\n",sum,ans);
    return 0;
}