【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. gsoap生成webservice调用客户端接口

    1.下载gsoap2.8 2.运行 wsdl2h.exe -o XXX.h XXX.wsdl wsdl文件可以是本地文件,也可以是服务器的wsdl,比如http://192.168.0.122:333 ...

  2. GCN实现3

    参考 : 首先看两篇论文,大概了解一下原理性的东西: GRAPH CONVOLUTIONAL NETWORKS THOMAS KIPF, 30 SEPTEMBER 2016 http://tkipf. ...

  3. GCN 简单numpy实现

    `#参考:https://blog.csdn.net/weixin_42052081/article/details/89108966 import numpy as np import networ ...

  4. Java数据库连接组件C3P0和DBCP

    C3P0连接池组件 开源数据库连接池组件 jar包:c3p0-0.9.2.jar和mchange-commons-java-0.2.2.3.jar 支持JDBC3规范和JDBC2的标准扩展 使用项目H ...

  5. Centos 7 下yum搭建lnmp环境(yum安装方式)

    我们都知道linux下安装软件主要有三种方式: 1.源码编译安装,即下载软件源代码,利用gcc g++ make 等编译工具进行编译安装: 此方式的优点:可以指定软件版本,可选择性好:编译时可以手动指 ...

  6. 4.Spark环境搭建和使用方法

    一.安装Spark spark和Hadoop可以部署在一起,相互协作,由Hadoop的HDFS.HBase等组件复制数据的存储和管理,由Spark负责数据的计算. Linux:CentOS Linux ...

  7. mybatis多对多关联查询

    多对多关系 一个学生可以选多门课程,而一门课程可以由多个学生选择,这就是一个典型的多对多关联关系.所谓多对多关系,其实是由两个互反的一对多关系组成.即多对多关系都会通过一个中间表来建立,例如选课表.学 ...

  8. 创建war类型项目(六)

    1. 创建maven project时选择packaging为war: 2. 在webapp文件夹下新建META-INF和WEB-INF/web.xml web.xml模板 <?xml vers ...

  9. Django简介(MVC、MTV)

    Django简介 MVC Model(模型)- 应用程序中处理数据逻辑部分且与数据库交互,用于存取数据的部分 View(视图)- 用于处理后的数据界面展示,且视图通常是由模型数据创建的,是用户看到并与 ...

  10. JDOJ 1044 Span

    JDOJ 1044 Span https://neooj.com/oldoj/problem.php?id=1044 Description 某国有N个村子,M条道路,为了实现“村村通工程”现在要”油 ...