Description
Secsa最近对最小生成树问题特别感兴趣。他已经知道如果要去求出一个n个点、m条边的无向图的最小生成树有一个Krustal算法和另一个Prim的算法。另外,他还知道,某一个图可能有多种不同的最小生成树。例如,下面图 3中所示的都是图 2中的无向图的最小生成树:
当然啦,这些都不是今天需要你解决的问题。Secsa想知道对于某一条无向图中的边AB,至少需要多少代价可以保证AB边在这个无向图的最小生成树中。为了使得AB边一定在最小生成树中,你可以对这个无向图进行操作,一次单独的操作是指:先选择一条图中的边 P1P2,再把图中除了这条边以外的边,每一条的权值都减少1。如图 4所示就是一次这样的操作:
Input
输入文件的第一行有3个正整数n、m、Lab分别表示无向图中的点数、边数、必须要在最小生成树中出现的AB边的标号。
接下来m行依次描述标号为1,2,3…m的无向边,每行描述一条边。每个描述包含3个整数x、y、d,表示这条边连接着标号为x、y的点,且这条边的权值为d。
输入文件保证1<=x,y<=N,x不等于y,且输入数据保证这个无向图一定是一个连通图。
Output
输出文件只有一行,这行只有一个整数,即,使得标号为Lab边一定出现最小生成树中的最少操作次数。
Sample Input
4 6 1
1 2 2
1 3 2
1 4 3
2 3 2
2 4 4
3 4 5
1 2 2
1 3 2
1 4 3
2 3 2
2 4 4
3 4 5
Sample Output
1
HINT
第1个样例就是问题描述中的例子。
1<=n<=500,1<=M<=800,1<=D<10^6
题解:
考虑一条边一定在最小生成树中意味着去掉这条边之后,权值小于等于它的边不能使这条边的两端点联通
固定一条边不变把剩下的边都减一可以看成把这条边加一.
那割掉一条边的代价就是v[x]-v[i]+1;
然后直接最小割即可.
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 2010
#define M 100010
#define inf 210000000
using namespace std;
int point[N],next[M<<1],x,y,v,n,m,S,T,V,cnt=1;
int cur[N],dis[N],pre[N],gap[N],pos;
struct use{
int st,en,v,id;
}e[M<<1],b[M];
void add(int x,int y,int v){
//cout<<x<<' '<<y<<' '<<v<<endl;
next[++cnt]=point[x];point[x]=cnt;
e[cnt].st=x;e[cnt].en=y;e[cnt].v=v;
next[++cnt]=point[y];point[y]=cnt;
e[cnt].st=y;e[cnt].en=x;e[cnt].v=v;
}
int read(){
int x(0);char ch=getchar();
while (ch<'0'||ch>'9') ch=getchar();
while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x;
}
bool cmp(use a,use b){
return a.v<b.v;
}
int isap(int ss,int tt){
int i,ans(0),mn,u=ss;
gap[0]=n;
for (i=1;i<=n;i++) cur[i]=point[i];
while (dis[ss]<n){
bool f=false;
for (i=cur[u];i;i=next[i])
if (dis[e[i].en]+1==dis[u]&&e[i].v){f=true;cur[u]=i;break;}
if (f){
pre[u=e[i].en]=i;
if (u==tt){
mn=inf;
for (i=tt;i!=ss;i=e[pre[i]].st) mn=min(mn,e[pre[i]].v);
ans+=mn;
for (i=tt;i!=ss;i=e[pre[i]].st) e[pre[i]].v-=mn,e[pre[i]^1].v+=mn;
u=ss;
}
}
else{
--gap[dis[u]];if (!gap[dis[u]]) return ans;
for (mn=n,i=point[u];i;i=next[i]) if (e[i].v) mn=min(mn,dis[e[i].en]);
gap[dis[u]=mn+1]++;cur[u]=point[u];if (u!=ss) u=e[pre[u]].st;
}
}
return ans;
}
int main(){
//freopen("a.in","r",stdin);
n=read();m=read();pos=read();
for (int i=1;i<=m;i++){
x=read();y=read();v=read();
b[i].st=x;b[i].en=y;b[i].v=v;b[i].id=i;
if (i==pos) S=x,T=y,V=v;
}
sort(b+1,b+m+1,cmp);
for (int i=1;i<=m;i++){
if (b[i].v>V) break;
if (b[i].id==pos) continue;
add(b[i].st,b[i].en,V-b[i].v+1);
}
cout<<isap(S,T);
}