dp_Zju_3688_The Review Plan II_不愿若菜

• F[n]:n个限位排列的方案数;
  g[n]:n个限位排列中有一项映射到自己，剩下所有项依然遵循限制所获得的方案数。//不难推得，k映射到k和k映射到k+1所得到的g[n]是相等的。
• 则：F[n] = F[n-1]*(n-3)+[2*(n-3)+1]*g[n-1]+g[n]+g[n-1]+g[n-2];
•     g[n] = F[n-1]+g[n-1];
• 考虑n的情况：1不能映射到{1,2}，...，n-1不能映射到{n-1,n}，n不能映射到{n,1}.
• 首先考虑一点：1不能映射到{n,2}，...，n-1不能映射到{n-1,n}的方案数也是F[n-1](这里只不过是把值域里的1换成n而已)。
• F[n]分两类考虑

    1. 考虑1不能映射到n的情况

• 那么对于n，我们先分配1...n-1，这n-1个映射，并且按照：1不能映射到{n,2}，...，n-1不能映射到{n-1,n}的规则映射，那么有方案数F[n-1]，然后n没有映射到任何数，并且1没有"被映射"到，这时候我们把k及k的象取出，改为 ,其中（k!=1,f(k)!=k），对于F[n-1]的每个方案，我们可以取(n-3)个这样的k，于是得到F[n-1]*(n-3);
• 注意到上述策略是不完全的。因为F[n]中可以有^和^的情况，但f(k)!=k且f(k)!=k+1。所以上述策略无法生成^这样的映射。对于这种情况，我们需要在1...n-1中安排一个映射使其违反限制。那么这个k可以在2...n-1中选择(2...n-2可以有两种违反限制的方案，n-1只有一种）；

    2. 然后我们考虑1映射到n的情况，等价于k映射到k-1的情况。这时候，有k-2不能映射到{k-2},k-1不能映射到{k},其他限制照旧。不难推出，这时候的方案数为 g[n]+g[n-1]+g[n-2];

• 然后对于g[n]的递推就更简单了。如果前面的都看得懂，这里就更不用说了:g[n] = f[n-1]+g[n-1];

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#define clr(A) memset(A,0,sizeof(A))
using namespace std;
typedef long long LL;
const int MM = 100005;
const int Mod = 1000000007;
int F[MM],g[MM];
int main()
{
    //freopen("Pr_temp.in","r",stdin);
    F[3] = 1,g[4] = 1;
    for(int i = 4;i<=100000;i++){
        F[i] = ((LL)F[i-1]*(i-2)+(2*i-3)*(LL)g[i-1]+g[i-2])%Mod;
        g[i] = (F[i-1]+g[i-1]) % Mod;
    }
    int x;
    while(~scanf("%d",&x))
    printf("%d\n",F[x]);
    return 0;
}