题目 对于正整数n，定义f(n)为n所含质因子的最大幂指数。例如f(1960)=f(2^3 5^1 7^2)=3, f(10007)=1, f(1)=0。 给定正整数a,b，求sigma(sigma(f(gcd(i,j)))) (i=1..a, j=1..b)。 输入格式 第一行一个数T，表示询问数 ...
分类：
其他好文 时间：
2018-04-30 21:13:27
阅读次数：
154
题目大意 给出一段序列，求其中最大公约数为1的四元组的个数。 思路 我们要用到反演、正难则反的思想。对于每一个大于1的数字$x$，求出最大公约数为$x$的四元组的个数$g(x)$，然后用排列中所有四元组的组合个数减去$\sum g(x)$即可。 直接求$g(x)$没有什么思路，但是求公约数中存在$x ...
分类：
其他好文 时间：
2018-04-28 01:19:49
阅读次数：
165
题目链接：https://vjudge.net/problem/SPOJ-PGCD 题目大意： 给定 \(N\) 和 \(M\)，求满足 \((1 \le x \le N), (1 \le y \le M)\)，且 \(gcd(x,y)\) 为素数的 \((x,y)\) 的对数。 知识点： 莫比乌斯 ...
分类：
其他好文 时间：
2018-04-27 02:12:51
阅读次数：
187
BZOJ_4176_Lucas的数论_杜教筛+莫比乌斯反演 Description 去年的Lucas非常喜欢数论题，但是一年以后的Lucas却不那么喜欢了。 在整理以前的试题时，发现了这样一道题目“求Sigma(f(i)),其中1<=i<=N”，其中 表示i的约数个数。他现在长大了，题目也变难了。 ...
分类：
其他好文 时间：
2018-04-22 13:07:26
阅读次数：
162
BZOJ_3994_[SDOI2015]约数个数和_莫比乌斯反演 Description 设d(x)为x的约数个数，给定N、M，求 设d(x)为x的约数个数，给定N、M，求 Input 输入文件包含多组测试数据。 第一行，一个整数T，表示测试数据的组数。 接下来的T行，每行两个整数N、M。 输入文件 ...
分类：
其他好文 时间：
2018-04-22 12:58:35
阅读次数：
148
题面 "传送门" Sol 莫比乌斯反演一波就是求 $$\sum_{k=1}^{min(a,b)}\lfloor\frac{a}{k}\rfloor\lfloor\frac{b}{k}\rfloor\sum_{d|k}f(d)\mu(\frac{k}{d})$$ 前面的分块就好了 考虑求 $$\sum ...
分类：
其他好文 时间：
2018-04-12 22:39:37
阅读次数：
226
【BZOJ3309】DZY Loves Math（莫比乌斯反演） 题面 求 $$\sum_{i=1}^a\sum_{j=1}^bf(gcd(a,b))$$ 其中，$f(x)$表示$x$分解质因数之后，最高的幂次 题解 完全不会莫比乌斯反演了。 先来推式子 $$\sum_{d=1}^a\sum_{i= ...
分类：
其他好文 时间：
2018-04-12 15:26:25
阅读次数：
145
"题目链接" cpp //928kb 56ms //注意d！ include include typedef long long LL; const int N=1 3]; LL n; bool Not_p[N+3]; void Make_Table(int N) { for(int i=2; i1 ...
分类：
其他好文 时间：
2018-04-08 18:03:36
阅读次数：
158
$f[i]$为$Fibonacci$的第$i$项，求$$\prod_{i=1}^n\prod_{j=1}^mf[\gcd(i,j)]\mod (10^9+7)$$ ...
分类：
其他好文 时间：
2018-04-04 21:00:38
阅读次数：
142
4816: [Sdoi2017]数字表格 Description Doris刚刚学习了fibonacci数列。用f[i]表示数列的第i项，那么 f[0]=0 f[1]=1 f[n]=f[n-1]+f[n-2],n>=2 Doris用老师的超级计算机生成了一个n×m的表格，第i行第j列的格子中的数是f ...
分类：
其他好文 时间：
2018-04-04 20:53:26
阅读次数：
154
