PossibleOrders TopCoder - 1643 hqg_ac

并查集维护，类似第二类斯特林数

$DP[i][j]$表示i个元素构成j个集合的方案数

$dp[i][j]=(dp[i-1][j]+d[i-1][j-1])*j$

#include <bits/stdc++.h>
using namespace std ;
typedef long long ll ;
struct PossibleOrders{
    int f[110],l[110],r[110],x,y ;
    int find(int x){return f[x]==-1?x:f[x]=find(f[x]);}
    ll work(int n){
        ll dp[20][20] ;
        memset(dp,0,sizeof(dp)) ;
        dp[0][0]=1;
        for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++) 
        dp[i][j]=(dp[i-1][j]+dp[i-1][j-1])*j ;
        ll ans=0;
        for (int i=1;i<=n;i++) ans+=dp[n][i] ;
        return ans ;
    } 
    int change(string s){
        int ans=0,k=0;
        while (k<s.length()) ans=ans*10+s[k++]-'0' ;
        return ans ;
    }
    ll howMany(int n,vector<string> a){
        for (int i=0;i<a.size();i++){
            int t=a[i].find("=") ;
            l[i]=change(a[i].substr(0,t)) ;
            r[i]=change(a[i].substr(t+1,a[i].length()-t-1)) ;
        }
        for (int i=0;i<n;i++) f[i]=-1 ;
        for (int i=0;i<a.size();i++) 
        if ((x=find(l[i]))!=(y=find(r[i]))) f[x]=y ;
        int m=0 ;
        for (int i=0;i<n;i++) if (f[i]<0) m++ ;
        return work(m) ; 
    }
};