2521: [Shoi2010]最小生成树

2521: [Shoi2010]最小生成树

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 463   Solved: 266
[ Submit][ Status][ Discuss]

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

Sample Output

1

HINT

第1个样例就是问题描述中的例子。

1<=n<=500,1<=M<=800,1<=D<10^6

Source

day2

[ Submit][ Status][ Discuss]



假设这条边是(x,y)，那这条边一定在最小生成树的充要条件为，

所有权值小于等于它的边无法使(x,y)连通

题目中的操作等价于选定一条边使它的权值+1

那么建一张图，只取原图中权值不超过关建边的边，容量为差值+1，代表删除它的代价

跑一遍最小割就行了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
 
const int maxn = 505;
const int maxm = 4040;
const int INF = 1E9 + 233;
 
struct E{
    int to,cap,flow; E(){}
    E(int to,int cap,int flow): to(to),cap(cap),flow(flow){}
}edgs[maxm];
 
int n,m,s,t,Lab,cnt,cur[maxn],L[maxn],A[maxm],B[maxm],w[maxm];
 
vector <int> v[maxn];
queue <int> Q;
 
void Add(int x,int y,int cap)
{
    v[x].push_back(cnt); edgs[cnt++] = E(y,cap,0);
    v[y].push_back(cnt); edgs[cnt++] = E(x,0,0);
}
 
bool BFS()
{
    for (int i = 1; i <= n; i++) L[i] = 0;
    Q.push(s); L[s] = 1;
    while (!Q.empty())
    {
        int k = Q.front(); Q.pop();
        for (int i = 0; i < v[k].size(); i++)
        {
            E e = edgs[v[k][i]];
            if (e.cap == e.flow || L[e.to]) continue;
            L[e.to] = L[k] + 1; Q.push(e.to);
        }
    }
    return L[t];
}
 
int Dinic(int x,int a)
{
    if (x == t) return a; int flow = 0;
    for (int &i = cur[x]; i < v[x].size(); i++)
    {
        E &e = edgs[v[x][i]];
        if (e.cap == e.flow || L[e.to] != L[x] + 1) continue;
        int f = Dinic(e.to,min(a,e.cap - e.flow));
        if (!f) continue; e.flow += f; flow += f;
        edgs[v[x][i]^1].flow -= f; a -= f;
        if (!a) return flow;
    }
    if (!flow) L[x] = -1; return flow;
}
 
int main()
{
    #ifdef DMC
        freopen("DMC.txt","r",stdin);
    #endif
     
    cin >> n >> m >> Lab;
    for (int i = 1; i <= m; i++)
        scanf("%d%d%d",&A[i],&B[i],&w[i]);
    s = A[Lab]; t = B[Lab];
    for (int i = 1; i <= m; i++)
        if (i != Lab && w[i] <= w[Lab])
        {
            Add(A[i],B[i],w[Lab] - w[i] + 1);
            Add(B[i],A[i],w[Lab] - w[i] + 1);
        }
    int MaxFlow = 0;
    while (BFS())
    {
        for (int i = 1; i <= n; i++) cur[i] = 0;
        MaxFlow += Dinic(s,INF);
    }
    cout << MaxFlow << endl;
    return 0;
}