Description

Input

题解:

这次是多组询问,上次推出来的式子是 $O(n)$ 的,我们需要更快的算法.
 
先定义几个后面可能会用到的函数:(在本题弱化版中推出来的)
 
$Sum(n,m)=\frac{n(n+1)}{2}\times \frac{m(m+1)}{2}$
 
$calc(n,m)=\sum_{d=1}^{n}\mu(d)d^2Sum(\frac{n}{d},\frac{m}{d})$
 
则 $ans=\sum_{d=1}^{n}d\times calc(\frac{n}{d},\frac{m}{d})$
 
把式子展开:
 
$=\sum_{d=1}^{n}d\sum_{b=1}^{\frac{n}{d}}\mu(b)b^2Sum(\frac{n}{db},\frac{m}{db})$
 
令 $D=db$ ,得:
 
$\sum_{D=1}^{n}Sum(\frac{n}{D},\frac{m}{D})\sum_{d|D}d\times\mu(\frac{D}{d})\times(\frac{D}{d})^2$ 
 
看上去有点别扭,再稍微改动一下:
 
$\sum_{D=1}^{n}Sum(\frac{n}{D},\frac{m}{D})\sum_{d|D}\frac{D}{d}\times\mu(d)\times d^2$ 
 
前面 $\sum_{D=1}^{n}Sum(\frac{n}{D},\frac{m}{D})$ 这个部分好像没法优化了,只能$\sqrt{n}$ 硬算.
 
好在 $f(D)=\sum_{d|D}\frac{D}{d}\times\mu(d)\times d^2$ 是一个积性函数,可以线筛.
 
只需分三种情况讨论即可:
 
(1) $x$ 为质数.$\Rightarrow$ $f(x)=x^{2}-x$ (只有 $d=1,d=D$ 时有贡献)
 
(2) $a,b$ 互质.$\Rightarrow$ $f(ab)=f(a)f(b)$ (积性函数定义)
 
(3) $p|i$ ,且 $p$ 为质数, $i$ 中 $p$ 的指数为 $1$
 
$f(i)=\sum_{d|i}\frac{i}{d}\times\mu(d)\times d^2$
 
$f(i\times p)=\sum_{d|i\times p}\frac{i\times p}{d}\times\mu(d)\times d^2$
 
由 $\mu$ 的定义可知,当一个数有质因子次幂大于 $1$ 时值为 $0$ ,即无贡献.
 
那么,$f(i\times p)$ ,$d$ 中 $p$ 的幂次必小于等于 $1$.
 
而 $f(i)$ 中所有 $d$ 中 $p$ 的次幂本来就小于等于 $1$.
 
故 $f(i)$ 与 $f(i\times p)$ 中有贡献的因子是完全相同的.
 
唯一不同就是 $f(i\times p)=f(i)\times p$
 
既然讨论好上面 $3$ 种情况,我们就可以愉快地线筛 $f(x)$ 啦!!
  1. #include<bits/stdc++.h>
  2. #define ll long long
  3. #define M 10001000
  4. #define maxn 10200100
  5. #define MOD 100000009
  6. using namespace std;
  7. int cnt, tot;
  8. int vis[maxn],mu[maxn], prime[maxn];
  9. ll h[maxn], sumv[maxn];
  10. void init()
  11. {
  12. int i,j;
  13. h[1]=1;
  14. for(i=2;i<M;++i)
  15. {
  16. if(!vis[i]) prime[++tot]=i, h[i]=(i-(ll)i*i)%MOD;
  17. for(j=1;j<=tot&&prime[j]*i<M;++j)
  18. {
  19. vis[prime[j]*i]=1;
  20. if(i%prime[j]==0)
  21. {
  22. h[prime[j]*i]=(prime[j]*h[i])%MOD;
  23. break;
  24. }
  25. else h[prime[j]*i]=(h[prime[j]]*h[i])%MOD;
  26. }
  27. }
  28. for(i=1;i<M;++i) sumv[i]=(sumv[i-1]+h[i])%MOD;
  29. }
  30. ll SUM(ll x,ll y)
  31. {
  32. x%=MOD, y%=MOD;
  33. ll r1=(x*(x+1)>>1)%MOD;
  34. ll r2=(y*(y+1)>>1)%MOD;
  35. return (r1*r2)%MOD;
  36. }
  37. ll query(ll n,ll m)
  38. {
  39. int i,last,re=0,j;
  40. if(n>m) swap(n,m);
  41. for(i=1;i<=n;i=j+1)
  42. {
  43. j=min(n/(n/i), m/(m/i));
  44. re += SUM(n/i, m/i) * (sumv[j]-sumv[i-1])%MOD;
  45. re%=MOD;
  46. }
  47. return (re+MOD)%MOD;
  48. }
  49. int main()
  50. {
  51. init();
  52. int n,m,T;
  53. scanf("%d",&T);
  54. while(T--)
  55. {
  56. scanf("%d%d",&n,&m);
  57. printf("%lld\n",query(n,m));
  58. }
  59. return 0;
  60. }

  

BZOJ 2693: jzptab 莫比乌斯反演 + 积性函数 +筛法的更多相关文章

  1. bzoj 2693: jzptab 线性筛积性函数

    2693: jzptab Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 444  Solved: 174[Submit][Status][Discus ...

  2. BZOJ 2694: Lcm 莫比乌斯反演 + 积性函数 + 线性筛 + 卡常

    求 $\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)\mu(gcd(i,j))^2$   $\Rightarrow \sum_{d=1}^{n}\mu(d)^2\sum_{i ...

  3. BZOJ 2693: jzptab [莫比乌斯反演 线性筛]

    2693: jzptab Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 1194  Solved: 455[Submit][Status][Discu ...

  4. BZOJ 2693 jzptab ——莫比乌斯反演

    同BZOJ 2154 但是需要优化 $ans=\sum_{d<=n}d*\sum_{i<=\lfloor n/d \rfloor} i^2 *\mu(i)* Sum(\lfloor \fr ...

  5. BZOJ 2693: jzptab( 莫比乌斯反演 )

    速度居然#2...目测是因为我没用long long.. 求∑ lcm(i, j) (1 <= i <= n, 1 <= j <= m) 化简之后就只须求f(x) = x∑u( ...

  6. 【BZOJ】2693: jzptab 莫比乌斯反演

    [题意]2154: Crash的数字表格 莫比乌斯反演,多组询问,T<=10000. [算法]数论(莫比乌斯反演) [题解]由上一题, $ans=\sum_{g\leq min(n,m)}g\s ...

  7. bzoj2693--莫比乌斯反演+积性函数线性筛

    推导: 设d=gcd(i,j) 利用莫比乌斯函数的性质 令sum(x,y)=(x*(x+1)/2)*(y*(y+1)/2) 令T=d*t 设f(T)= T可以分块.又由于μ是积性函数,积性函数的约束和 ...

  8. $BZOJ$2818 $gcd$ 莫比乌斯反演/欧拉函数

    正解:莫比乌斯反演/欧拉函数 解题报告: 传送门$QwQ$ 一步非常显然的变形,原式=$\sum_{d=1,d\in prim}^{n}\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd ...

  9. 【莫比乌斯反演】关于Mobius反演与lcm的一些关系与问题简化(BZOJ 2154 crash的数字表格&&BZOJ 2693 jzptab)

    BZOJ 2154 crash的数字表格 Description 今天的数学课上,Crash小朋友学习了最小公倍数(Least Common Multiple).对于两个正整数a和b,LCM(a, b ...

随机推荐

  1. STM32学习之路-感觉自己走到了一个天大的坑里了!

    先前兴致勃勃的来弄16位并口驱动LCD.本以为就须要改下LCD IC的初始化即可了,没想到弄了这么多天最终发现自己走进了一个深坑了 T T 原因是我的开发板是奋斗V5的, 它确实有MCU外扩IO口, ...

  2. BZOJ 1055 区间DP

    1055: [HAOI2008]玩具取名 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 1144  Solved: 668[Submit][Statu ...

  3. 以太网接口TCP/IP协议介绍,说的很容易懂了

      以太网接口TCP/IP协议介绍,说的很容易懂了  TCP/IP协议,或称为TCP/IP协议栈,或互联网协议系列. TCP/IP协议栈(按TCP/IP参考模型划分) 应用层 FTP SMTP HTT ...

  4. Weex 和 React Native 的比较看这里

    写在前面 目前主流的应用大体分成三类:Native App, Web App, Hybrid App. Native App 特点: 性能好 完美的用户体验 开发成本高,无法跨平台 升级困难 (审核) ...

  5. BZOJ_3280_小R的烦恼_最小费用最大流

    BZOJ_3280_小R的烦恼_最小费用最大流 Description 小R最近遇上了大麻烦,他的程序设计挂科了.于是他只好找程设老师求情.善良的程设老师答应不挂他,但是要 求小R帮助他一起解决一个难 ...

  6. uva10828

    https://vjudge.net/problem/UVA-10828 裸高斯消元... 但是要判无解和无穷解. 当出现一个环时会无解,环上每个点只有一个出边. #include<bits/s ...

  7. Tomcat优化和JVM分析工具

    Tomcat的常见优化和JVM常见分析工具 Tomcat的常用优化配置 (1) 内存空间: /etc/sysconfig/tomcat JAVA_OPTS="-server -Xms32g ...

  8. 盘点国内网站常用的一些 CDN 公共库加速服务(转载)

    百度CND jQuery 地址:<script src="http://libs.baidu.com/jquery/1.9.0/jquery.js"></scri ...

  9. 【转】20道Spring Boot面试题

    面试了少量人,简历上都说自己熟习 Spring Boot, 或者者说正在学习 Spring Boot,一问他们时,都只停留在简单的使用阶段,很多东西都不清楚,也让我对面试者大失所望. 下面,我给大家总 ...

  10. golang——字符串与编码

    1.字符编码 (1)ASCII码 一个字节表示的英文.数字.标点符号等字符. 国际标准ASCII码为0-127即128个字符,二进制最高位为0,其余为扩展ASCII码. (2)GB2312 两字节,主 ...