题目描述
科学家们在 Samuel 星球上的探险得到了丰富的能源储备，这使得空间站中大型计算机 Samuel II 的长时间运算成为了可能。由于在去年一年的辛苦工作取得了不错的成绩，小联被允许用 Samuel II 进行数学研究。
小联最近在研究和约数有关的问题，他统计每个正数 N 的约数的个数，并以 f(N) 来表示。例如 12 的约数有 1,2,3,4,6,12，因此 f(12)=6。下表给出了一些 f(N) 的取值：
| N | 1 | 2 | 3 | 4 | 5 | 6 |
|---|---|---|---|---|---|---|
| f(N) | 1 | 2 | 2 | 3 | 2 | 4 |
现在请你求出：
i= ∑(1到n) f(i)
输入格式
输入一个整数 n。
输出格式
输出答案。
输入输出样例
输入 #1复制
3
输出 #1复制
5
说明/提示
- 对于 20% 的数据，N≤5000；
- 对于 100% 的数据，1≤N≤10^6。
根据约数的一个基本性质，每个正约数 在 1 ~ n 中出现的次数为
，
so问题转换成求 的值,套整除分块板子；
#include<bits/stdc++.h>
using namespace std;
//#pragma GCC optimize(2)
#define endl '\n'
#define lowbit(x) ((x)&-(x))
const int mod=998244353;
typedef long long ll;
ll ans=0,n1,m1;
ll t,s1,s2,s3,s4,max1=0,min1=100000000,sum=0,n,m,i,j,k,l,r;
ll u[200005],w,v[200005];
inline int read() {bool sym=0;int res=0;char ch=getchar();while(!isdigit(ch))sym |=(ch =='-'),ch=getchar();while(isdigit(ch)) res =(res<<3)+(res<<1)+(ch^48),ch=getchar();return sym ? -res : res;
}
void print(int x) {if(!x)return;print(x/10);putchar(x%10+'0');
}
int isPrime(int n) {float n_sqrt;if(n==1) return 0;if(n==2 || n==3) return 1;if(n%6!=1 && n%6!=5) return 0;n_sqrt=floor(sqrt((float)n));for(int i=5; i<=n_sqrt; i+=6) {if(n%(i)==0 | n%(i+2)==0) return 0;}return 1;
}/*ll cal(ll n) {ll res = 0;for(ll l = 1, r; l <= n; l = r + 1) {r = n / (n / l);res += n / l *(r-l+1)+ (l+r) * (r - l + 1)/2 ;// 等差数列 求和}return res;
}*/int main() {ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);cin>>n;for(l=1; l<=n; l=r+1) {r=n/(n/l);ans+=n/l*(r-l+1);}cout <<ans<< endl;return 0;
}//mio lover
