bzoj 1004 Cards & poj 2409 Let it Bead —— 置换群

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1004

关于置换群：https://www.cnblogs.com/nietzsche-oier/p/6883880.html

https://files-cdn.cnblogs.com/files/HocRiser/Burnside.pdf

原来 burnside 引理中的“不动点”是指一种不变化的方案啊；

这道题就用 burnside 引理，但给出的 m 个置换还不是置换群，需要再加一个单位元，即 \( a[i] = i \) 的置换；

用三维DP求出每种置换的“不动点”个数，枚举可以通过记录总数而减少一维。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const xn=,xm=;
int n,m,sr,sb,sg,a[xn],f[xm][xm][xm],mod,tot;
bool vis[xn];
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return f?ret:-ret;
}
ll pw(ll a,int b)
{
  ll ret=;
  for(;b;b>>=,a=(a*a)%mod)if(b&)ret=(ret*a)%mod;
  return ret;
}
int upt(int x){while(x>=mod)x-=mod; while(x<)x+=mod; return x;}
void dfs(int x)
{
  tot++; vis[x]=;
  if(!vis[a[x]])dfs(a[x]);
}
int main()
{
  sr=rd(); sb=rd(); sg=rd(); m=rd(); mod=rd();
  n=sr+sb+sg; m++;
  for(int i=;i<=m;i++)
    {
      memset(vis,,sizeof vis);
      memset(f,,sizeof f); f[][][]=;//
      if(i==m)for(int j=;j<=n;j++)a[j]=j;
      else for(int j=;j<=n;j++)a[j]=rd(),vis[j]=;
      tot=; int res=;
      for(int j=;j<=n;j++)
    {
      if(vis[j])continue;
      tot=; dfs(j); res+=tot;
      for(int j=sr;j>=;j--)
        for(int k=sb;k>=;k--)
          //for(int l=sg;l>=0;l--)
          {
        int l=res-j-k; if(l<)continue;
        if(j>=tot)f[j][k][l]=upt(f[j][k][l]+f[j-tot][k][l]);
        if(k>=tot)f[j][k][l]=upt(f[j][k][l]+f[j][k-tot][l]);
        if(l>=tot)f[j][k][l]=upt(f[j][k][l]+f[j][k][l-tot]);
          }
    }
    }
  printf("%lld\n",(ll)f[sr][sb][sg]*pw(m,mod-)%mod);
  return ;
}

题目：http://poj.org/problem?id=2409

Polya 定理裸题；

gcd 可以从模的剩余系的角度来看；

注意加上翻转后有 2n 个置换。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int const xn=;
int n,m;
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int pw(int a,int b)
{
  int ret=;
  for(;b;b>>=,a=a*a)if(b&)ret=ret*a;
  return ret;
}
int main()
{
  while()
    {
      scanf("%d%d",&m,&n);
      if(!m&&!n)return ;
      int ans=;
      for(int i=;i<=n;i++)ans+=pw(m,gcd(n,i));
      if(n%)ans+=n*pw(m,n/+);
      else ans+=n/*pw(m,n/)+n/*pw(m,n/+);
      printf("%d\n",ans//n);//2n
    }
  return ;
}