CodeForces 285 D.Permutation Sum（dfs）

Description

$a$和$b$为两个$1$~$n$的排列，问有多少种可能使得序列$x_i=(a_1-1+b_i-1)\ mod \ n+1$也是一个排列

Input

一个整数$n(1\le n\le 16)$

Output

输出方案数，结果模$10^9+7$

Sample Input

3

Sample Output

18

Solution

首先固定$a_i=i$，找到合法的$b_i$方案数后答案乘$n!$即可，因为$n$很小所以直接暴搜合法$b$序列本地打表即可，每次用$B$($16$位$01$)表示$b$序列没用的数字，以及用$C$(也是$16$位$01$)表示 c <script type="math/tex" id="MathJax-Element-20">c</script>序列已经有的数字，以此剪枝会快很多

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
#define mod 1000000007
int ans[]={0,1,0,18,0,1800,0,670320,0,734832000,0,890786230,0,695720788,0,150347555,0};
int fact[22];
int num,m;
void dfs(int pos,int B,int C)
{
    if(pos==m+1)
    {
        num++;
        if(num>=mod)num-=mod;
        return ;
    }
    for(int i=1;i<=m;i++)
        if(!((1<<i)&B))
        {
            int j=(i-1+pos-1)%m+1;
            if(!((1<<j)&C))
            dfs(pos+1,B^(1<<i),C^(1<<j));
        }
}
void init()
{
    fact[0]=1;
    for(int i=1;i<=16;i++)fact[i]=(ll)i*fact[i-1]%mod;
    for(int i=1;i<=16;i++)
    {
        m=i,num=0;
        dfs(1,0,0);
        ans[i]=(ll)num*fact[i]%mod;
        printf("%d,",ans[i]);
    }
}
int main()
{
    //init();
    int n;
    while(~scanf("%d",&n))printf("%d\n",ans[n]);   
    return 0;
}