The Review Plan I ZOJ 3687

The Review Plan I  

题意：  你要用n 天复习  n 门课， 每天一门，但是有限制条件，某一天不能复习某几门课。

禁位排列的模板题，看了好长时间别人的博客，看懂了棋子多项式，但是代码看不懂（orz）

棋子多项式与禁位排列：点击打开链接         点击打开链接

定理：

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define MOD 55566677
#define N 55
using namespace std;


typedef long long ll;
ll f[N],res;
int n,m;
int a[N][N];
int vis1[N],vis2[N];
int ban[N][N];
void init()
{
    f[0]=1;
    for(int i=1;i<N;i++)
        f[i]=(f[i-1]*i)%MOD;
}
void dfs(int k,int num)
{
    if(k>m)
    {
        if(num&1)  res=((res-f[n-num])%MOD+MOD)%MOD;
        else   res=(res+f[n-num])%MOD;
        cout<<" res "<<res<<endl;
        return ;
    }
    dfs(k+1,num);
    if((!vis1[a[k][0]])&&(!vis2[a[k][1]]))
    {
        vis1[a[k][0]]=vis2[a[k][1]]=1;
        dfs(k+1,num+1);
        vis1[a[k][0]]=vis2[a[k][1]]=0;
    }

}
int main()
{
    init();
    int i;
    while(~scanf("%d%d",&n,&m))
    {
        memset(vis1,0,sizeof vis1);
        memset(vis2,0,sizeof vis2);
        memset(ban,0,sizeof ban);
        for(i=1;i<=m;i++)
        {
            scanf("%d%d",&a[i][0],&a[i][1]);
            if(ban[a[i][0]][a[i][1]])
            {
                i--;
                m--;
            }
            else
                ban[a[i][0]][a[i][1]]=1;
        }
        res=0;
        dfs(1,0);
        res=(res%MOD+MOD)%MOD;
        printf("%lld\n",res);
    }
    return 0;
}