传送门

表示这种\(Burnside\)定理之类的东西一用就忘qwq

题目要求不同染色方案数,因为变换方式只有旋转,所以只有\(n\)个置换,然后可能会出现某种方案有循环节,这个循环节长度显然要是\(\gcd(n,m)\)的因数,我们枚举循环节个数,直接套个polya然后可以得到答案为\(\frac{\sum_{d|\gcd(n,m)}f(\frac{n}{d},\frac{m}{d})\varphi(d)}{n}\),其中\(f(n,m)\)为有\(n\)个珠子,要把\(m\)个染色,并且满足最长连续的黑色个数以及首尾黑色个数\(\le k\)的方案

这个东西,可以考虑把白色的珠子提出来,然后枚举首尾放了多少个黑色\(i\),这样有\(i+1\)种方案,剩下的黑色就要放在\(n-m-1\)个空位中,并且满足一个空位不超过\(k\)个黑色.这个可以容斥,枚举一些位置已经超过\(k\)个,然后剩下随便放.这个式子长这样:$$\sum_{i=0}{k}(i+1)\sum_{j=0}{n-m-1}(-1)^j\binom{n-m-1}{j}g(m-i-j*(k+1),n-m-1)$$

其中\(g(n,m)\)表示\(n\)个一样的物品放进\(m\)个不一样的盒子的方案数,也就是\(\binom{n+m-1}{m-1}\)

代码巨丑

#include<bits/stdc++.h>
#define LL long long
#define db long double
#define il inline
#define re register using namespace std;
const int N=1e5+10,mod=998244353;
il LL rd()
{
LL x=0,w=1;char ch=0;
while(ch<'0'||ch>'9') {if(ch=='-') w=-1;ch=getchar();}
while(ch>='0'&&ch<='9') {x=(x<<3)+(x<<1)+(ch^48);ch=getchar();}
return x*w;
}
int gcd(int a,int b){return b?gcd(b,a%b):a;}
int fac[N],iac[N];
int C(int n,int m){return m<0||n<m?0:1ll*fac[n]*iac[m]%mod*iac[n-m]%mod;}
int fpow(int a,int b){int an=1;while(b){if(b&1) an=1ll*an*a%mod;a=1ll*a*a%mod,b>>=1;} return an;}
int n,m,nn,kk,prm[20],nm[20],tt,ans;
int wk(int n,int m)
{
if(m<=kk) return C(n,m);
int an=0,l=n-m-1;
for(int i=0;i<=kk;++i)
{
int nw=0;
for(int j=0;j<=l&&m-i-j*(kk+1)>=0;++j)
{
int x=1ll*C(l,j)*C(m-i-j*(kk+1)+l-1,l-1)%mod;
nw=(nw+((j&1)?mod-x:x))%mod;
}
an=(an+1ll*nw*(i+1)%mod)%mod;
}
return an;
}
void dfs(int o,int d,int phi)
{
if(o>tt)
{
ans=(ans+1ll*wk(n/d,m/d)*phi%mod)%mod;
return;
}
dfs(o+1,d,phi);
for(int i=1;i<=nm[o];++i)
{
d*=prm[o];
dfs(o+1,d,phi*(prm[o]-1)*fpow(prm[o],i-1));
}
} int main()
{
fac[0]=1;
for(int i=1;i<=N-10;++i) fac[i]=1ll*fac[i-1]*i%mod;
iac[N-10]=fpow(fac[N-10],mod-2);
for(int i=N-10;i;--i) iac[i-1]=1ll*iac[i]*i%mod;
n=rd(),m=rd(),kk=rd();
if(n==m) return printf("%d\n",(bool)!kk),0;
nn=gcd(n,m);
int x=nn,sqt=sqrt(nn);
tt=0;
for(int i=2;i<=sqt&&x>1;++i)
if(x%i==0)
{
prm[++tt]=i,nm[tt]=0;
while(x%i==0) ++nm[tt],x/=i;
}
if(x>1) prm[++tt]=x,nm[tt]=1;
ans=0,dfs(1,1,1);
printf("%lld\n",1ll*ans*fpow(n,mod-2)%mod);
return 0;
}

luogu P4916 魔力环的更多相关文章

  1. [jzoj 6084] [GDOI2019模拟2019.3.25] 礼物 [luogu 4916] 魔力环 解题报告(莫比乌斯反演+生成函数)

    题目链接: https://jzoj.net/senior/#main/show/6084 https://www.luogu.org/problemnew/show/P4916 题目: 题解: 注: ...

  2. GOOD BYE OI

    大米饼正式退役了,OI给我带来很多东西 我会的数学知识基本都在下面了 博客园的评论区问题如果我看到了应该是会尽力回答的... 这也是我作为一个OIer最后一次讲课的讲稿 20190731 多项式乘法 ...

  3. Luogu 魔法学院杯-第二弹(萌新的第一法blog)

    虽然有点久远  还是放一下吧. 传送门:https://www.luogu.org/contest/show?tid=754 第一题  沉迷游戏,伤感情 #include <queue> ...

  4. luogu p1268 树的重量——构造,真正考验编程能力

    题目链接:http://www.luogu.org/problem/show?pid=1268#sub -------- 这道题费了我不少心思= =其实思路和标称毫无差别,但是由于不习惯ACM风格的题 ...

  5. [luogu P2170] 选学霸(并查集+dp)

    题目传送门:https://www.luogu.org/problem/show?pid=2170 题目描述 老师想从N名学生中选M人当学霸,但有K对人实力相当,如果实力相当的人中,一部分被选上,另一 ...

  6. [luogu P2647] 最大收益(贪心+dp)

    题目传送门:https://www.luogu.org/problem/show?pid=2647 题目描述 现在你面前有n个物品,编号分别为1,2,3,--,n.你可以在这当中任意选择任意多个物品. ...

  7. Luogu 考前模拟Round. 1

    A.情书 题目:http://www.luogu.org/problem/show?pid=2264 赛中:sb题,直接暴力匹配就行了,注意一下读入和最后一句话的分句 赛后:卧槽 怎么只有40 B.小 ...

  8. luogu P2580 于是他错误的点名开始了

    luogu  P2580 于是他错误的点名开始了 https://www.luogu.org/problem/show?pid=2580 题目背景 XS中学化学竞赛组教练是一个酷爱炉石的人. 他会一边 ...

  9. CJOJ 1331 【HNOI2011】数学作业 / Luogu 3216 【HNOI2011】数学作业 / HYSBZ 2326 数学作业(递推,矩阵)

    CJOJ 1331 [HNOI2011]数学作业 / Luogu 3216 [HNOI2011]数学作业 / HYSBZ 2326 数学作业(递推,矩阵) Description 小 C 数学成绩优异 ...

随机推荐

  1. 斯坦福大学公开课机器学习: neural networks learning - autonomous driving example(通过神经网络实现自动驾驶实例)

    使用神经网络来实现自动驾驶,也就是说使汽车通过学习来自己驾驶. 下图是通过神经网络学习实现自动驾驶的图例讲解: 左下角是汽车所看到的前方的路况图像.左上图,可以看到一条水平的菜单栏(数字4所指示方向) ...

  2. error:crosses initialization of ...

    switch(c) { case 0x01: int temp = a + b; .... break; case 0x02: break; default:break; } 此时会报如题所示错误 原 ...

  3. typescript和coffeescript简介

    typescript 1.什么是typescript?和JavaScript有什么关系? 1.它是一门独立的语言,实现了自己的编译器   2.它的目标代码是JavaScript(很多语言的目标代码是机 ...

  4. JS模块化开发(三)——seaJs+grunt

    1.seaJs直接构建存在的问题 由于模块之间的依赖require引用的是模块名,当多个js模块被合并成一个时,会由于找不到模块名而报错 2.seaJs+grunt开发 用到的插件:grunt-cmd ...

  5. Error:Failed to resolve: :Base:

    这个问题是变通了一下,原来是: //implementation(name: "Base", ext: "aar") 修改成: implementation f ...

  6. Vue(基础三)_监听器与计算属性

    一.前言 本文主要涉及:     1.watch()监听单个属性 2.computed可以监听多个属性 3.点击li标签切换音乐案例 二.主要内容 1.watch()监听器(监听单一数据) (1)监听 ...

  7. linux driver ------ 字符设备驱动 之 “ 创建设备节点流程 ”

    在字符设备驱动开发的入门教程中,最常见的就是用device_create()函数来创建设备节点了,但是在之后阅读内核源码的过程中却很少见device_create()的踪影了,取而代之的是device ...

  8. ECharts基础

    echarts: js引用:<script type="text/javascript" src="js/echarts.js"></scri ...

  9. go etcd

    etcd介绍 GitHub:https://github.com/coreos/etcd 官网:https://coreos.com/etcd/ 下载:https://github.com/coreo ...

  10. ipv4转化为ipv6

    十進制轉換成十六進位 IPV6為十六進位,所以十進制轉換成十六進位192=c0 168=a8192.168.1.1 轉成 16 進制為 c0.a8.01.01可以使用 Windows 工程版或是程式設 ...