题面

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

令 \(g(n)\) 表示 \(f(1),f(2)…f(n)\) 的最小公倍数,给定两个正整数 \(n\) 和 \(p\) ,其中 \(p\) 是质数,并且保证 \(f(1),f(2)…f(n)\) 在模 \(p\) 意义下均不为 \(0\) .

请计算\(\displaystyle \sum_{i=1}^{n} i*g(i)\) .其在模 \(p\) 的值。

\(T \le 210, \sum n \le 3 \times 10^6\)

题解

其实它就是给你了 \(f(n)\) 的特征方程,其实可以解出来,具体来说 递推式到通相式 直接待定系数,解回去就把系数代回去就行了。

但是没有必要,\(f(n)\) 是个能打表的东西。

\[f(n) = \begin{cases}
0, &n = 0 \\
1, &n = 1 \\
2f(n - 1) + f(n-2) &\mathrm{otherwise}
\end{cases}
\]

然后它具有和 \(fib(i)\) 一样的性质, \(\gcd (f(i), f(j)) = f(\gcd(i, j))\) 。

至于证明,可以类似 \(fib(i)\) 的性质证明。

然后考虑最值反演:

不会的话可以看我 这篇博客

然后我们可以把这个对于 \(\mathrm{lcm}\) 放到指数上来就是所有指数的 \(\max\) ,那么我们可以用最值反演变成 \(\min\) ,也就是 \(\gcd\) 。

形式如下

\[\mathrm{lcm}(S) = \prod_{T \subseteq S, T \not = \varnothing} \gcd(T)^{(-1)^{|T|-1}}
\]

然后我们利用这个性质进行推导。

\[\begin{align}
g(n) &= \prod_{T \subseteq S, T \not = \varnothing} \gcd_{x \in T}\{f(x)\}^{(-1)^{|T|+1}} \\
&= \prod_{T \subseteq S, T \not = \varnothing} f(\gcd(T))^{(-1)^{|T|+1}}\\
\end{align}
\]

然后考虑构造一个函数 \(h\) 满足

\[f(n) = \prod_{d|n} h(d)
\]

由于在模素数意义下,所以必定存在 \(h\) 满足要求。

然后我们考虑每个 \(h(d)\) 的贡献:

\[\begin{align}
\Rightarrow
g(n) &= \prod_{d = 1}^{n} h(d)^{\sum_{T\subseteq S, T \not = \varnothing}[d | \gcd(T)](-1)^{|T|+1}} \\
&= \prod_{d = 1}^{n} h(d)
\end{align}
\]

至于上面那个系数为什么 \(=1\) 呢?可以考虑 \(T\) 集合中所有为 \(d\) 倍数的数构成的所有集合,刚好是 \(2^k - 1\) 的,除去空集,那么剩下的一一抵消,剩下一个 \(1\) 的贡献。

考虑求出 \(h\) 是可以枚举倍数来优化到 \(O(n \ln n)\) 的。

具体来说假设我们 \(<n\) 的所有 \(h\) 都求好了,那么当前的 \(\displaystyle h(n) = f(n) \prod_{d|n, d \not = n} \frac{1}{h(d)}\) 。

其实可以用莫比乌斯反演推回去的,但没必要筛 \(\mu\) 。

那么能不能用线性筛优化到 \(O(n)\) 呢?似乎是不行的,因为 \(h\) 不是积性的,我们知道 \(h(i)\) 与 \(h(p)\) 无法推出 \(h(i \times p)\) 。

然后最后复杂度就是 \(O(\sum n \ln n)\) 的。

总结

对于有些看起来不太好求的 \(\mathrm{lcm}\) 题,考虑最值反演。然后考虑构造积性函数,最后求出积性函数求回来就行了。

代码

/**************************************************************
Problem: 4833
User: zjp_shadow
Language: C++
Result: Accepted
Time:5908 ms
Memory:13024 kb
****************************************************************/ #include <bits/stdc++.h> #define For(i, l, r) for(register int i = (l), i##end = (int)(r); i <= i##end; ++i)
#define Fordown(i, r, l) for(register int i = (r), i##end = (int)(l); i >= i##end; --i)
#define Set(a, v) memset(a, v, sizeof(a))
#define Cpy(a, b) memcpy(a, b, sizeof(a))
#define debug(x) cout << #x << ": " << (x) << endl
#define DEBUG(...) fprintf(stderr, __VA_ARGS__) using namespace std; typedef long long ll; template<typename T> inline bool chkmin(T &a, T b) {return b < a ? a = b, 1 : 0;}
template<typename T> inline bool chkmax(T &a, T b) {return b > a ? a = b, 1 : 0;} inline int read() {
int x(0), sgn(1); char ch(getchar());
for (; !isdigit(ch); ch = getchar()) if (ch == '-') sgn = -1;
for (; isdigit(ch); ch = getchar()) x = (x * 10) + (ch ^ 48);
return x * sgn;
} void File() {
#ifdef zjp_shadow
freopen ("4883.in", "r", stdin);
freopen ("4883.out", "w", stdout);
#endif
} const int N = 1e6 + 1e3; int Mod, h[N], g[N], ans[N]; ll fpm(ll x, int power) {
ll res = 1;
for (; power; power >>= 1, (x *= x) %= Mod)
if (power & 1) (res *= x) %= Mod;
return res;
} void Solve(int maxn) {
h[1] = 1;
For (i, 2, maxn)
h[i] = (h[i - 1] * 2ll + h[i - 2]) % Mod;
For (i, 2, maxn) {
int inv = fpm(h[i], Mod - 2);
for (int j = i * 2; j <= maxn; j += i)
h[j] = 1ll * h[j] * inv % Mod;
}
g[0] = 1; int ans = 0;
For (i, 1, maxn) {
g[i] = 1ll * g[i - 1] * h[i] % Mod;
ans = (ans + 1ll * g[i] * i) % Mod;
}
printf ("%d\n", ans);
} int main () { File(); int cases = read();
while (cases --) {
int n = read(); Mod = read(); Solve(n);
} return 0; }

BZOJ 4833: [Lydsy1704月赛]最小公倍佩尔数(数论 + 最值反演)的更多相关文章

  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. 【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 ...

  4. 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 ...

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

    [BZOJ4833]最小公倍佩尔数(min-max容斥) 题面 BZOJ 题解 首先考虑怎么求\(f(n)\),考虑递推这个东西 \((1+\sqrt 2)(e(n-1)+f(n-1)\sqrt 2) ...

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

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

  7. bzoj 4831 [Lydsy1704月赛]序列操作 dp

    [Lydsy1704月赛]序列操作 Time Limit: 1 Sec  Memory Limit: 128 MBSubmit: 203  Solved: 69[Submit][Status][Dis ...

  8. bzoj 4836 [Lydsy1704月赛]二元运算 分治FFT+生成函数

    [Lydsy1704月赛]二元运算 Time Limit: 8 Sec  Memory Limit: 128 MBSubmit: 577  Solved: 201[Submit][Status][Di ...

  9. BZOJ.4832.[Lydsy1704月赛]抵制克苏恩(期望DP)

    题目链接 \(f[s][i][j][k]\)表示还剩\(s\)次攻击,分别有\(i,j,k\)个血量为\(1,2,3\)的奴隶主时,期望受到伤害. 因为期望是倒推,所以这么表示从后往前求,注意\(a, ...

随机推荐

  1. 16-使用Selenium模拟浏览器抓取淘宝商品美食信息

    淘宝由于含有很多请求参数和加密参数,如果直接分析ajax会非常繁琐,selenium自动化测试工具可以驱动浏览器自动完成一些操作,如模拟点击.输入.下拉等,这样我们只需要关心操作而不需要关心后台发生了 ...

  2. Survey项目总结

    1.Ioc深入理解 Inverse of control org.springframework.scheduling.quartz.SchedulerFactoryBean org.mybatis. ...

  3. net core 小坑杂记之配置文件读取(不定期更新)

    其实很早就想写了,原想等积累差不多了再写的,但是发现遇到一个当时记下效果会比较好,所以就不定期更新这个系列了,后面获取会整个整理一下. 此篇记载net core入门时踩的一些坑,网上教程太少了,也不规 ...

  4. Python之字符串格式化

    1)     占位符%s: %s是通用的占位符,所有类型不管是string还是int还是float全都代表. 如果使用%d,则只能代表整数:如果是%f,则只能代表小数: 2)     直接用加号+连接 ...

  5. 输入input

    用input接收到的类型全部都是字符串!!! 要查看变量类型,可以使用type()模块: 字符串不能和数字进行比较,因此如果输入是以input方式输入的,需要先转换成数字格式:

  6. MYSQL业务数据简单脱敏方案

    removesensitiveinformationplan.sh #!/bin/sh #!在模拟库上运行如下脚本 #!生产库crm-db #!模拟库crm-mock #!.备份生产库 mysqldu ...

  7. C# 和 c++的语法不同点

    GC  Garbage Collection 垃圾回收器 自动释放资源 关键字: new 1.创建对象 2.隐藏从父类继承的同名函数 using 1.引用命名空间 2. using(FileStrea ...

  8. 如何入门vue之一

    入门vue 首先从vue的指令开始学起. vue的指令: v-if  根据是否得到的布尔值进行是否显示. v-show:根据是否得到的布尔值是否显示.不同的地方在于隐藏是style隐藏. v-on 监 ...

  9. Linux df 与du用法

    df 查看一级目录的大小,但是不能查看文件的大小.du 可以查看目录或者文件大小. 1 df的基本用法 df命令可以显示目前所有文件系统的总空间及当前可用空间,用法如下: -a 全部文件系统列表-h ...

  10. 一条SQL语句执行得很慢的原因有哪些?(转)

    一条 SQL 语句执行的很慢,那是每次执行都很慢呢?还是大多数情况下是正常的,偶尔出现很慢呢?所以我觉得,我们还得分以下两种情况来讨论. 1.大多数情况是正常的,只是偶尔会出现很慢的情况. 2.在数据 ...