bzoj 1016 最小生成树计数

    首先能发现一个规律，就是重构最小生成树的时候，一定不可能用一条权值较大的边和一条权值较小的边去替换他们中间的两条边。简而言之，就是只能权值相同的边相互替换。再进一步说，就是每种权值的边的数目是一定的。
　　题目上说值相同的边最多10条，那么我们可以DFS选哪些，然后用并查集来判断是否成环。这里要注意，我们先kruskal把每种边选几条求出来，再按kruskal的顺序dfs每一种，每dfs一种，我们就要把这些边连起来（也就是边dfs边kruskal），这样可以确保求出每种边的组合是合法的。最后把它们乘起来就行了。
　　看一些题解上说如果相同的边多了要用矩阵，表示不太懂......

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#define maxn 120
#define maxm 1200
#define ms 31011
using namespace std;
struct et
{
        int s,t,val;
}e[maxm*2];
int f[maxm],v[maxm],c[maxm],l[maxm],r[maxm],q[maxm];
bool vis[maxm];
int n,m,tot,num,now,cnt;
long long sum[maxm],ans;

int find(int i)
{
        if (!f[i]) return i;
        return find(f[i]);
}

void dfs(int col,int k,int ll)
{
        if (k>v[col]) {sum[col]++;return;}
        for (int i=ll+1;i<=r[col];i++)
        {
                int fx=find(e[i].s);
                int fy=find(e[i].t);
                if (fx!=fy) 
                {
            int tmp=f[fx];
                        f[fx]=e[i].t;
                        dfs(col,k+1,i);
                        f[fx]=tmp;
                }
        }
}

void add(int x,int y,int z)
{
        e[++tot].s=x; e[tot].t=y; e[tot].val=z;
}

bool cmp(et a,et b)
{
        return a.val<b.val;
}

int main()
{
        int x,y,z;
        scanf("%d %d",&n,&m);
        for (int i=1;i<=m;i++)
        {
                scanf("%d %d %d",&x,&y,&z);
                add(x,y,z);
        }
        sort(e+1,e+m+1,cmp);
        int j=0;
        for (int i=1;i<=m;i++)
        {
                if (e[i].val!=e[i-1].val) num++,l[num]=i;
                r[num]=i,c[i]=num;
        }//lisan
        for (int i=1;i<=m;i++)
        {
                int fx=find(e[i].s),fy=find(e[i].t);
                if (fx!=fy)
                {
                        f[fx]=e[i].t;
                        v[c[i]]++;
                }
        }//kruskal
        memset(f,0,sizeof(f));
    now=0;
        for (int i=1;i<=m;i++)
        {
        if (now!=c[i]) now=c[i],dfs(now,1,l[now]-1);
                int fx=find(e[i].s),fy=find(e[i].t);
                if (fx!=fy) f[fx]=e[i].t,cnt++;
        }//dfs
        if (cnt<n-1) ans=0; else ans=1;
        for (int i=1;i<=num;i++) if  (sum[i]) ans*=sum[i],ans%=ms;
        printf("%d\n",ans);
        return 0;
}