2023百度之星第三场【染色游戏】题解

题目：

小度在公园玩一个染色游戏，染色板为一个长为 n，宽为 m 的长方形网格，一开始它们的颜色都是白色。
小度的颜料可以将其中的 k 个的格子染成黑色，颜料需要用完，且不能重复染色。
最终的要求是任意相邻两行或任意相邻两列要么保证完全一致，要么完全不一致。

完全一致指相邻行/列中相邻的格子要么同为白色，要么同为黑色。
完全不一致指相邻行/列中相邻的格子一个为白色，一个为黑色。

请计算有多少种染色方案。

这是比赛第二题，赛时最后半小时第一题想不出来了，来看这道题，有了一点想法，但是后面卡住了，赛后看了别人题解发现自己题意理解错了，这题算是能补的范围内，故在此叙述一下补题的成果吧。

对于题意的理解，可以发现以下的性质，对于已经涂好色的一行，他的下一行要么和它相同，要么和它相反，所以颜色都是由前一行确定，第二个最重要的性质是只要行满足了这个条件，列也会满足这个条件，所以我们只要考虑行。我赛时理解错的地方在于，我以为必须全部都相同，或者必须交错地都相反。其实可以部分行相同，部分行相反。而我那种理解有一种很严重的问题就是无解的情况很多，全部相同，全部相反而言都对k有严格要求，才能有解。正确的题意的话，可以适当地选择相同相反的行数，使得最终染色的总数为k。

因此我们不妨枚举第一行中染色的块数为i，那么第一行的染色方案数就是C(m,i)，在m个列中选i个列涂色。然后为了让染色总和为k，假设n行中总共有y个行和第一行染色方案相同，则有(n-y)个行的染色方案和第一行不同，不同的行染色的块数就是(m-i)。我们就可以得到以下的方程：

y*i+(n-y)*(m-i)==k

我们可以从中判断是否存在这样的y，如果有解，则答案可以加上C(m,i)*C(n,y)。

这里有个重要的理解地方，因为我们的i的意义是枚举第一行中染色的块，但是C(n,y)是任意选y行，可能第一行并没有被选中。所以我们需要修改一下想法。

对于最终的染色好的成果而言，每一行的染色成果，要么是超半的块被染，要么不到一半的块被选，所以我们枚举这两种方案中较少染色块的块数，另外超半的染色方案也被对应确立下来。所以枚举i从0到n/2，如果n是偶数，再去特判i==n/2时的答案，累加到答案中。

那么最终的代码就如下所示了：

//#pragma GCC optimize(2)
#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define pii pair<int,int>
#define pll pair<long long,long long>
#define int long long
using namespace std;
int read(){
    int x=0,f=0,ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=1;ch=getchar();}
    while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
    return f?-x:x;
}
const int maxn=1e7+10;//需要的最大组合数n，m的范围
const int mod=998244353;
ll fac[maxn], ifac[maxn];
inline ll ksm(ll a,ll b){
    ll res=1;
    for(;b;b>>=1,a=a*a%mod){
        if(b&1) res=res*a%mod;
    }
    return res;
}
inline ll C(int n, int m){//m是小的数字，n是大的
    if (n < m || m < 0)return 0;
    return 1ll*fac[n] * ifac[m] % mod * ifac[n - m] % mod;
}

inline ll A(int n, int m){//m是小的数字，n是大的
    if (n < m || m < 0)return 0;
    return fac[n] * ifac[n - m] % mod;
}
void init(){
    fac[0] = 1;//阶乘预处理
    for (int i = 1; i <= maxn - 5; i++)
        fac[i] = fac[i - 1] * i % mod;
    ifac[maxn - 5] = ksm(fac[maxn - 5], mod - 2);
    for (int i = maxn - 5; i; i--)
        ifac[i - 1] = ifac[i] * i % mod;
}
signed  main(){
    init();
    int n=read(),m=read(),k=read();
    int ans=0;
    for(int i=0;i<m/2;i++){//枚举较小的染色格数中选取i个格子染色
    //y*i+(m-i)*(n-y)==k
        int t=(m-i)*n-k,b=m-2*i;
        if(t%b) continue;
        int y=t/b;
        ans+=C(m,i)*C(n,y);
        ans%=mod;
    }
    if(m%2==0&&(n*m/2==k)){
        ans+=C(m,m/2)*ksm(2,n-1);
        ans%=mod;
    }
    cout<<ans<<endl;
    system("pause");
    return 0;
}