51nod 1188 最大公约数之和 V2

Posted lokiii

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了51nod 1188 最大公约数之和 V2相关的知识,希望对你有一定的参考价值。


第二个\\( O(T\\sqrt(n)) \\)复杂度T了..T了..T了...天地良心,这能差多少?!
于是跑去现算(。

\\[\\sum_{i=1}^{n-1}\\sum_{j=i+1}^{n}gcd(i,j) \\]

\\[\\sum_{d=1}^{n}d\\sum_{i=1}^{n-1}\\sum_{j=i+1}^{n}[gcd(i,j)==d] \\]

\\[\\sum_{d=1}^{n}d(\\sum_{j=1}^{n}\\sum_{i=1}^{j}[gcd(i,j)==d]-\\sum_{j=1}^{n}[gcd(j,j)==d]) \\]

\\[\\sum_{d=1}^{n}d(\\sum_{j=1}^{\\left \\lfloor \\frac{n}{d} \\right \\rfloor}\\sum_{i=1}^{j}[gcd(i,j)==1]-1) \\]

\\[\\sum_{d=1}^{n}d(\\sum_{j=1}^{\\left \\lfloor \\frac{n}{d} \\right \\rfloor}\\phi(j)-1) \\]

然后与\\( O(nlnn) \\)处理出所有答案。

#include<iostream>
#include<cstdio>
using namespace std;
const int N=5000005,m=5000000;
int T,n,phi[N],q[N],tot;
long long s[N],ans[N],con;
bool v[N];
int main()
{
	phi[1]=1;
	for(int i=2;i<=m;i++)
	{
		if(!v[i])
		{
			q[++tot]=i;
			phi[i]=i-1;
		}
		for(int j=1;j<=tot&&i*q[j]<=m;j++)
		{
			int k=i*q[j];
			v[k]=1;
			if(i%q[j]==0)
			{
				phi[k]=phi[i]*q[j];
				break;
			}
			phi[k]=phi[i]*(q[j]-1);
		}
	}
	for(int i=1;i<=m;i++)
		for(int j=2;j<=m/i;j++)
			ans[i*j]+=phi[j]*i;
	for(int i=1;i<=m;i++)
		ans[i]+=ans[i-1];
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d",&n);
		printf("%lld\\n",ans[n]);
	}
	return 0;
}

以上是关于51nod 1188 最大公约数之和 V2的主要内容,如果未能解决你的问题,请参考以下文章

51nod 1188 最大公约数之和 V2

51Nod 最大公约数之和V1,V2,V3;最小公倍数之和V1,V2,V3

1188 最大公约数之和 V2

51nod 1096 距离之和最小 1108 距离之和最小 V2

51nod 1190 最小公倍数之和 V2

51nod_1040:最大公约数之和