HDU 4966 GGS-DDU (最小树形图-朱刘算法)

题目地址：HDU 4966
刚开始没看清总级别只有500这一条件，看成了每一个都是500。。然后建图思路就想歪了。。。。。后来才发现是总共只有500。。那么建图就很简单了。。把每个科目的每个等级都设为一个点，把所有的0等级设为同一个树根。然后把所有科目的高等级向低等级连一条权值为0的有向边，第一个作用是保证最后的最小树形图是所有点都可达，第二个作用是保证每节课的的所需等级，只要达到高等级，那么使低等级也符合。然后再对每节课连边即可。
代码如下：

#include <iostream>
#include <string.h>
#include <math.h>
#include <queue>
#include <algorithm>
#include <stdlib.h>
#include <map>
#include <set>
#include <stdio.h>
using namespace std;
#define LL __int64
#define pi acos(-1.0)
//#pragma comment(linker, "/STACK:1024000000")
const int mod=1e9+7;
const int INF=0x3f3f3f3f;
const double eqs=1e-9;
const int MAXN=40000+10;
int cnt, NV;
int in[510], color[510], id[510], pre[510], a[60], num[60];
struct node
{
        int u, v, w;
}edge[20000];
void add(int u, int v, int w)
{
        edge[cnt].v=v;
        edge[cnt].u=u;
        edge[cnt++].w=w;
}
int d_mst(int root)
{
        int ans=0, i, u, v;
        while(1){
                memset(in,INF,sizeof(in));
                for(i=0;i<cnt;i++){
                        u=edge[i].u;
                        v=edge[i].v;
                        if(u!=v&&edge[i].w<in[v]){
                                in[v]=edge[i].w;
                                pre[v]=u;
                        }
                }
                for(i=1;i<=NV;i++){
                        if(in[i]==INF&&i!=root) return-1;
                }
                int tot=0;
                memset(color,-1,sizeof(color));
                memset(id,-1,sizeof(id));
                in[root]=0;
                for(i=1;i<=NV;i++){
                        ans+=in[i];
                        u=i;
                        while(color[u]!=i&&id[u]==-1&&u!=root){
                                color[u]=i;
                                u=pre[u];
                        }
                        if(id[u]==-1&&u!=root){
                                id[u]=++tot;
                                for(v=pre[u];v!=u;v=pre[v]){
                                        id[v]=tot;
                                }
                        }
                }
                if(!tot) break;
                for(i=1;i<=NV;i++){
                        if(id[i]==-1) id[i]=++tot;
                }
                for(i=0;i<cnt;i++){
                        v=edge[i].v;
                        edge[i].u=id[edge[i].u];
                        edge[i].v=id[edge[i].v];
                        if(edge[i].u!=edge[i].v){
                                edge[i].w-=in[v];
                        }
                }
                NV=tot;
                root=id[root];
        }
        return ans;
}
int main()
{
        int n, m, i, j, k, c, d, l1, l2, w, ans;
        while(scanf("%d%d",&n,&m)!=EOF&&n+m){
                num[0]=1;
                for(i=1;i<=n;i++){
                        scanf("%d",&a[i]);
                        num[i]=num[i-1]+a[i];
                }
                cnt=0;
                NV=num[n];
                for(i=1;i<=n;i++){
                        for(j=1;j<=a[i];j++){
                                add(num[i-1]+j,1,0);
                                for(k=1;k<j;k++){
                                        add(num[i-1]+j,num[i-1]+k,0);
                                }
                        }
                }
                while(m--){
                        scanf("%d%d%d%d%d",&c,&l1,&d,&l2,&w);
                        add(l1?num[c-1]+l1:1,l2?num[d-1]+l2:1,w);
                }
                ans=d_mst(1);
                printf("%d\n",ans);
        }
        return 0;
}