容斥原理。

求出f(m) ，f(m)指代至少有m个位置不合法的方案数。

怎么求？

注意到位置为id，权值为v ,不合法的情况，当且仅当 v = id+k或 v= id-k

因此，我们把每一个位置和权值抽象成点 ，不合法的情况之间连一条边，可以构成二分图。

借用大佬的图。

由此可知，当选了n条边，就恰好n个位置不合法，限制条件是：连的边不能相邻，

把二分图展开成k条链，进行dp。

还是借用大佬ez_lcw的图

由此总共有2n 个点 k 条链，链与链之间无边 互不干涉。

dp(i,j,pd)表示考虑到第i号点， 连了j条边，是否有连接i 到 i-1号点。

转移方程

则可得 

简单的乘法原理罢了。

#include<bits/stdc++.h>
using namespace std;
long long n,k;
long long fac[4010];
long long dp[4010][4010][2];
int mod = 1e9+7;
int pd[4999];//判断是否为链头 0 表示是头 头不能连接上一个 
int main(){
	freopen("neverk.in","r",stdin);
	freopen("neverk.out","w",stdout);
	cin>>n>>k;
	fac[0] = 1;
	for(int i = 1;i <= n;i++){
		fac[i] = (fac[i-1]*(long long)i)%mod;
	}
	int tot = 0;
	for(int i=1;i<=k;i++)
	{
		for(int t=0;t<2;t++)
		{
			for(int j=i;j<=n;j+=k)
			{
				tot++;
				if(i!=j) pd[tot]=1;
			}
		}
	}
	dp[0][0][0] = 1;
	for(int i = 1;i <= 2*n;i++){
		for(int j = 0;j <= n;j++){
			dp[i][j][0] = (dp[i-1][j][0] + dp[i-1][j][1])%mod;
			if(j&&pd[i]) dp[i][j][1] = dp[i-1][j-1][0];
		}
	}
	long long cnt = 0;
	for(int i = 0;i <= n;i++){
		//cout<<fac[n-i]<<" "<<dp[2*n][i][0]+dp[2*n][i][1]<<endl;
		long long t = fac[n-i]*(dp[2*n][i][0]+dp[2*n][i][1]);
		t%=mod;
		//cout<<t<<endl;
		if(i%2 == 0){
			cnt = (cnt +t)%mod;
		}
		else{
			cnt = (cnt - t + mod)%mod;
		}
	}
	cout<<cnt;
	return 0;
}