AtCoder Beginner Contest 172 D Sum of Divisors 与约数个数有关的公式推导

AtCoder Beginner Contest  172   比赛人数10142  比赛开始后6分钟看到A题,之后,比赛开始后8分钟看到所有题

AtCoder Beginner Contest   172   D   Sum of Divisors    与约数个数有关的公式推导

总目录详见https://blog.csdn.net/mrcrack/article/details/104454762

在线测评地址https://atcoder.jp/contests/abc172/tasks/abc172_d

题目大意：给定数据n,找到1的约数个数f(1),找到2的约数个数f(2),找到3的约数个数f(3),......,找到n的约数个数f(n),计算1*f(1)+2*f(2)+3*f(3)+......+n*f(n)的值.

样例模拟如下：

4

23

f(1)=1的约数个数=(约数1(共计1个))
f(2)=2的约数个数=(约数1,约数2(共计2个))
f(3)=3的约数个数=(约数1,约数2,约数3(共计3个))
f(4)=4的约数个数=(约数1,约数2,约数4(共计4个))

1*f(1)+2*f(2)+3*f(3)+4*f(4)
=(1+2+3+4)*约数1的个数1+(2+4)*约数2的个数1+3*约数3的个数1+4*约数4的个数1
=10*1+6*1+3*1+4*1
=23

(1+2+3+4)*约数1的个数1
(1+2+3+4)是等差数列，开始位置数值是1,结束位置数值是4/1=4

(2+4)*约数2的个数1
(2+4)是等差数列，开始位置数值是2,结束位置数值是2+(4/2-1)*2=4

3*约数3的个数1
3是等差数列，开始位置数值是3,结束位置数值是3+(4/3-1)*3=3

4是等差数列，开始位置数值是4,结束位置数值是4+(4/4-1)*4=4

AC代码如下：

#include <stdio.h>
#define LL long long
LL ans;
int main(){
	int n,i,m,a,b;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		m=n/i;//等差数列元素个数
		a=i,b=a+(m-1)*i;//a等差数列起始元素，b等差数列末了元素
		ans+=(LL)(a+b)*m/2;//等差数列求和
	}
	printf("%lld\n",ans);
	return 0;
}

 

以下为比赛时提交的代码，留存此处。

#include <stdio.h>
#define LL long long
LL ans,b;
int n,prime[2000],not_prime[10010],tot;
void linear_shaker(int x){
	int i,j;
	for(i=2;i<=x;i++){
		if(!not_prime[i])prime[++tot]=i;
		for(j=1;prime[j]*i<=x;j++){
			not_prime[prime[j]*i]=1;
			if(i%prime[j]==0)break;
		}
	}
}
int main(){
	int i,n,cnt,x,j;
	linear_shaker(10000);
	ans=0;
	scanf("%d",&n);
	for(i=1;i<=n;i++){
		x=i,b=1;
		for(j=1;prime[j]*prime[j]<=x;j++){
			cnt=0;
			if(x%prime[j]==0){
				while(x%prime[j]==0)x/=prime[j],cnt++;
				b*=(cnt+1);
			}
		}
		if(x>1)b*=2;
		ans+=b*i;
	}
	printf("%lld\n",ans);
	return 0;
}