图论：朱刘算法

有向图中点点连通边权之和最小

算法过程不研究了，以后能看懂再说。。

直接贴一道以前写过的题

Openjudge的题面是地震之后，实则为一道POJ题目

#include<cstdio>  
#include<iostream>  
#include<cstring>  
#include<cmath>  
using namespace std;  
const int maxn = 105;  
const int INF = 0x7fffffff;  
template<typename Type>  
struct Directed_MST
{  
    void init(int _n)
    {  
        n=_n;   
        ans = 0;  
        memset(vis, 0, sizeof(vis));  
        memset(inc, 0, sizeof(inc));  
        for(int i=0; i<=n; ++i)
        {  
            w[i][i] = INF;  
            for(int j=i+1; j<=n; ++j)  
                w[i][j]=w[j][i]=INF;  
        }  
    }  
    void insert(int u, int v, Type _w)
    {  
        if(w[u][v]>_w) w[u][v] = _w;  
    }  
    Type directed_mst(int u)
    {    
        dfs(u);  
        for(int i=1; i<=n; ++i)  
            if(!vis[i]) { return -1; }  
        memset(vis, 0, sizeof(vis));  
        while(true)
        {    
            for(int i=1; i<=n; ++i)if(i!=u&&!inc[i])
            {  
                w[i][i]=INF, pre[i] = i;  
                for(int j=1; j<=n; ++j)if(!inc[j] && w[j][i]<w[pre[i]][i])
                {  
                    pre[i] = j;  
                }  
            }   
            int i;  
            for(i=1; i<=n; ++i)if(i!=u&&!inc[i])
            {  
                int j=i, cnt=0;  
                while(j!=u && pre[j]!=i && cnt<=n) j=pre[j], ++cnt;  
                if(j==u || cnt>n) continue;  
                break;  
            }                 
            if(i>n)
            {    
                for(int i=1; i<=n; ++i)if(i!=u && !inc[i]) ans+=w[pre[i]][i];  
                return ans;  
            }     
            int j=i;  
            memset(vis, 0, sizeof(vis));  
            do{  
                ans += w[pre[j]][j], j=pre[j], vis[j]=inc[j]=true;  
            }while(j!=i);  
            inc[i] = false; 
            for(int k=1; k<=n; ++k)if(vis[k])
            {
                for(int j=1; j<=n; ++j)if(!vis[j])
                {
                    if(w[i][j] > w[k][j]) w[i][j] = w[k][j];  
                    if(w[j][k]<INF && w[j][k]-w[pre[k]][k] < w[j][i])  
                        w[j][i] = w[j][k] - w[pre[k]][k];  
                }    
            }  
        }  
        return ans;        
    }  
    void dfs(int u)
    {  
        vis[u] = true;  
        for(int i=1; i<=n; ++i)if(!vis[i]&&w[u][i]<INF)
        {  
            dfs(i);  
        }  
    }   
    Type ans;
    int n; 
    int pre[maxn];
    bool vis[maxn];
    bool inc[maxn];
    Type w[maxn][maxn];
};  
Directed_MST<double>G;  
double x[maxn],y[maxn];  
double getDist(double x1,double y1,double x2,double y2
)
{  
    return sqrt(pow(x1-x2,2)+pow(y1-y2,2));  
}  
  
int main()
{  
    int n,m;  
    while(~scanf("%d%d",&n,&m))
    {  
        G.init(n);  
        for(int i=1; i<=n; ++i)  
            scanf("%lf%lf",&x[i],&y[i]);  
        for(int i=0; i<m; ++i)
        {  
            int a,b;  
            scanf("%d%d",&a,&b);  
            if(a==b)continue;  
            G.insert(a,b,getDist(x[a],y[a],x[b],y[b]));  
        }  
        double ans = G.directed_mst(1);  
        if(ans < 0) puts("NO");  
        else printf("%.2f
", ans);  
    }  
    return 0;  
} 

裸的最小树形图