bzoj 1004
Pros
給若干個置換, 和三種顏色, 問有多少種本質不同的染色方案?
Sol
Burnside:
首先由 \(burnside\) 引理
\[Ans = \frac{\sum_{\:i = 1}^{\:m}\;f_i}{m}\]
其中 \(f_i\) 表示第 \(i\) 個置換的不動點的個數.
不動點的個數 就是在該置換下只有一種同構體的方案數.
置換的循環分解:
然后,前面 lyp 學長在這里講,
關于一個置換能夠分解成若干個環,
\(example:\)
\[ \begin{equation} \left( \begin{array}{ccccc} 1 & 2 & 3 & 4 & 5\\ 3 & 2 & 5 & 1 & 4\\ \end{array} \right) \end{equation} \]
包含環:\(1 \rightarrow 3 \rightarrow 5 \rightarrow 4 \rightarrow 1\) 和 \(2 \rightarrow 2\)
對于任意一個環, 作用在上面的一次置換能夠讓這個環上面的元素循環移位一格.
不動點的條件
那么要讓這個環能夠成為不動點的充要條件就是,每個環上的顏色都是相同的.
然后,我們就能夠 DP 了. 預處理環是學的別人的,感覺挺漲見識的.
#include<cstdlib> #include<cstdio> #include<algorithm> #include<cstring> using namespace std; const int maxn = 100; int n,m,sr,sb,sg,p,ans; int f[maxn][maxn][maxn]; int g[maxn],vis[maxn],cir[maxn]; inline int dfs(int x){return vis[x] ? 0 : (vis[x] = 1, dfs(g[x]) + 1); } inline void init(){for(int i = 0; i < maxn; ++i) vis[i] = cir[i] = 0;memset(f,0,sizeof(f)); ans = 0; } inline void add_mod(int &a,int b){a = (a + b) % p; } inline int exgcd(int a,int b,int p){return p % a ?(exgcd(b % a, a, (-p % a + a) % a) * b + p) / a : p / a; } int main() {freopen("card.in","r",stdin);freopen("card.out","w",stdout);scanf("%d %d %d %d %d\n",&sr,&sb,&sg,&m,&p);n = sr + sb + sg; m++;for(int i = 1; i <= m; ++i){for(int j = 1; j <= n; ++j)i == m ? g[j] = j : scanf("%d",&g[j]);init();for(int j = 1; j <= n; ++j) if(!vis[j]) cir[++cir[0]] = dfs(j);f[0][0][0] = 1;for(int j = 1; j <= cir[0]; ++j)for(int k = 0; k <= sr; ++k)for(int l = 0; l <= sb; ++l){int &F = f[j][k][l], t = n - k - l;if(k >= cir[j]) add_mod(F, f[j-1][k - cir[j]][l]);if(l >= cir[j]) add_mod(F, f[j-1][k][l - cir[j]]);if(t >= cir[j]) add_mod(F, f[j-1][k][l]);}ans = (ans + f[cir[0]][sr][sb]) % p;}ans = ans * exgcd(m, p, 1) % p;printf("%d\n",ans);return 0; }轉載于:https://www.cnblogs.com/Mr-ren/p/4297761.html
總結
- 上一篇: Android开发系列(二十八):使用S
- 下一篇: http协议之request