POJ 1639度限制生成树

题目就是给出一些英文名称， 并且要求源点为park的入度不能超过k
给出思路：
1.首先抹去park， 在剩下的连通分量中求出最小生成树加入答案， 此时的度是最小的， 即如果现在的度大于k， 则无解（题目好像没这个情况）
2.由于我们一开始删去了和park相邻的每一条边， 所以我们要加上每个联通子图到park的最小边， 但这样并不是最好的， 我们要把度加到k， 这样会使答案更优（至于具体的操作， 留一点思考的空间， 也可以看我的代码）

#include <map>#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>using namespace std;const int N = 150;const int oo = 0x3f3f3f;#define FILL(a, b) memset(a, b, sizeof(a))#define rep(i, s, t) for(int i = s; i <= t; ++i)map<string, int> id;int insert(string s) {    if(!id.count(s)) id[s] = id.size();    return id[s];}struct Edge {    int u, v, w;}p[N];struct Part_Prim {    int G[N][N], n, m, limit, fa[N];    int dis[N], vis[N], cnt, used[N][N];    void init() {        id.clear(); id["Park"] = 1; cnt = 0;        rep(i, 1, N-1) rep(j, 1, N-1) G[i][j] = oo;        FILL(vis, 0), FILL(used, 0), FILL(p, 0);    }    int prim(int s, int ret = 0, int add = oo) {        int to_root = 0; ++cnt;        rep(i, 0, n) dis[i] = oo;        dis[s] = 0;        rep(i, 1, n-1) {            int k = 0;            rep(j, 2, n) if(dis[k] > dis[j] && !vis[j]) k = j;            if(!k) break;            ret += dis[k]; vis[k] = cnt;            if(add > G[k][1]) to_root = k, add = G[k][1];            rep(j, 2, n)                 if(G[k][j] < dis[j] && !vis[j])                     dis[j] = G[k][j], fa[j] = k;            if(fa[k]) used[k][fa[k]] = used[fa[k]][k] = 1;        }        used[1][to_root] = used[to_root][1] = 1;        return ret + add;    }    void dfs(int u, int f) {        rep(v, 2, n)            if(used[u][v] && v ^ f) {                p[v] = p[u];                if(G[u][v] ^ oo && G[u][v] > p[v].w)                    p[v] = (Edge) {u, v, G[u][v]};                dfs(v, u);            }    }    int parking(int ret = 0, int t = 0) {        rep(i, 2, n) if(!vis[i]) ret += prim(i);//      cout << ret << endl;        while(cnt++ < limit) {            dfs(1, 0);            rep(i, 2, n)                 if(G[1][i] ^ oo && !used[1][i] && p[i].w-G[1][i] > p[t].w-G[1][t])                     t = i;            ret -= p[t].w-G[1][t]; //          cout << p[t].w << "!" << endl;            used[1][t] = used[t][1] = 1;            used[p[t].u][p[t].v] = used[p[t].v][p[t].u] = 0;        }        return ret;    }    void solve(int _ = 0) {        while(~scanf("%d", &_)) {            init();            while(_--) {                string s1, s2; int w, u, v;                cin >> s1 >> s2 >> w;                u = insert(s1), v = insert(s2);                G[u][v] = G[v][u] = w;            }            n = id.size();            scanf("%d", &limit);            printf("Total miles driven: %d\n", parking());        }    }}WORK;int main() {#ifndef ONLINE_JUDGE    freopen("IP.in", "r", stdin);    freopen("OP.out", "w", stdout);#endif    WORK.solve();    return 0;}