CodeTON Round 8 (Div. 1 + Div. 2, Rated, Prizes!) E. Farm Game（博弈论+组合数学）

题目

FJ和FN（Farmer John和Farmer Nhoj）在农场玩游戏。

一条数轴，0和l+1两点有墙，[1,l]这l个点需要放2n头奶牛，其中n头属于FJ，n头属于FN

且两个人的奶牛，交错分布，即：

如果 𝑎1,𝑎2,...,𝑎𝑛 代表一个农场主的奶牛的位置，

而 𝑏1,𝑏2,...,𝑏𝑛 代表另一个农场主的奶牛的位置，

那么，有0<𝑎1<𝑏1<𝑎2<𝑏2<…<𝑎𝑛<𝑏𝑛<𝑙+1成立。

现在两个人开始下棋，

每次一个人选一个数字k(1<=k<=n)，并选择k头奶牛，

将它们同时向左或同时往右移动一个单位，

并且不能将自己的奶牛移动到墙上或其他农场主的奶牛上。

如果一位农夫无法移动任何奶牛，那么他就输了。

FJ先手，给定l(l<=1e6)和n(n<=l/2)，

求如果两个农夫都以最佳方式下棋，FJ可能获胜的棋局数，答案对998244353取模。

如果博弈无限期进行下去，此时认为没有农民获胜

实际t(t<=1e4)组样例，但保证suml不超过1e6

思路来源

菜菜园子群、propane

题解

考虑ai与bi之间的距离，即(a1,b1)、(a2,b2)，...

把奶牛贴上就有主动权了，后续不管对手怎么移动都可以效仿

首先第一点每个人只会向一个方向移动，

因为往反方向移动的话对手可以顺着该方向移动没有任何区别，

然后可以发现每次移动只可能有n个间距会改变，

所以相当于有n堆石子，每次操作可以选择若干个不空的堆取走一个

那么，FJ输当且仅当这n堆石子的距离都是奇数，

因为终态n堆奶牛贴上的时候距离都是奇数，

不然FJ可以选取所有距离为偶数堆的，将其变为奇数，从而必胜

然后将问题转化成，n个间隔，此外还有n+1个空隙，

每个间隔都是奇数，每个空隙都是>0的数，相加之和为n+1的方案数

a[1]+...+a[n]+b1+...+b[n+1]=l+1，

其中间隔a[i]均为奇数

每个间隔是奇数的方案数不好约束，

但是每个间隔是偶数的方案数是易求的，

用偶数限制并且同时除以2限制即可

所以，对每个间隔都加1，得到一个大于0的偶数，

令c[i]=a[i]+1=2*d[i]，

2(d[1]+...+d[n])+b1+...+b[n+1]=l+n+1

枚举n个间隔总共占了i（i为偶数）的位置，那么每个间隔分的就是i/2插挡板

剩下n+1个空隙，每个空隙分的是l+n+1-i插挡板

求出这部分必输的答案后，用C(l,2n)减去这部分答案

再乘以2即可，表示二人奶牛顺序可互换

代码

#include<iostream>
#include<cstdio>
#include<vector>
#include<map>
using namespace std;
#define rep(i,a,b) for(int i=(a);i<=(b);++i)
#define per(i,a,b) for(int i=(a);i>=(b);--i)
typedef long long ll;
typedef double db;
typedef pair<int,int> P;
#define fi first
#define se second
#define pb push_back
#define dbg(x) cerr<<(#x)<<":"<<x<<" ";
#define dbg2(x) cerr<<(#x)<<":"<<x<<endl;
#define SZ(a) (int)(a.size())
#define sci(a) scanf("%d",&(a))
#define pt(a) printf("%d",a);
#define pte(a) printf("%d\n",a)
#define ptlle(a) printf("%lld\n",a)
#define debug(...) fprintf(stderr, __VA_ARGS__)
using namespace std;
const int mod=998244353,N=1e6+10;
int Finv[N],fac[N],inv[N];
int t,l,n;
int modpow(int x,int n,int mod){
	int res=1;
	for(;n;x=1ll*x*x%mod,n>>=1)
	if(n&1)res=1ll*res*x%mod;
	return res;
}
void init(int n){ //n<N
    inv[1]=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)fac[i]=1ll*fac[i-1]*i%mod,Finv[i]=1ll*Finv[i-1]*inv[i]%mod;
	//Finv[n]=modpow(fac[n],mod-2,mod);
	//for(int i=n-1;i>=1;--i)Finv[i]=1ll*Finv[i+1]*(i+1)%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 x,int y){//x个分给y人插板
    return C(x-1,y-1);
}
int main(){
    init(N-5);
    sci(t);
    while(t--){
        sci(l),sci(n);
        int ans=0;
        for(int i=2;i<=l+1;i+=2){//所有的都是奇数 2(a1+...+a[n])+b1+...+b[n+1]=l+1
            ans=(ans+1ll*f(l+1+n-i,n+1)%mod*f(i/2,n)%mod)%mod;
        }
        //pte(ans);
        ans=((C(l,2*n)-ans)%mod+mod)%mod;
        ans=2ll*ans%mod;
        pte(ans);
    }
    return 0;
}

小飞侠的打表

#include <bits/stdc++.h>
using namespace std;

vector<vector<int>> gened;
map<vector<int>, int> f;

string tostr(vector<int>& x){
    std::string s = "{";
    for(int y: x) s += to_string(y) + ", ";
    s += "}";
    return s;
}

int dfs(vector<int>& cur, int l=0){
    //printf("cur==%s\n", tostr(cur).c_str());
    int n = cur.size();
    if(f.count(cur)){
        return f[cur];
    }
    int mn = 1e9;
    set<int> mex;
    for(int i = 1; i < 1<<n; ++i){
        mn = 1e9;
        for(int j = 0; j < n; ++j){
            if(i>>j&1){
                mn = min(mn, cur[j]);
            }
        }
        //printf("i==%d, mn==%d\n", i, mn);
        for(int x = 1; x <= min(1, mn); ++x){
            for(int j = 0; j < n; ++j){
                if(i>>j&1){
                    cur[j] -= x;
                }
            }
            int tmp = dfs(cur, l+1);
            //printf("l==%d, i==%d, x==%d, cur==%s, tmp==%d\n", l, i, x, tostr(cur).c_str(), tmp);
            mex.insert(tmp);
            for(int j = 0; j < n; ++j){
                if(i>>j&1){
                    cur[j] += x;
                }
            }
        }
    }
    int fuck = 0;
    while(mex.count(fuck)) {
        ++fuck;
    }
    return f[cur] = fuck;
}

void gen(vector<int>& t, int id, int lim){
    if(id >= t.size()){
        gened.push_back(t);
        return;
    }
    for(int i = 0; i <= lim; ++i){
        t[id] = i;
        gen(t, id+1, lim);
    }
}

int main(void){
    vector<int> t(4);
    gen(t, 0, 4);
    for(auto& x: gened){
        printf("dfs %s==%d\n", tostr(x).c_str(), dfs(x));
    }
    gened.clear();
    vector<int> t2(2);
    gen(t2, 0, 10);
    for(auto& x: gened){
        printf("dfs %s==%d\n", tostr(x).c_str(), dfs(x));
    }
    vector<int> tmp=  {1, 2};
    printf("hint dfs(1, 2)==%d\n", dfs(tmp));
}