Description

给出一个数字N

Input

第一行为一个正整数T,表示数据组数。
接下来T行为询问,每行包含一个正整数N。
T<=5000,N<=10^7
题解: 求 $\sum_{i=1}^{n}\sum_{j=1}^{n}\varphi(gcd(i,j))$
 
$=\sum_{d=1}^{n}\varphi(d)\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd(i,j)==d]$
 
$=\sum_{d=1}^{n}\varphi(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{n}{d}}[gcd(i,j)==1]$
 
其中 $\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd(i,j)==1]=2\times\sum_{i=1}^{n}\varphi(i)-1$
 
可以用数形结合来思考:
 
$\sum_{i=1}^{n}\sum_{j=1}^{i}[gcd(i,j)==1]=\sum_{i=1}^{n}\varphi(i)$
 
这个式子可以看作是对角线一半的欧拉函数值(以 $i$ 为横坐标,$j$ 为纵坐标)
 
而 $\sum_{i=1}^{n}\sum_{j=1}^{n}[gcd(i,j)==1]$ 正好是对角线折叠过去的结果.
 
不过还要减去 $(1,1)$ 多出来的贡献(只有 $1$ 是和自己本身互质的数)
 
令 $Sum(n)=\sum_{i=1}^{n}\varphi(i)$
 
答案为 $\sum_{d=1}^{n}\varphi(d)\sum_{i=1}^{\frac{n}{d}}\sum_{j=1}^{\frac{n}{d}}[gcd(i,j)==1]$
 
$\Rightarrow\sum_{d=1}^{n}\varphi(d)\times(2\times\sum_{i=1}^{\frac{n}{d}}\varphi(i)-1)$
 
$\Rightarrow\sum_{d=1}^{n}\varphi(d)\times(2\times Sum(\frac{n}{d})-1)$
 
$\Rightarrow2\times\sum_{d=1}^{n}\varphi(d)Sum(\frac{n}{d})-Sum(n)$
 
提前与处理一下欧拉函数前缀和即可.
#include<bits/stdc++.h>
#define maxn 10000003
#define M 10000000
#define ll long long
using namespace std;
void setIO(string s)
{
string in=s+".in";
freopen(in.c_str(),"r",stdin);
}
int cnt;
bool vis[maxn];
int prime[maxn], phi[maxn];
ll sumv[maxn];
inline void Linear_shaker()
{
phi[1]=1;
int i,j;
for(i=2;i<=M;++i)
{
if(!vis[i]) prime[++cnt]=i, phi[i]=i-1;
for(j=1;j<=cnt&&1ll*i*prime[j]<=M;++j)
{
vis[i*prime[j]]=1;
if(i%prime[j]==0)
{
phi[i*prime[j]]=phi[i]*prime[j];
break;
}
phi[i*prime[j]]=phi[i]*(prime[j]-1);
}
}
for(i=1;i<=M;++i) sumv[i]=sumv[i-1]+phi[i];
}
inline void solve()
{
int n,i,j;
scanf("%d",&n);
ll re=0;
for(i=1;i<=n;i=j+1)
{
j=(n/(n/i));
re+=(sumv[j]-sumv[i-1])*(sumv[n/i]*2-1);
}
printf("%lld\n",re);
}
int main()
{
// setIO("input");
Linear_shaker();
int T;
scanf("%d",&T);
while(T--) solve();
return 0;
}

  

最新文章

  1. Arduino下LCD1602综合探究(上)——1602的两种驱动方式,如何使LCD的控制编程变得更简单
  2. Git学习笔记(10)——搭建Git服务器
  3. 黄聪:PHP页面跳转如何实现延时跳转
  4. 构建VIM下的C++编程环境
  5. sed中求公共前缀
  6. Python dir()函数
  7. CSS.01 -- 选择器及相关的属性文本、文字、字体、颜色、
  8. Android开发学习之路--Broadcast Receiver之初体验
  9. python_unittest学习小结
  10. [Swift]LeetCode804. 唯一摩尔斯密码词 | Unique Morse Code Words
  11. 用word发CSDN blog,免去插图片的烦恼
  12. vue中原生file上传图片
  13. [UnityShader基础]05.模板测试
  14. kubernetes集群搭建(4):node节点安装
  15. .NET持续集成与自动化部署之路第三篇——测试环境到生产环境的一键部署策略(Windows)
  16. openpyxl 实现excel字母列号与数字列号之间的转换
  17. BZOJ 4806 - 4809 象棋四题
  18. MyBatis #{} 和 ${} 引用值的用法
  19. window下安装mongodb3.6
  20. SQL 文件以及文件组

热门文章

  1. VMware 虚拟化编程(6) — VixDiskLib 虚拟磁盘库详解之二
  2. Linux负载均衡实现
  3. jmeter逻辑控制详解(1)
  4. 关于Vue的理解以及与React框架的对比
  5. python+selenium的frame表单切换
  6. 红黑树的删除操作---以JDK源码为例
  7. Tornado 的核是什么??
  8. P1828 香甜的黄油 (spfa)
  9. git常用相关操作
  10. ASE Alpha Sprint - backend scrum 5