bzoj 1877: [SDOI2009]晨跑

1877: [SDOI2009]晨跑

Time Limit: 4 Sec   Memory Limit: 64 MB
Submit: 1747   Solved: 924
[ Submit][ Status][ Discuss]

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。

Source

Day1

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define N 20003
using namespace std;
int n,m,point[N*2],next[N*4],v[N*4],remain[N*4],c[N*4],tot=-1;
int minn,flow,dis[N],can[N],laste[N];
const int inf=1e9;
void add(int x,int y,int z,int k)
{
  tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; remain[tot]=z; c[tot]=k;
  tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; remain[tot]=0; c[tot]=-k;
}
int addflow(int s,int t)
{
  int now=t; int ans=inf;
  while(now!=s)
  {
  	ans=min(ans,remain[laste[now]]);
  	now=v[laste[now]^1];
  }
  now=t;
  while (now!=s)
   {
   	 remain[laste[now]]-=ans;
   	 remain[laste[now]^1]+=ans;
   	 now=v[laste[now]^1];
   }
  return ans;
}
bool spfa(int s,int t)
{
  memset(dis,0x7f,sizeof(dis));
  dis[s]=0; can[s]=1;
  queue<int> p;
  p.push(s);
  while (!p.empty())
  {
   int now=p.front(); p.pop();
   for (int i=point[now];i!=-1;i=next[i])
    if (dis[v[i]]>dis[now]+c[i]&&remain[i])
     {
       dis[v[i]]=dis[now]+c[i];
       laste[v[i]]=i;
       if (!can[v[i]])
       {
       	 can[v[i]]=1;
       	 p.push(v[i]);
       }
     }
    can[now]=0;
  }
  if (dis[t]>inf) return false;
  int tl=addflow(s,t);
  flow+=tl;
  minn+=tl*dis[t];
  return true;
}
void maxflow(int s,int t)
{
  flow=0; minn=0;
  while (spfa(s,t));
}
int main()
{
 memset(next,-1,sizeof(next));
 memset(point,-1,sizeof(point));
 scanf("%d%d",&n,&m);
 for (int i=2;i<=n-1;i++)
  add(i,i+n,1,0);
 for (int i=1;i<=m;i++)
  {
    int x,y,z; scanf("%d%d%d",&x,&y,&z);
    add(x+n,y,1,z);
  }
 add(1,n+1,inf,0); add(n,n+n,inf,0);
 maxflow(1,n+n);
 printf("%d %d\n",flow,minn);
}