n的约数和以及前n个数的约数和

• 求n的正约数和
  唯一分解定理：
  $N=p_1^{a_1}\ast p_2^{a_1}\ast ...\ast p_n^{a_n}\ast$
  ${\sigma }_1=(1+p_1+p_1^2+...+p_1^{a_1})\ast (1+p_2+p_2^2+...+p_2^{a_2})\ast ...\ast (1+p_n+p_n^2+...+p_n^{a_n})$

故直接写就好了。

• 求前n的个数的约数和

参考博客：
神犇1

约数和函数 ：${\sigma }_1(n)={\sum }_{d|n}d$
现在求 ${\sum }_{i=1}^n={\sigma }_1(i)$

${\sum }_{i=1}^n={\sigma }_1(i)={\sum }_{i=1}^{n}i\ast \lfloor \frac{n}{i}\rfloor$

代码1：

#include<bits/stdc++.h>
using namespace std;

typedef long long LL;

int main()
{
    
    LL n;
    scanf("%lld",&n);

    int num=0;
    LL ed=1,sum=0;
    
    ///整数分块
    for(LL a=1;a<=n;a=ed+1)
    {
        ed=(n/(n/a));
        sum+=(a+ed)*(ed-a+1)/2*(n/a);
//        num++;
    }
    printf("%lld\n",sum);
}

神犇2
代码2参考上面的链接，写这个是为了以后自己好回顾：

#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std;

typedef long long LL:

LL n,i,a,b;
LL sum;

int main(){
	while(~scanf("%lld",&n)){
	sum = 0;
	for(i=1;i*i<=n;i++){
            ///计算前sqrt(n)的因子i贡献
		sum += i * (n / i);

       ///另一个因子n/i，j属于区间[a,b]内，n/j都相同，故直接算它们的和
		b = n/i;
		a = n/(i+1)+1;
//		printf("a=%lld b=%lld\n",a,b);
		///区间[a,b]和的贡献为它们的i倍
		sum += (a+b) * (b-a+1) / 2 * i;
	}
	i--;
	if (i == n / i)
		sum -= i * (n / i);
	printf("%lld\n",sum);
	}
	return 0;
}