题目
Alice和Bob在一个圆格上玩游戏,格子共有1,...,n个,
格子i(i<n)与i+1相邻,n与1相邻,每个格子初始都是空的,
Alice先手,每次可以选一个空的格子,填A或者B,
特别地,填的时候,相邻的格子不能有相同的字母,否则无法填
一个玩家输,当且仅当该他填,但是他无法填了,
问二人都最优填的时候,终局有多少种,答案对1e9+7取模
思路来源
[CF1536F]Omkar and Akmar-CSDN博客
https://www.cnblogs.com/PPLPPL/p/14878872.html
题解
首先,应当推出结论,后手必胜,即二人的填数不存在策略
反证法可以证明,去掉空格后如果是奇数个位置,
一定存在两个颜色相同,并且这俩由于规则限制不能相邻,一定还可以填
所以,所有终态局面,都需要被统计,
枚举最后有k个空位,则去掉空位一定是ABABAB的串或者BABABA的串,且串长为偶数
可以把空位插入到这些字母里去,其中两个AB之间最多插一个空位,两个BA之间也是
可以看成是形如「球球空球空球」这样的一个长为n的成环的串,
那么由于两个空不会相邻,所以当n>=3时,每个空左右侧各有一个球,
统计「左空右球」的方案和统计「左球右空」的方案是等价的,
因为每个方案都能对应到一个「左空右球」上,一个「左球右空」上
所以,不妨认为是n个点的环上,选k对「左球右空」的方案
也可以认为是选k对相邻点匹配(每个点只能出现在一个匹配里)的方案,
拆环,
1. 如果第一个点是空,那么第n个点是球,长为n-2的链上选k-1对匹配的方案,
2. 如果第一个点不是空,长为n的链选k对匹配的方案,
然后可以给剩下的球涂色,
从第一个球开始涂A或涂B对应两种方案
k个空位对应n-k个字母,终态局面定下来之后,乘上俩人填n-k个字母的顺序,
即为所求
代码1
#include <iostream>
//#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll modpow(ll x,ll n,ll mod){ll res=1;for(;n;n>>=1,x=x*x%mod)if(n&1)res=res*x%mod;return res;}
const int mod=1e9+7;
const int N=1e6+10;
bool ok[N];
int prime[N],phi[N],two[N],cnt;
int n,k,Finv[N],fac[N],inv[N];
void init(int n){
inv[1]=1;two[0]=1;
for(int i=2;i<=n;++i)inv[i]=(1ll*(mod-mod/i)*inv[mod%i])%mod;
fac[0]=Finv[0]=1;
for(int i=1;i<=n;++i)two[i]=2ll*two[i-1]%mod,fac[i]=1ll*fac[i-1]*i%mod,Finv[i]=1ll*Finv[i-1]*inv[i]%mod;
}
int C(int n,int m){
if(m<0||m>n)return 0;
return 1ll*fac[n]*Finv[n-m]%mod*Finv[m]%mod;
}
int cal(int n){
int res=0;
for(int k=0;k<=n/2;++k){//枚举有多少个空位
if((n-k)%2==0){
int v=1ll*(C(n-k,k)+C(n-k-1,k-1))%mod%mod;
res=(res+2ll*v*fac[n-k])%mod;
}
}
return res;
}
int main(){
init(N-5);
scanf("%d",&n);
printf("%d\n",cal(n));
return 0;
}
代码2
#include <iostream>
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll modpow(ll x,ll n,ll mod){ll res=1;for(;n;n>>=1,x=x*x%mod)if(n&1)res=res*x%mod;return res;}
const int mod=1e9+7;
const int N=1e6+10;
bool ok[N];
int prime[N],phi[N],two[N],cnt;
int n,k,Finv[N],fac[N],inv[N];
void init(int n){
inv[1]=1;two[0]=1;
for(int i=2;i<=n;++i)inv[i]=(1ll*(mod-mod/i)*inv[mod%i])%mod;
fac[0]=Finv[0]=1;
for(int i=1;i<=n;++i)two[i]=2ll*two[i-1]%mod,fac[i]=1ll*fac[i-1]*i%mod,Finv[i]=1ll*Finv[i-1]*inv[i]%mod;
}
int C(int n,int m){
if(m<0||m>n)return 0;
return 1ll*fac[n]*Finv[n-m]%mod*Finv[m]%mod;
}
int f(int n,int k){
return C(n-k+1,k);
}
int cal(int n){
int res=0;
for(int k=0;k<=n/2;++k){//枚举有多少个空位
if((n-k)%2==0){
int v=1ll*(f(n,k)-f(n-4,k-2)+mod)%mod%mod;
res=(res+2ll*v*fac[n-k])%mod;
}
}
return res;
}
int main(){
init(N-5);
scanf("%d",&n);
printf("%d\n",cal(n));
return 0;
}