【bzoj2693】jzptab 莫比乌斯反演+线性筛

题目描述

输入

一个正整数T表示数据组数

接下来T行每行两个正整数表示N、M

输出

T行每行一个整数表示第i组数据的结果

样例输入

1
4 5

样例输出

122

题解

莫比乌斯反演+线性筛

由于要处理多组询问，所以 bzoj2154 的做法就不好用了，但是这个结论可以套用过来。

然后推公式：

（UPD：上面公式最后一行请自行把 $k$ 改成 $n$ ... 由于这里是图片形式就不改了）

设f1(n)=n²mu(n)，f2(n)=n，则显然f2是积性函数，f1为两个积性函数的乘积，也是积性函数。

那么f为f1和f2的狄利克雷卷积，也是积性函数。

所以可以尝试快筛f(n)。

当n为质数时，显然f(n)=n-n^2。

当n不为质数时，即n=i*p，p|i，p是质数，那么观察f(n)化简之后的式子，n新增加出来的约数一定是包含p^2的，它的mu值一定是0，所以f(n)的改变只是从i*...变为了n*...，所以此时f(n)=f(i)*p。

这样我们就可以快筛出f函数及其前缀和，然后对于每个询问分块求解即可。

#include <cstdio>
#include <algorithm>
#define mod 100000009
using namespace std;
const int n = 10000000;
typedef long long ll;
int prime[n + 10] , tot;
ll g[n + 10] , sum[n + 10];
bool np[n + 10];
ll s(int x)
{return (ll)x * (x + 1) / 2 % mod;
}
ll cal(ll a , ll b)
{int i , last;ll ans = 0;for(i = 1 ; i <= a && i <= b ; i = last + 1) last = min(a / (a / i) , b / (b / i)) , ans = (ans + (sum[last] - sum[i - 1] + mod) % mod * s(a / i) % mod * s(b / i) % mod) % mod;return ans;
}
int main()
{int i , j , T , a , b;g[1] = sum[1] = 1;for(i = 2 ; i <= n ; i ++ ){if(!np[i]) g[i] = ((i - (ll)i * i) % mod + mod) % mod , prime[++tot] = i;for(j = 1 ; j <= tot && i * prime[j] <= n ; j ++ ){np[i * prime[j]] = 1;if(i % prime[j] == 0){g[i * prime[j]] = g[i] * prime[j] % mod;break;}else g[i * prime[j]] = g[i] * g[prime[j]] % mod;}sum[i] = (sum[i - 1] + g[i]) % mod;}scanf("%d" , &T);while(T -- ) scanf("%d%d" , &a , &b) , printf("%lld\n" , cal(a , b));return 0;
}

转载于:https://www.cnblogs.com/GXZlegend/p/7000042.html

【bzoj2693】jzptab 莫比乌斯反演+线性筛相关推荐

【bzoj2694】Lcm 莫比乌斯反演+线性筛
题目描述求$\sum\limits_{i=1}^n\sum\limits_{j=1}^m|\mu(gcd(i,j))|lcm(i,j)$,即$gcd(i,j)$不存在平方因子的$lcm(i,j)$之 ...
BZOJ3309 DZY Loves Math（莫比乌斯反演+线性筛）
一通正常的莫比乌斯反演后,我们只需要求出g(n)=Σf(d)*μ(n/d)的前缀和就好了. 考虑怎么求g(n).当然是打表啊.设n=∏piai,n/d=∏pibi .显然若存在bi>1则这个d没 ...
bzoj4407 于神之怒加强版（莫比乌斯反演+线性筛）
4407: 于神之怒加强版 Time Limit: 80 Sec Memory Limit: 512 MB Submit: 355 Solved: 174 [Submit][Status][Dis ...
【bzoj3309】DZY Loves Math 莫比乌斯反演+线性筛
Description 对于正整数n,定义f(n)为n所含质因子的最大幂指数.例如f(1960)=f(2^3 * 5^1 * 7^2)=3, f(10007)=1, f(1)=0. 给定正整数a,b, ...
bzoj 2820 YY的GCD - 莫比乌斯反演 - 线性筛
Description 神犇YY虐完数论后给傻×kAc出了一题给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对kAc这种傻×必 ...
Bzoj3309-DZY Loves Math【莫比乌斯反演,线性筛】
正题 bzoj没了,在darkbzoj交吧题目链接:https://darkbzoj.tk/problem/3309 题目大意定义f(x)f(x)f(x)表示xxx所有质因数中最大的指数幂. 求∑ ...
南昌邀请赛网络赛 G. tsy's number（莫比乌斯反演+线性筛）
题目链接: tsy's number 题意: 求 (mod ) 思路: 原式 = 设设原式 = 设T = id ,将 i 的求和变为 T 的求和: 原式 = :这个是与的 ...
洛谷P2257 YY的GCD 莫比乌斯函数反演+线性筛
洛谷P2257 YY的GCD 标签莫比乌斯反演线性筛前言这题貌似和莫反没多大关系,就是用到了一个莫比乌斯函数的性质了,其他就是推公式,优化和式. 我的第一道懵逼反演-真的好难好难-而且套路特别 ...
【专题】莫比乌斯反演
问题引入: 给定整数 N 和 M.求满足且为质数的点对的个数. 数据范围: 接下来会见到以下内容: 莫比乌斯函数莫比乌斯函数的线性筛迪利克雷卷积介绍莫比乌斯反演整除分块杜教筛介绍 ...

【bzoj2693】jzptab 莫比乌斯反演+线性筛

一个正整数T表示数据组数

【bzoj2693】jzptab 莫比乌斯反演+线性筛相关推荐

最新文章

热门文章