zoj 3874 Permutation Graph (cdq分治+NTT)

   因为做做题学会了NTT ,比FFT的精度高了很多,收货很大。
/*
code by islands

y1=a[0]+a[1]x^1+a[2]x^2.....a[n]x^n
y2=b[0]+b[1]x^1+b[2]x^2.....b[m]x^m
y=y1*y2;  在O(nlgn)的复杂度求出y的系数

*/
#include<cstdio>
#include<cmath>
#include<queue>
#include<stack>
#include<string>
#include<cstring>
#include<iostream>
#include<map>
#include<vector>
#include<algorithm>
#include<stdlib.h>
#include<set>
#include<ctime>
#include<cmath>
using namespace std;
typedef long long LL;
const int mmax  = 1<<17;
const int mod = 786433;
const int g = 11;  //原根
LL quick_mod(LL a,LL b)
{
    LL ans=1;
    for(;b;b/=2)
    {
        if(b&1)
            ans=ans*a%mod;
        a=a*a%mod;
    }
    return ans;
}
int rev(int x,int r)  //蝴蝶操作
{
    int ans=0;
    for(int i=0; i<r; i++)
    {
        if(x&(1<<i))
        {
            ans+=1<<(r-i-1);
        }
    }
    return ans;
}
void NTT(int n,LL A[],int on) // 长度为N (2的次数)
{
    int r=0;
    for(;; r++)
    {
        if((1<<r)==n)
            break;
    }
    for(int i=0; i<n; i++)
    {
        int tmp=rev(i,r);
        if(i<tmp)
            swap(A[i],A[tmp]);
    }
    for(int s=1; s<=r; s++)
    {
        int m=1<<s;
        LL wn=quick_mod(g,(mod-1)/m);
        for(int k=0; k<n; k+=m)
        {
            LL  w=1;
            for(int j=0; j<m/2; j++)
            {
                LL t,u;
                t=w*(A[k+j+m/2]%mod)%mod;
                u=A[k+j]%mod;
                A[k+j]=(u+t)%mod;
                A[k+j+m/2]=((u-t)%mod+mod)%mod;
                w=w*wn%mod;
            }
        }
    }
    if(on==-1)
    {
        for(int i=1;i<n/2;i++)
            swap(A[i],A[n-i]);
        LL inv=quick_mod(n,mod-2);
        for(int i=0;i<n;i++)
            A[i]=A[i]%mod*inv%mod;
    }

}
LL A[mmax+10],B[mmax+10];
LL dp[mmax+10];
LL jie[mmax+10];
void cdq(int l,int r)
{

    if(l==r)
    {
        dp[l]=jie[l]-dp[l];
        dp[l]=(dp[l]%mod+mod)%mod;
        return ;
    }
    int mid=(l+r)>>1;
    cdq(l,mid);
    int n=r-l+1;
    for(int i=0,j=l;i<n;i++,j++)
    {
        if(j<=mid)
            A[i]=dp[l+i];
        else
            A[i]=0;
    }
    for(int i=0;i<n;i++)
        B[i]=jie[i+1];

    NTT(n,A,1);
    NTT(n,B,1);
    for(int i=0;i<n;i++)
        A[i]=A[i]*B[i]%mod;
    NTT(n,A,-1);
    for(int i=r,j=n-2;i>mid;i--,j--)
    {
        LL tmp=A[j];
        dp[i]+=tmp;
        dp[i]=(dp[i]%mod+mod)%mod;
    }
    cdq(mid+1,r);
}
void pre()
{
    jie[0]=1;
    for(int i=1;i<mmax;i++)
        jie[i]=jie[i-1]*i%mod;
    int n=1;
    while(n<=mmax) n<<=1;
    n/=2;
    memset(dp,0,sizeof dp);
    cdq(1,n);
}
int c[mmax+10];
int main()
{
    pre();
    int T;
    cin>>T;
    while(T--)
    {
        int n,m;
        scanf("%d %d",&n,&m);
        LL ans=1;
        while(m--)
        {
            int cnt;
            scanf("%d",&cnt);
            for(int i=0;i<cnt;i++)
                scanf("%d",&c[i]);
            sort(c,c+cnt);
            bool fg=1;
            for(int i=1;i<cnt;i++)
                if(c[i]!=c[i-1]+1)
                {
                    fg=0;
                    break;
                }
            if(!fg)
                ans=0;
            else
            {
                ans*=dp[cnt];
                ans%=mod;
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}

你可能感兴趣的:(分治,NTT算法)