其实之前写过一个板子,但是一点印象都没有,所以今天重写了一下,顺便把这个题当成板子就行了。
其实费用流就是把bfs换成spfa,但是中间有一个原则,就是费用优先,在费用(就是c)上跑spfa,顺便求出流量。
其实理解起来还算简单,就是先spfa找最小费用路径,然后在路径上找能跑的最大流。
题干:
Description
Elaxia最近迷恋上了空手道,他为自己设定了一套健身计划,比如俯卧撑、仰卧起坐等 等,不过到目前为止,他
坚持下来的只有晨跑。 现在给出一张学校附近的地图,这张地图中包含N个十字路口和M条街道,Elaxia只能从 一
个十字路口跑向另外一个十字路口,街道之间只在十字路口处相交。Elaxia每天从寝室出发 跑到学校,保证寝室
编号为1,学校编号为N。 Elaxia的晨跑计划是按周期(包含若干天)进行的,由于他不喜欢走重复的路线,所以
在一个周期内,每天的晨跑路线都不会相交(在十字路口处),寝室和学校不算十字路 口。Elaxia耐力不太好,
他希望在一个周期内跑的路程尽量短,但是又希望训练周期包含的天 数尽量长。 除了练空手道,Elaxia其他时间
都花在了学习和找MM上面,所有他想请你帮忙为他设计 一套满足他要求的晨跑计划。
Input
第一行:两个数N,M。表示十字路口数和街道数。
接下来M行,每行3个数a,b,c,表示路口a和路口b之间有条长度为c的街道(单向)。
N ≤ ,M ≤ 。
Output
两个数,第一个数为最长周期的天数,第二个数为满足最长天数的条件下最短的路程长 度。
Sample Input Sample Output HINT
Source
代码:
#include<iostream>
#include<cstdio>
#include<cmath>
#include<ctime>
#include<queue>
#include<algorithm>
#include<cstring>
using namespace std;
#define duke(i,a,n) for(int i = a;i <= n;i++)
#define lv(i,a,n) for(int i = a;i >= n;i--)
#define clean(a) memset(a,0,sizeof(a))
const int INF = << ;
typedef long long ll;
typedef double db;
template <class T>
void read(T &x)
{
char c;
bool op = ;
while(c = getchar(), c < '' || c > '')
if(c == '-') op = ;
x = c - '';
while(c = getchar(), c >= '' && c <= '')
x = x * + c - '';
if(op) x = -x;
}
template <class T>
void write(T x)
{
if(x < ) putchar('-'), x = -x;
if(x >= ) write(x / );
putchar('' + x % );
}
const int N = ;
struct node
{
int l,r,c,f,nxt,other;
}a[N * ];
int len = ,lst[N],st,ed;
void add(int x,int y,int f,int c)
{
int k1,k2 = ;
a[++len].l = x;
a[len].r = y;
a[len].c = c;
a[len].f = f;
a[len].nxt = lst[x];
lst[x] = len;
k1 = len;
a[++len].l = y;
a[len].r = x;
a[len].c = -c;
a[len].f = ;
a[len].nxt = lst[y];
lst[y] = len;
k2 = len;
a[k1].other = k2;
a[k2].other = k1;
}
int n,m,dis[N],flow[N],vis[N];
int pre[N],head[N],ans1 = ,ans2 = ;
queue <int> qu;
bool spfa()
{
clean(vis);
memset(dis,,sizeof(dis));
memset(flow,,sizeof(flow));
memset(pre,,sizeof(pre));
qu.push(st);
vis[st] = ;
dis[st] = ;
pre[st] = ;
while(!qu.empty())
{
int x = qu.front();
qu.pop();
for(int k = lst[x];k;k = a[k].nxt)
{
if(a[k].f == )continue;
int y = a[k].r;
if(dis[y] > dis[x] + a[k].c)
{
dis[y] = dis[x] + a[k].c;
pre[y] = x;head[y] = k;
flow[y] = min(flow[x],a[k].f);
if(!vis[y])
vis[y] = ,qu.push(y);
}
}
vis[x] = ;
}
return pre[ed];
}
void mvp()
{
// cout<<"QAQ"<<endl;
int x = flow[ed];
int i = head[ed];
ans1 ++;
/*duke(j,1,2 * n)
printf("%d ",head[j]);*/
i = head[ed];
while(i)
{
// cout<<i<<endl;
ans2 += x * a[i].c;
a[i].f -= x;
a[a[i].other].f += x;
i = head[a[i].l];
// cout<<a[i].l<<endl;
}
}
int main()
{
read(n);read(m);
st = ;ed = n * ;
duke(i,,m)
{
int u,v,w;
read(u);read(v);read(w);
add(u + n,v,,w);
}
head[a[].l] = ;
duke(i,,n - )
add(i,i + n,,);
add(,st + n,INF,);
add(n,ed,INF,);
while(spfa()) mvp();
printf("%d %d\n",ans1,ans2);
return ;
}
/*
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
*/