标签：style class blog code width 2014
因子和 |
||
| Accepted : 47 | Submit : 289 | |
| Time Limit : 4000 MS | Memory Limit : 65536 KB | |
题目描述如果b能整除a，我们称b为a的因子。现在假设n的所有因子和为f(n)； 输入不超过2000个样例，每个样例占一行，为两个整数a，b；a=b=0时表示输入结束。 输出每个样例输出一个整数，占一行。 样例输入1 1 6 6 1 6 0 0 样例输出1 12 33 |
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#define ll __int64
using namespace std;
ll a[5000001];
int main()
{
int n,m;
ll int sum;
int i=0,j=0;
a[0] = 0;
int e = (int)(sqrt(0.0 + 5000000)+1);
for(i=2;i<=e;i++)
{
a[i*i] += i;
for(j=i+1;i*j<=5000000;j++)
{
a[i*j] += i+j;
}
}
while(scanf("%d%d",&n,&m)!=EOF)
{
if(n==0 && m== 0)
break;
sum = 0;
for(int i = n;i<=m;i++) <- 这里浪费时间较多，要重新再遍历一次
{
sum+= a[i];
if(i>1)
sum += i;
}
if(n==0)
cout<<sum+m-n;
else
cout<<sum+m-n+1;
}
return 0;
}#include <cstdio>
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cmath>
#define N 5000000
#define ll __int64
using namespace std;
ll a[N+1]={0};
int m,n;
int main()
{
a[0]=0; a[1]=1;
for(int i=2; i<=N; i++)
a[i]=1+i;
for(int i=2; i<=N; i++) //以i来暴力枚举当前数的所有可能因子
{
a[i] += a[i-1];
//往后覆盖i的2倍的公差为j的元素
for(int j=i*2; j<=N;j+=i) //<-熟悉吧 ，纯素数筛
a[j]+=i;
}
while(cin>>n>>m)
{
if(m==0 && n==0)
break;
if(n==0)
cout<<a[m]<<endl;
else
cout<<a[m]-a[n-1]<<endl;
}
return 0;
}
标签：style class blog code width 2014
原文地址：http://blog.csdn.net/wjw0130/article/details/30794997
