Homework of PE dp 正难则反

http://116.56.140.75:8000/JudgeOnline/problem.php?cid=1075&pid=3

Problem D: Homework of PE
Time Limit: 1 Sec Memory Limit: 128 MB
Submit: 117 Solved: 37
[Submit][Status][Web Board]
Description
At the begining of PE class, the math teacher comes to the class,
and says that the PE teacher has a sick.
So the math teacher gives Nobita a math homework.
The teacher says that s(n) indicates the number of primes in range [1,n].
And the teacher also says that for a permutation of 1-n, if the number x is in the xth position from left,
we say that x is in a correct position.
And then teacher wants to know.
For all permutation of 1-n, how many of them can satisfy that,
for all s(n) primes in range [1,n], there are exactly m of them are in the correct position.
And because the answer is so big, so Nobita needs to know the answer % 1000000007.
Input
The first line contains two integers n and m.
(1<=n<=200, 1<=m<=s(n)).
Output
Output the answer % 1000000007.
Sample Input
5 1
Sample Output
42

F[i][j]表示有i个数（由j个质数和（i-j）个质数）组成，j个质数不在自己位置上的方案数，

正难则反，
f[i][j]=jie[i]-sigma(f[i-k][j-k])(0

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int MAXN = 205;
const int P = (int)1e9 + 7;
LL jie[MAXN];
LL c[MAXN][MAXN];
LL f[MAXN][MAXN];
bool is_prime(int x){
    for(int i=2;i*i<=x;++i){
        if(x % i == 0)
            return false;
    }
    return true;
}
void init(const int n){
    jie[0] = 1;
    for(int i=1;i<=n;++i)
        jie[i] = jie[i - 1] * i % P;
    for(int i=0;i<=n;++i){
        for(int j=0;j<=i;++j){
            c[i][j] = j == 0 ? 1 : (c[i - 1][j] + c[i - 1][j - 1]) % P;
        }
    }
}
LL solve(const int n,const int m){
    init(n);
    int s = 0;
    for(int i=2;i<=n;++i){
        if(is_prime(i))
            ++s;
    }
    memset(f,0,sizeof f);
    for(int i=1;i<=n;++i){
        for(int j=0;j<=s;++j){
            LL tmp = 0;
            for(int k=1;k<=j;++k)
                (tmp += c[j][k] * f[i - k][j - k] % P) %= P;
            f[i][j] = (jie[i] - tmp + P) % P;
        }
    }
    return f[n - m][s - m] * c[s][m] % P;
}
int main(){
    int n,m;
    while(~scanf("%d %d",&n,&m)){
        printf("%d\n",(int)solve(n,m));
    }
    return 0;
}