CodeForces 997C Sky Full of Stars

本文参考自: 原文地址

题目：点击打开链接
题意：给n*n个格子，每个格子可以填3种颜色。问有多少种填色方案，至少有一列或一行是同色的。
分析：首先反过来考虑，找不幸运的组合，因为这样便于容斥。 
要把列和行分开考虑。 
首先排除掉列上出现相同的情况，也就是(3^n−3)^n种，这个也就是容斥里开始时最大的范围，显然还包括了一些行上相同的情况，所以下面开始容斥。 
在列上不出现相同情况下，行出现相同的情况要分两类讨论， 
1. 那些整行出现的同一颜色的行的颜色完全相同，比如有三行出现整行颜色一样，并且都是绿，这种情况列上的颜色不能和他们相同，有3∗C(i,n)∗(3^(n−i)−1)^n 种。 
2. 那些整行出现的同一颜色的行的颜色不完全相同，比如有三行出现整行颜色一样，但他们分别是红黄蓝，这种情况，列上就可以随便选了，有(3^i−3)∗C(i,n)∗3^(n−i)n。 
最终使用容斥。

ans=3n2−(3n−3)n−Σni=1(−1)iCin(3∗(3n−i−1)n+(3i−3)∗3(n−i)n)

代码：

#pragma comment(linker, "/STACK:102400000,102400000")///手动扩栈
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cstdio>
#include<bitset>
#include<vector>
#include<cmath>
#include<ctime>
#include<stack>
#include<queue>
#include<deque>
#include<list>
#include<set>
#include<map>
using namespace std;
#define debug test
#define mst(ss,b) memset((ss),(b),sizeof(ss))
#define ll long long
#define ull unsigned long long
#define pb push_back
#define mp make_pair
#define inf 0x3f3f3f3f
#define eps 1e-10
#define MOD 998244353
#define PI acos(-1.0)
const int N = 1e6+10;

ll gcd(ll p,ll q){return q==0?p:gcd(q,p%q);}
int to[4][2]={{-1,0},{1,0},{0,-1},{0,1}};

ll n,ans;

ll qp(ll a,ll b) {
    ll s=1;
    while(b) {
        if(b&1) s=s*a%MOD;
        a=a*a%MOD;
        b>>=1;
    }
    return s;
}

int main() {
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    cin>>n;
    ans=(qp(3,n*n)-qp(qp(3,n)-3,n)+MOD)%MOD;
    ll xs=-1,c1=1;
    for(int i=1;i<=n;i++) {
        xs*=-1;
        c1=c1*(n-i+1)%MOD;
        c1=c1*qp(i,MOD-2)%MOD;
        ans += xs*(3*c1%MOD*qp(qp(3,n-i)-1,n)%MOD + c1*(qp(3,i)-3)%MOD*qp(3,(n-i)*n)%MOD)%MOD;
    }
    ans = (ans%MOD+MOD)%MOD;
    cout<<ans<<endl;
    return 0;
}

参考博客：https://blog.csdn.net/z631681297/article/details/80889816

ans=3n2−(3n−3)n−Σni=1(−1)iCin(3∗(3n−i−1)n+(3i−3)∗3(n−i)n)