模是1e8+7啊卧槽真tm蛋疼……
因为题里要求循环同构的算同一种,所以我们考虑一下burnside,置换一共有n个,分别是往后窜1~n个位置
对于窜x个位置的置换,一共有gcd(n,x)个循环节,且第i个位置属于第i%(gcd(n,x))+1个循环节
那么对于一个窜x个位置的置换,合法的不动点数量就是长度为gcd(n,x)的环,不考虑循环同构,没有超过k个女生坐在一起的方案数
如果gcd(n,x)<=k,那么如果k<n,合法方案数为2^gcd(n,x)-1,如果k>=n,合法方案数为2^gcd(n,x)
否则
用f[i][j]表示长度为i的线段(不是环),最后有且仅有j个女生,没有超过k个女生坐在一起的方案数
F[i][j]表示长度为i的线段,最后有且仅有j个女生,且第一个是男生,没有超过k个女生坐在一起的方案数
简单n^2递推
那么长度为i的环的方案就是长度为i的线段的方案数减去两边接起来超过k个了的方案数
枚举两边的女生个数和简单算一下即可
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<queue>
#include<map>
#include<bitset>
#include<stack>
#include<vector>
#include<set>
using namespace std;
#define MAXN 2010
#define MAXM 1010
#define INF 1000000000
#define MOD 100000007
#define ll long long
#define eps 1e-8
int n,k;
int mi[MAXN];
int f[MAXN][MAXN],F[MAXN][MAXN],S[MAXN],s[MAXN],g[MAXN];
int gcd(int x,int y){
return !y?x:gcd(y,x%y);
}
ll MI(ll x,ll y){
ll re=1;
while(y){
if(y&1){
(re*=x)%=MOD;
}
(x*=x)%=MOD;
y>>=1;
}
return re;
}
int main(){
int i,j;
int tmp;
mi[0]=1;
for(i=1;i<=2000;i++){
mi[i]=(mi[i-1]<<1)%MOD;
}
scanf("%d",&tmp);
while(tmp--){
scanf("%d%d",&n,&k);
if(k>n){
k=n;
}
f[0][0]=s[0]=1;
F[0][0]=S[0]=1;
F[1][0]=S[1]=1;
for(i=2;i<=n;i++){
F[i][0]=S[i-1];
S[i]=F[i][0];
for(j=1;j<=k;j++){
F[i][j]=F[i-1][j-1];
(S[i]+=F[i][j])%=MOD;
}
}
for(i=1;i<=n;i++){
f[i][0]=s[i-1];
s[i]=f[i][0];
for(j=1;j<=k;j++){
f[i][j]=f[i-1][j-1];
(s[i]+=f[i][j])%=MOD;
}
g[i]=s[i];
for(j=k+1;j<=min(2*k,i-1);j++){
(g[i]+=MOD-(ll)F[i-j][0]*max(0,min(k,j-1)-max(1,j-k)+1)%MOD)%=MOD;
}
}
ll ans=0;
for(i=1;i<=n;i++){
int l=gcd(i,n);
if(l<=k){
if(k<n){
(ans+=mi[l]-1+MOD)%=MOD;
}else{
(ans+=mi[l])%=MOD;
}
}else{
(ans+=g[l])%=MOD;
}
}
(ans*=MI(n,MOD-2))%=MOD;
printf("%lld\n",ans);
}
return 0;
}
/*
1 9 8
*/