【題解】JSOI2015染色問題

  好像這個容斥仍是明顯的。一共有三個要求,能夠用組合數先知足一個,再用容斥解決剩下的兩個維。(反正這題數據範圍這麼小,隨便亂搞均可以)。用 \(a[k][i]\) 表示使用 \(k\) 種顏色,至少有 \(i\) 列沒有染色的方案數,能夠容斥預處理獲得使用 \(k\) 種顏色染色使得每行每列均被染色的方案數。而後再容斥一下保證每種顏色都用上就能夠了。c++

#include <bits/stdc++.h>
using namespace std;
#define maxn 500
#define CNST 450
#define int long long
#define mod 1000000007
int n, m, K, ans, f[maxn];
int S[maxn], C[maxn][maxn];

int read()
{
    int x = 0, k = 1;
    char c; c = getchar();
    while(c < '0' || c > '9') { if(c == '-') k = -1; c = getchar(); }
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * k;
}

int Qpow(int x, int timer)
{
    int base = 1; if(timer < 0) return 1;
    for(; timer; timer >>= 1, x = x * x % mod)
        if(timer & 1) base = base * x % mod;
    return base;
}

void Up(int &x, int y) { x = (x + y) % mod; }
void Pre()
{
    for(int i = 0; i < CNST; i ++) C[i][0] = 1;
    for(int i = 1; i < CNST; i ++)
        for(int j = 1; j < CNST; j ++)
            Up(C[i][j], (C[i - 1][j - 1] + C[i - 1][j]) % mod);
}

int Get(int X)
{
    int ret = 0;
    for(int i = 0; i <= m; i ++) 
        S[i] = Qpow((Qpow(X + 1, m - i) - 1), n) % mod;
    for(int i = 0; i <= m; i ++)
        Up(ret, C[m][i] * ((i & 1) ? -S[i] : S[i]) % mod);
    return ret;
}

signed main()
{
    n = read(), m = read(), K = read();
    Pre(); for(int i = 0; i <= K; i ++) f[K - i] = Get(i);
    for(int i = 0; i <= K; i ++)
        Up(ans, C[K][i] * ((i & 1) ? -f[i] : f[i]) % mod);
    printf("%lld\n", (ans + mod) % mod);
    return 0;
}
相關文章
相關標籤/搜索