BZOJ3444 最后的晚餐

在合法的情况下，最后肯定是所有人的喜欢关系组成了若干个链，然后这些链可以任意排列，长度大于2的链可以翻转，所以答案就是链数量的阶乘乘以2的长度大于1的链个数次方

不合法的情况：分叉了，有环

当然两个人互相喜欢这种情况是不算环的>_<

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 500010
#define MAXM 1010
#define ll long long
#define INF 1000000000
#define MOD 989381
#define eps 1e-8
int f[MAXN],siz[MAXN];
int d[MAXN];
int n,m;
int tot;
ll ans=1;
bool vis[MAXN];
map<int,bool>love[MAXN];
int fa(int x){
    return f[x]==x?x:f[x]=fa(f[x]);
}
int main(){
    int i,x,y;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++){
        f[i]=i;
        siz[i]=1;
    }
    for(i=1;i<=m;i++){
        scanf("%d%d",&x,&y);
        love[x][y]=1;
        if(fa(x)!=fa(y)){
			d[x]++;
	        d[y]++;
	        if(d[x]>2||d[y]>2){
	            printf("0\n");
	            return 0;
	        }
	        siz[fa(x)]+=siz[fa(y)];
	        f[fa(y)]=fa(x);
        }else if(!love[y][x]){
        	printf("0\n");
        	return 0;
        }
    }
    for(i=1;i<=n;i++){
        if(!vis[fa(i)]){
            vis[fa(i)]=1;
            if(siz[fa(i)]>=2){
                (ans*=2)%=MOD;
            }
            tot++;
        }
    }
    for(i=tot;i;i--){
        (ans*=i)%=MOD;
    }
    printf("%lld\n",ans);
    return 0;
}
 
/*
 
*/