b2324/luogu4542 [ZJOI2011]营救皮卡丘

题意：

       中文题

思路：

首先，如果要探寻新的点，假设i->j，一定是走的路径上比他们小的点过去的，并且在这个条件下，一定是走的最短路。这个过程可以使用floyd完成。对于i<j，假如能够到达，那么就建立一条长度为跑完floyd的dist[i][j]。
那么现在问题就变成了最多使用k条路径覆盖整个dag，并且要求花费最小
很容易联想到路径覆盖问题，在那个问题中，我们将一个点拆成两个点，一个表示入度，一个表示出度，答案就是n-maxflow，因为匹配数每增大1,就意味着有两个路径合二为一（假设开始每个点是单个点的路径）。
在现在这个问题中，出现了一些变化，首先我们要求每个路径从0出发，同时要求结果最小。我们使用费用流，用和最小路径覆盖一样的建图，即一个点拆成两个，一个表示入度，一个表示出度，S与左端点连线，流量为1,花费为0,右端的点用类似的方法与T连边。对于i,j(i<j)，假设i,j存在一条路径，那么就让左端的i与右端的j连一条流量为1，花费为dist[i][j]。这样跑出的结果就是答案。

为什么呢？首先我们从是否是全部覆盖与是否合法上开想，首先我们能够保证答案在全覆盖的情况下是合法的，因为我们跑floyd的时候，保证只能从小的点转移。那么怎么保证全覆盖?题目保证了一定有解，而mcmf保证了最大流，故流量一定是n，即每个点一定在某个路径中，所以一定是全覆盖。

错误及反思：

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<iostream>
#include<stack>
#include<set>
#include<string>
#include<map>
#include<cmath>
#include<bitset>
using namespace std;
#define pb push_back
#define pf push_front
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define mst(x,a) memset(x,a,sizeof(x))
#define all(x) (x).begin(),(x).end()
#define CLOSE ios::sync_with_stdio(false)
typedef pair<int,int> pii;
typedef long long ll;
typedef unsigned long long ull;
typedef vector<int> vi;
#define fi first
#define se second
#define sz(x) ((int)x.size())
#define cl(x) x.clear()
const int N = 100100;
const int M = 200010;
const int inf=0x3f3f3f3f;
const int mod=1000000007;
ll qpow(ll a,ll b){ll ans=1;while(b){if(b&1)ans=ans*a%mod;a=a*a%mod;b/=2;}return ans;}
ll qpow(ll a,ll b,ll c){ll ans=1;while(b){if(b&1)ans=ans*a%c;a=a*a%c;b/=2;}return ans;}

int dist[155][155];

int head[N];
int d[N], incf[N], pre[N], v[N];
int n, m, tot, s, t, maxflow, ans, k;

struct EDGE{
    int to,nex,c,v;
}e[M*2];


void add(int x, int y, int z, int c) {
    e[tot].to=y; e[tot].c=z; e[tot].v=c; e[tot].nex=head[x]; head[x]=tot++;
    e[tot].to=x; e[tot].c=0; e[tot].v=-c; e[tot].nex=head[y]; head[y]=tot++;
}

bool spfa() {
	queue<int> q;
	memset(d, 0x3f, sizeof(d));
	memset(v, 0, sizeof(v));
	q.push(s); d[s] = 0; v[s] = 1;
	incf[s] = 1 << 30;
	while (q.size()) {
		int x = q.front(); v[x] = 0; q.pop();
		for (int i = head[x]; i!=-1; i = e[i].nex) {
			if (!e[i].c) continue;
			int y = e[i].to;
			if (d[y]>d[x] + e[i].v) {
				d[y] = d[x] + e[i].v;
				incf[y] = min(incf[x], e[i].c);
				pre[y] = i;
				if (!v[y]) v[y] = 1, q.push(y);
			}
		}
	}
	if (d[t] == 0x3f3f3f3f) return false;
	return true;
}

void update() {
	int x = t;
	while (x != s) {
		int i = pre[x];
		e[i].c -= incf[t];
		e[i^1].c += incf[t];
		x=e[i^1].to;
	}
	maxflow += incf[t];
	ans += d[t]*incf[t];
}

void init(){
    for(int i=1;i<=3*n;i++) head[i]=-1;
    tot=0; ans=0;
}

int main(){
    mst(dist,0x3f);
    scanf("%d%d%d",&n,&m,&k);
    for(int i=0;i<=n;i++) dist[i][i]=0;
    for(int i=1;i<=m;i++){
        int a,b,c; scanf("%d%d%d",&a,&b,&c);
        dist[a][b]=min(dist[a][b],c); dist[b][a]=min(c,dist[b][a]);;
    }
    for(int k=0;k<=n;k++)
        for(int i=0;i<=n;i++)
            for(int j=0;j<=n;j++){
                if(k<=max(i,j)) dist[i][j]=min(dist[i][j],dist[i][k]+dist[k][j]);
            }
    s=2*(n+1)+1; t=s+1;
    init();

    add(s,1,k,0);

    for(int i=1;i<=n;i++) add(s,i+1,1,0),add(n+1+i+1,t,1,0);

    for(int i=0;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(dist[i][j]!=0x3f3f3f3f)
                add(i+1,j+1+n+1,1,dist[i][j]);
        }
    }
    while (spfa()) update();
    printf(" %d %d\n",maxflow,ans);
}