【BZOJ3202】项链(莫比乌斯反演,Burnside引理)
【BZOJ3202】项链(莫比乌斯反演,Burnside引理)
题面
题解
首先读完题目,很明显的感觉就是,分成了两个部分计算。
首先计算本质不同的珠子个数,再计算本质不同的项链个数。
前面一个部分和\(gcd\)相关,一种莫比乌斯反演的感觉。
后面一个部分出现了旋转操作,要求本质不同的方案数,不难想到Burnside引理。
首先先考虑怎么求本质不同的珠子个数。
我们直接考虑无序的三元组\((x,y,z)\),满足\(x,y,z\le a,gcd(x,y,z)=1\)
容斥考虑最终答案,就是\(\frac{1}{6}(S3-S2*3+S1*2)\),其中\(S3\)是无序的三元组个数,\(S2\)是满足互质的二元组个数,而\(S1\)则是三个数都相同并且\(gcd=1\)的三元组个数。
\(S1\)显然只有一个,即可以构成唯一合法三元组\((1,1,1)\)
考虑怎么计算\(S2\)和\(S3\)。
看起来\(s2\)很熟悉,也就是\(\sum_{i=1}^a\sum_{j=1}^a[gcd(i,j)=1]\)
然后直接莫比乌斯反演。
设\(F(x)=\sum_{i=1}^a\sum_{j=1}^a[gcd(i,j)=x]\)
\(G(x)=\sum_{x|d}F(d)=\sum_{i=1}^a\sum_{j=1}^a[x|gcd(i,j)]=[\frac{a}{x}]^2\)
然后又因为\(F(x)=\sum_{x|d}G(d)*\mu(\frac{d}{x})\)
要求的是\(F(1)\),即\(F(1)=\sum_{i=1}^aG(i)*\mu(i)=\sum_{i=1}^a[\frac{a}{i}]^2\mu(i)\)。
然后是\(S3\),即三元组个数,和上面一样的推法,可以得到贡献是\(\sum_{i=1}^a[\frac{a}{i}]^3\mu(i)\)
好了,那么计算本质不同的珠子的个数就算完了。
接下来考虑如何计算本质不同的项链的个数。
因为只有旋转计算本质不同,那么置换一共\(n\)个,是顺时针旋转\(i\)位。
那么答案就是所有置换的不动点个数除以\(n\)。那么不难往Polya定理方面靠。
只需要知道每个循环的大小即可。
对于顺时针旋转\(i\)位的循环长度显然就是\(\frac{n}{gcd(n,i)}\),循环个数\(gcd(n,i)\)也很显然。
因为是不动点的个数,因此每一个循环内的所有珠子的颜色必须相同。同时,题目限制相邻两个位置的珠子颜色不能相同。因此,这里可以等价的看做有\(gcd(n,i)\)个珠子,收尾相连,要求相邻的珠子不同色的方案数,假装这个的方案数是\(f(x)\),表示\(x\)个珠子收尾相连时候的方案数。那么现在要求的答案就是\(\sum_{i=1}^nf(gcd(i,n))\)。\(n\)太大了,考虑优化。
我们枚举\(gcd\),即\(\sum_{d=1}^n\sum_{i=1}^n[gcd(i,n)=d]f(d)\),那么显然\(d\)是\(n\)的因数。所以可以接着写成\(\sum_{d|n}\sum_{i=1}^n[gcd(i,n)=d]f(d)\),进一步推,也就是\(\sum_{d|n}\sum_{i=1,d|i}^n[gcd(\frac{i}{d},\frac{n}{d})=1]f(d)\),即需要知道所有和\(\frac{n}{d}\)互质的\(\frac{i}{d}\),考虑\(\frac{i}{d}\)的最大取值只有\(\frac{n}{d}\),所以这个值显然就是\(\varphi(\frac{n}{d})\)。因此,总的贡献就是\(\sum_{d|n}f(d)\varphi(\frac{n}{d})\)。
现在来考虑怎么计算\(f(x)\),假设一共有\(m\)中不同的珠子
首先我们可以在两个珠子之间插入一个不同的珠子,即\(f(x-1)*(m-2)\),亦或者强制让首位两个珠子颜色相同,然后插入一个不同颜色的珠子将他们隔开,这个操作等价于我现有\(f(x-2)\)个珠子构成了一个环,然后插入一个和首位相同的珠子,再插入一个颜色不同的珠子隔开,即\(f(x-2)*(m-1)\),那么转移就是\(f(x)=f(x-1)*(m-2)+f(x-2)*(m-1)\)。
显然可以矩阵乘法直接算。然而这种东西一般也可以推递推式的。
我们尽可能让他往一个等比的方向上靠,
那么我们最好能把式子化简为这种形式:\(f(x)+af(x-1)=bf(x-1)+abf(x-2)\)
那么列出方程组:
\]
所以有\(b=a+m-2\),所以有\(a(a+m-2)-(m-1)=0\),也就是\(a^2+(m-2)a-(m-1)=0\)
因式分解一下就是\((a+(m-1))(a-1)=0\)
解出来\(a=1\)或者\(a=1-m\),看着\(a=1\)好做些,
那么式子写成\(f(x)+f(x-1)=(m-1)*(f(x-1)+f(x-2))\)
所以得到\(f(x)+f(x-1)=(m-1)^{x-2}(f(1)+f(2))\),
显然有\(f(1)=0,f(2)=m*(m-1)\),为了方便,设\(S(x)=f(x)+f(x-1)\)。
递推式换成\(S(x)=(m-1)S(x-1)=(m-1)^{x-2}S(2)=(m-1)^{x-1}m\)
考虑一下\(f(x)\)怎么计算,\(f(x)=S(x)-f(x-1)=m(m-1)^{x-1}-f(x-1)\)
一般这样子都可以构建一个递推式。
f(x)&=m(m-1)^{x-1}-f(x-1)\\
f(x)+A(m-1)&=-(f(x-1)+A)\\
mA&=m(m-1)^{x-1}\\
A&=(m-1)^{x-1}
\end{aligned}
\]
所以,我们可以构建出这样一个等式:
\]
在考虑一下边界情况,可以得到\(f(x)-(m-1)^x=(-1)^x(m-1)\)
即\(f(x)=(-1)^x(m-1)+(m-1)^x\),这样子就可以快速幂计算\(f(x)\)了。
然而发现一些很不舒服的东西,因为\(n\)太大,导致\(n\)可能是模数的倍数,那就不能够直接除\(n\)了。
怎么办呢?那么首先我们直接模\(mod^2\),这样子最终的结果会是\(mod\)的倍数,所以可以先除掉\(mod\)然后再考虑逆元就行了。
离线一下就跑得飞快了,目前是洛谷rk1,bzoj rk3
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
#define ll long long
const ll mod=1000000007;
const ll MOD=mod*mod;
const ll inv6 = 833333345000000041ll;
#define MAX 10000001
inline ll read()
{
ll x=0;bool t=false;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=true,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return t?-x:x;
}
void add(ll &x,ll y){x+=y;if(x>=MOD)x-=MOD;}
ll Multi(ll x,ll y){return (x*y-(ll)(((long double)x*y+0.5)/(long double)MOD)*MOD+MOD)%MOD;}
ll fpow(ll a,ll b)
{
ll s=1;
while(b){if(b&1)s=Multi(s,a);a=Multi(a,a);b>>=1;}
return s;
}
ll qpow(ll a,ll b)
{
ll s=1;
while(b){if(b&1)s=s*a%mod;a=a*a%mod;b>>=1;}
return s;
}
ll sqr(ll x){return Multi(x,x);}
ll cub(ll x){return Multi(sqr(x),x);}
bool zs[MAX];
int pri[MAX/10],tot,mu[MAX];
ll n,m,ans;int a,cnt;
ll nn[50],aa[50];
ll p[MAX/10],q[MAX/10];
void pre(ll N)
{
mu[1]=1;
for(int i=2;i<=N;++i)
{
if(!zs[i])pri[++tot]=i,mu[i]=-1;;
for(int j=1;j<=tot&&i*pri[j]<=N;++j)
{
zs[i*pri[j]]=true;
if(i%pri[j])mu[i*pri[j]]=-mu[i];
else break;
}
}
for(int i=1;i<=N;++i)mu[i]+=mu[i-1];
}
ll Calc(ll n)
{
ll S2=0,S3=0;
for(ll i=1,j;i<=n;i=j+1)
{
j=n/(n/i);
add(S3,Multi(cub(n/i),(mu[j]-mu[i-1]+MOD)%MOD));
add(S2,Multi(sqr(n/i),(mu[j]-mu[i-1]+MOD)%MOD));
}
add(S3,Multi(S2,3)),add(S3,2);
return Multi(S3,inv6);
}
ll F(ll n)
{
ll ret=fpow(m-1,n);
if(n&1)add(ret,(MOD-(m-1))%MOD);
else add(ret,m-1);
return ret;
}
void dfs(int i,ll d,ll phi)
{
if(i==cnt+1){add(ans,Multi(phi,F(n/d)));return;}
dfs(i+1,d,phi);
d*=p[i];phi*=p[i]-1;dfs(i+1,d,phi);
for(int x=2;x<=q[i];++x)
d*=p[i],phi*=p[i],dfs(i+1,d,phi);
}
int main()
{
int T=read();ll mx=0;
for(int i=1;i<=T;++i)nn[i]=read(),aa[i]=read();
for(int i=1;i<=T;++i)mx=max(mx,max((ll)sqrt(nn[i]),aa[i]));
pre(mx);
for(int TT=1;TT<=T;++TT)
{
n=nn[TT];a=aa[TT];
m=Calc(a);cnt=0;ll x=n;
for(int i=1;i<=tot&&1ll*pri[i]*pri[i]<=x;++i)
if(x%pri[i]==0)
{
p[++cnt]=pri[i];q[cnt]=0;
while(x%pri[i]==0)++q[cnt],x/=pri[i];
}
if(x>1)p[++cnt]=x,q[cnt]=1;
ans=0;dfs(1,1,1);
if(n%mod==0)ans=(ll)(ans/mod)*qpow(n/mod,mod-2)%mod;
else ans=(ans%mod)*qpow(n%mod,mod-2)%mod;
printf("%lld\n",ans);
}
return 0;
}
【BZOJ3202】项链(莫比乌斯反演,Burnside引理)的更多相关文章
- P3307-[SDOI2013]项链【Burnside引理,莫比乌斯反演,特征方程】
正题 题目链接:https://www.luogu.com.cn/problem/P3307 题目大意 \(n\)个珠子的一个环形项链,每个珠子有三个\(1\sim k\)的整数. 两个珠子不同当且仅 ...
- 洛谷P3307 [SDOI2013]项链 [polya定理,莫比乌斯反演]
传送门 思路 很明显的一个思路:先搞出有多少种珠子,再求有多少种项链. 珠子 考虑这个式子: \[ S3=\sum_{i=1}^a \sum_{j=1}^a\sum_{k=1}^a [\gcd(i,j ...
- 洛谷 P3307 - [SDOI2013]项链(Burnside 引理+数论)
题面传送门 看到题目我们显然可以将题目拆分成两部分:首先求出有多少个符合要求的珠子 \(c\),这样我们就可以将每种珠子看成一种颜色,题目也就等价于有多少种用 \(c\) 种颜色染长度为 \(n\) ...
- 【POJ 2154】 Color (置换、burnside引理)
Color Description Beads of N colors are connected together into a circular necklace of N beads (N< ...
- [模板] 容斥原理: 二项式反演 / Stirling 反演 / min-max 容斥 / 子集反演 / 莫比乌斯反演
//待更qwq 反演原理 二项式反演 若 \[g_i=\sum_{j=1}^i {\binom ij} f_j\] , 则有 \[ f_i=\sum_{j=1}^i (-1)^{i-j} {i \ch ...
- Burnside引理与Polya定理 学习笔记
原文链接www.cnblogs.com/zhouzhendong/p/Burnside-Polya.html 问题模型 有一个长度为 $n$ 的序列,序列中的每一个元素有 $m$ 种取值. 如果两个序 ...
- 【POJ2888】Magic Bracelet Burnside引理+欧拉函数+矩阵乘法
[POJ2888]Magic Bracelet 题意:一个长度为n的项链,有m种颜色的珠子,有k个限制(a,b)表示颜色为a的珠子和颜色为b的珠子不能相邻,求用m种珠子能串成的项链有多少种.如果一个项 ...
- 【uva 10294】 Arif in Dhaka (First Love Part 2) (置换,burnside引理|polya定理)
题目来源:UVa 10294 Arif in Dhaka (First Love Part 2) 题意:n颗珠子t种颜色 求有多少种项链和手镯 项链不可以翻转 手镯可以翻转 [分析] 要开始学置换了. ...
- Polya定理与Burnside引理
也许更好的阅读体验 \(Burnside引理\) 公式 \(\begin{aligned}L=\frac{1}{|G|}\sum_{i=1}^{|G|}D_{G_i}\end{aligned}\) 一 ...
随机推荐
- 学习angularjs的ng-hide和ng-disabled
一,页面上有一个checkbox和一个文本框.切换checkbox能对文本框输入文本与否: <input type="checkbox" ng-model="ckS ...
- ActiveMQ在C#中的应用
本文是在.NET Framework框架下的应用,截止到目前ActiveMQ还不支持.NET Core,而RabbitMQ已经支持.NET Core,希望ActiveMQ能尽快支持. ActiveMQ ...
- 运行supervisord -c /etc/supervisor/supervisord.conf 出错,解决办法
坑都让我踩了...... 1 supervisord -c /etc/supervisor/supervisord.conf 什么意思? 答:手动启动:supervisord 具体详见 ...
- SKINNY加密算法详解(无代码,仅加密)
原作者论文请参考<The SKINNY Family of Block Ciphers and Its Low-Latency Variant MANTIS> 地址为:https://li ...
- Supervisor (进程管理利器) 使用说明 - 运维笔记
一.Supervisor简单介绍supervisor是一个 Client/Server模式的系统,允许用户在类unix操作系统上监视和控制多个进程,或者可以说是多个程序.supervisor与laun ...
- Mysql读写分离方案-MySQL Proxy环境部署记录
Mysql的读写分离可以使用MySQL Proxy和Amoeba实现,其实也可以使用MySQL-MMM实现读写分离的自动切换.MySQL Proxy有一项强大功能是实现"读写分离" ...
- 第九次Scrum meeting
第九次Scrum meeting 任务及完成度: 成员 12.31 1.1 陈谋 任务1040:完成stackoverflow的数据处理后的json处理(99%) 任务1114-1:完成对网页数据的 ...
- Linux内核及分析 第四周 扒开系统调用的三层皮(上)
实验过程 选择20号系统调用getpid(取得进程识别码) 在网上查询getpid函数的C语言代码以及其嵌入式汇编语句 C语言代码: #include <stdio.h> #include ...
- JAVA链表中迭代器的实现
注:本文代码出自<java数据结构和算法>一书. PS:本文中类的名字定义存在问题,Link9应改为Link.LinkList9应该为LinkList.由于在同包下存在该名称,所以在后面接 ...
- 总结 推广app
扫一扫二维码即可安装使用我们的app,方便快捷. 电脑端下载地址:http://pan.baidu.com/s/1bocWPPX http://a.app.qq.com/o/simple.jsp?pk ...