【BZOJ4833】最小公倍佩尔数(min-max容斥)

题面

BZOJ

题解

首先考虑怎么求\(f(n)\),考虑递推这个东西

\((1+\sqrt 2)(e(n-1)+f(n-1)\sqrt 2)=e(n)+f(n)\sqrt 2\)

拆开之后可以得到:\(e(n)=e(n-1)+2f(n-1),f(n)=f(n-1)+e(n-1)\)。

把每一层的\(e\)都给展开,得到:\(\displaystyle f(n)=1+f(n-1)+2\sum_{i=1}^{n-2}f(i)\)

然后差分搞搞,\(\displaystyle f(n)-f(n-1)=f(n-1)-f(n-2)+2*f(n-2)\)。

得到\(f(n)=2f(n-1)+f(n-2)\),特殊的\(f(0)=0,f(1)=1\)。

然后我们发现要求\(lcm\),那么就先考虑\(f(a)\)和\(f(b)\)的\(gcd\)是什么。

这个东西显然可以类似斐波那契数列那样子利用辗转相减得到\(gcd(f(a),f(b))=f(gcd(a,b))\)。

接下来就可以考虑怎么求答案了。

然后\(lcm\)的式子是对于每个质因子,考虑其\(max\)。

考虑\(min-max\)容斥,把\(max\)变成\(min\),那么就可以从\(lcm\)变成\(gcd\)。

然后把\(min-max\)容斥的式子给写出来:

\[max(S)=\sum_{T\subset S}(-1)^{|T|+1}min(T)
\]

套到\(lcm\)上就是:

\[lcm(S)=\prod_{T\subset S}gcd(T)^{(-1)^{|T|+1}}
\]

那么就有

\[g(n)=\prod_{T\subset S}f_{gcd(T)}^{(-1)^{|T|+1}}=\prod_{i=1}^n f_i^{\sum_{T\subset S}[gcd(T)=i](-1)^{|T|+1}}
\]

上面那个指数看着就可以莫比乌斯反演一下之类的,然后令上面那一堆东西是\(a[i]\),然后令\(b[i]=\sum_{i|d}a[d]\)这个系数稍微推一下,得到:

\[b[i]=\sum_{i|d}a[d]=\sum_{T\subset S}[i|gcd(T)](-1)^{|T|+1}
\]

这个值显然之和是否存在\(i\)倍数的数相关,存在就是\(1\),没有就是\(0\)。

而莫比乌斯反演可以得到

\[a[i]=\sum_{i|d}\mu(\frac{d}{i})b[d]
\]

再把这个东西带回去

\[\begin{aligned}
g[n]&=\prod_{i=1}^n f_i^{a[i]}\\
&=\prod_{i=1}^n f_i^{\sum_{i|d}\mu(\frac{d}{i})b[d]}\\
&=\prod_{i=1}^n\prod_{i|d}f_i^{\mu(\frac{d}{i})b[d]}
\end{aligned}\]

因为\(d\)的范围在\(n\)以内,所以必定存在\(d\)的倍数,所以\(b[d]=1\),那么只需要提前一个\(log\)预处理后面一半就行了。

#include<iostream>
#include<cstdio>
using namespace std;
#define MAX 1000100
inline int read()
{
int 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;
}
int n,MOD;
bool zs[MAX];
int pri[MAX],mu[MAX],tot;
int f[MAX],g[MAX],s[MAX],inv[MAX];
int fpow(int a,int b){int s=1;while(b){if(b&1)s=1ll*s*a%MOD;a=1ll*a*a%MOD;b>>=1;}return s;}
void Sieve(int 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{mu[i*pri[j]]=0;break;}
}
}
}
int main()
{
Sieve(MAX-1);
int T=read();
while(T--)
{
n=read();MOD=read();
f[1]=1;for(int i=2;i<=n;++i)f[i]=(2ll*f[i-1]+f[i-2])%MOD;
for(int i=1;i<=n;++i)s[i]=1,inv[i]=fpow(f[i],MOD-2);
for(int i=1;i<=n;++i)
for(int j=i;j<=n;j+=i)
if(mu[j/i]==1)s[j]=1ll*s[j]*f[i]%MOD;
else if(mu[j/i]==-1)s[j]=1ll*s[j]*inv[i]%MOD;
g[0]=1;for(int i=1;i<=n;++i)g[i]=1ll*g[i-1]*s[i]%MOD;
int ans=0;for(int i=1;i<=n;++i)ans=(ans+1ll*g[i]*i)%MOD;
printf("%d\n",ans);
}
}

【BZOJ4833】最小公倍佩尔数(min-max容斥)的更多相关文章

  1. BZOJ4833: [Lydsy1704月赛]最小公倍佩尔数(min-max容斥&莫比乌斯反演)(线性多项式多个数求LCM)

    4833: [Lydsy1704月赛]最小公倍佩尔数 Time Limit: 8 Sec  Memory Limit: 128 MBSubmit: 240  Solved: 118[Submit][S ...

  2. [Lydsy1704月赛] 最小公倍佩尔数

    4833: [Lydsy1704月赛]最小公倍佩尔数 Time Limit: 8 Sec  Memory Limit: 128 MBSubmit: 202  Solved: 99[Submit][St ...

  3. BZOJ4833: [Lydsy1704月赛]最小公倍佩尔数

    Problem 传送门 Sol 容易得到 \[f_n=e_{n-1}+f_{n-1},e_{n-1}=f_{n-1}+e_{n-1},f_1=e_1=1\] 那么 \[f_n=2\times \sum ...

  4. BZOJ 4833: [Lydsy1704月赛]最小公倍佩尔数(数论 + 最值反演)

    题面 令 \({(1+\sqrt 2)}^n=e(n)+f(n)*\sqrt2\) ,其中 \(e(n),f(n)\) 都是整数,显然有 \({(1-\sqrt 2)}^n=e(n)-f(n)*\sq ...

  5. 【bzoj 4833】[Lydsy1704月赛]最小公倍佩尔数

    Description 令 $(1+\sqrt 2)^n=e(n)+\sqrt 2\cdot f(n)$ ,其中 $e(n),f(n)$ 都是整数,显然有 $(1-\sqrt 2)^n=e(n)-\s ...

  6. [bzoj 4833]最小公倍佩尔数

    传送门 Description   Let \((1+\sqrt2)^n=e(n)+f(n)\cdot\sqrt2\) , both \(e(n)\) and \(f(n)\) are integer ...

  7. BZOJ2440(全然平方数)二分+莫比乌斯容斥

    题意:全然平方数是指含有平方数因子的数.求第ki个非全然平方数. 解法:比較明显的二分,getsum(int middle)求1-middle有多少个非全然平方数,然后二分.求1-middle的非全然 ...

  8. YYHS-分数(二分+容斥)

    题目描述 KJDH是个十分善于探索的孩子,有一天他把分子分母小于等于n的最简分数列在了纸上,他想找到这些分数里第k小的数,这对于KJDH来说当然是非常轻易,但是KJDH最近多了很多妹子,他还要去找妹子 ...

  9. Luogu5206 【WC2019】数树 【容斥,生成函数】

    题目链接 第一问白给. 第二问: 设 \(b=y^{-1}\),且以下的 \(Ans\) 是除去 \(y^n\) 的. 设 \(C(T)\) 是固定了 \(T\) 中的边,再连 \(n-|T|-1\) ...

随机推荐

  1. i春秋暑期训练营丨渗透测试工程师开课啦

    每个人的夏天 都有专属的解锁方式 或来一次难忘的旅行 或躺在家里吹着空调吃西瓜 又或者是和小伙伴参加暑期训练营 i春秋暑期渗透测试工程师 报名通道已全部开启 为了保证课程质量,采取小班教学,每班仅限3 ...

  2. undefined reference to `BN_new'

    出现如下错误 undefined reference to `BN_CTX_new' undefined reference to `BN_new' undefined reference to `B ...

  3. Unity导出Gradle工程给Android Studio使用

    1 Unity导出Gradle项目 Unity打包时Build System选择Gradle,勾选Export Project 2 Android Studio导入Unity导出的Gradle项目 打 ...

  4. 章节十四、5- web页面的截图

    一.以雅虎网站为例,当我们在登录时,输入错误的用户名然后点击“下一步”,用户名输入框会提示红色字体,这个时候我们就将页面进行截图. http://commons.apache.org/proper/c ...

  5. LeetCode——Rank Scores

    Write a SQL query to rank scores. If there is a tie between two scores, both should have the same ra ...

  6. python 之serial、pyusb 使用开发

    说明:本次是在windows 系统操作实现的. serial 使用场景,获取得力扫码枪的扫码数据,该扫码枪支持三种通讯接口设置,如下图 即插即用的是 USB-KBW功能,插上去即可获取扫码数据,第二种 ...

  7. Python:tarxjb简单、安全文件拷贝、传输

    tarxjb 简单.安全文件拷贝.传输 描述 通过python paramiko库实现简易ssh.sftp执行操作,从而实现文件的远程传输 Github 优点: 可靠传输,文件不易受损 安全传输,避免 ...

  8. 基于VLC库C#开发可播放摄像头及任意格式视频的播放器

    前言 本文主要讲述,在WPF中,借助Vlc.DotNet调用VLC类库,实现视频播功能,下面我们先来做开发前的准备工作. 准备工作 首先,我们创建一个项目WpfVLC,然后,进入Neget搜索Vlc. ...

  9. 如何使用jmockit进行单元测试

    1. Jmockit简介 JMockit 是用以帮助开发人员编写测试程序的一组工具和API,它完全基于 Java 5 SE 的 java.lang.instrument 包开发,内部使用 ASM 库来 ...

  10. github配置密钥

    我们在githob创建项目后,本地使用git 克隆代码 需要在githob配置密钥,才可以 步骤: 下载git,进行安装,安装好后.点击桌面,右键,选择>>git  bash 在弹出的黑框 ...