Codeforces Round #733 (Div. 1 + Div. 2) F. Omkar and Akmar(组合数学+博弈思维题)

题目

Alice和Bob在一个圆格上玩游戏，格子共有1，...，n个，

格子i(i<n)与i+1相邻，n与1相邻，每个格子初始都是空的，

Alice先手，每次可以选一个空的格子，填A或者B，

特别地，填的时候，相邻的格子不能有相同的字母，否则无法填

一个玩家输，当且仅当该他填，但是他无法填了，

问二人都最优填的时候，终局有多少种，答案对1e9+7取模

思路来源

登录 - Luogu Spilopelia

[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;
}