[JSOI2015]染色问题

Posted Hs-black

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[JSOI2015]染色问题相关的知识,希望对你有一定的参考价值。

一道比较简单的容斥题(?

可以容斥一下, 设(i, j ,k)分别代表(i)行有颜色, (j)行有颜色, 出现(k)种颜色的方案数
[ ans = sum_{i=0}^nsum_{j=0}^msum_{k=0}^c(-1)^{(i+j+k+n+m+c)}C_n^iC_m^jC_c^k(k+1)^{(n*m)} ]
k+1是因为一个格子可以不填颜色, 反手一个预处理组合数, O(n^3)稳过

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
using namespace std;
const int P = 1e9+7;
const int N = 505;
ll C[N][N], fpw[N*N];
int n, m, c;
ll ans;
int main() {
    cin >> n >> m >> c;
    fpw[0] = C[0][0] = 1;
    for (int i = 1;i <= 400; i++) {
        C[i][0] = 1;
        for (int j = 1;j <= i; j++) {
            C[i][j] = C[i-1][j] + C[i-1][j-1];
            if (C[i][j] >= P) C[i][j] -= P;
        }
    }
    int p = (n ^ m ^ c) & 1;
    for (int k = 0;k <= c; k++) {
        for (int i = 1;i <= n * m; i++) fpw[i] = fpw[i-1] * (k+1) % P;
        for (int i = 0;i <= n; i++) {
            for (int j = 0;j <= m; j++) {
                ll res = C[n][i] * C[m][j] % P * C[c][k] % P * fpw[i*j] % P;
                if (((i ^ j ^ k) & 1) ^ p) ans -= res;
                else ans += res;
            }
        }
    }
    cout << (ans % P + P) % P << endl;
    return 0;
}

以上是关于[JSOI2015]染色问题的主要内容,如果未能解决你的问题,请参考以下文章

bzoj4487[Jsoi2015]染色问题 容斥原理

[Jsoi2015]染色问题

bzoj 4487: [Jsoi2015]染色问题

[JSOI2015]染色问题

2015 JSOI冬令营训练 彩色格子 题解

7.多关联双亲染色体片段代换系(CSSL)群体的代谢组分析