码迷,mamicode.com
首页 > 其他好文 > 详细

hunnu11546:Sum of f(x)

时间:2017-04-30 17:12:07      阅读:103      评论:0      收藏:0      [点我收藏+]

标签:name   ret   eof   queue   mil   init   printf   names   inf   

Problem description
  令f(x)为x的全部约数之和,x的约数即能够被x整除的数。如f(24)=1+2+3+4+6+8+12+24=60),求 f(l) + f(l + 1) + …… + f(r)
Input
  第一行为一个整数T(T<=100000),表示数据的组数。


接下来T行,每行有两个整数l。r(1 <= l <= r <= 200000)

Output
  对每组数据,输出f(l)+f(l+1)+……+f(r) 的和
Sample Input
2
3 5
10 20
Sample Output
17
270
Problem Source
  HUNNU Contest 


開始用这样的方法来求

 for(i = 3; i<=200000; i++)
    {
        a[i] = 1+i;
        for(j = 2; j*j<=i; j++)
        {
            if(i%j==0)
            {
                a[i]+=j;
                if(i/j!=j)
                    a[i]+=(i/j);
            }
        }
    }

超时。改成

    for(i = 1; i<=200000; i++)
    {
        for(j = 1; i*j<=200000; j++)
        {
            a[i*j]+=i;
        }
    }

才AC

开了一个变量去级数,发现上面的运算了5千多W次

以下的才两百多W次

看来区别还是蛮大的


#include <iostream>
#include <stdio.h>
#include <string.h>
#include <stack>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <math.h>
#include <bitset>
#include <algorithm>
#include <climits>
#include <time.h>
using namespace std;

#define LS 2*i
#define RS 2*i+1
#define UP(i,x,y) for(i=x;i<=y;i++)
#define DOWN(i,x,y) for(i=x;i>=y;i--)
#define MEM(a,x) memset(a,x,sizeof(a))
#define W(a) while(a)
#define gcd(a,b) __gcd(a,b)
#define LL long long
#define N 200005
#define MOD 1000000007
#define INF 0x3f3f3f3f
#define EXP 1e-8

LL a[N];
LL sum[N];

void init()
{
    LL i,j,k;
    MEM(a,0);
    MEM(sum,0);
    for(i = 1; i<=200000; i++)
    {
        for(j = 1; i*j<=200000; j++)
        {
            a[i*j]+=i;
        }
    }
    for(i = 1; i<=200000; i++)
        sum[i] = sum[i-1]+a[i];
}
int main()
{
    LL i,j,k;
    int l,r;
    init();
    int t;
    scanf("%d",&t);
    while(t--)
    {
       scanf("%d%d",&l,&r);
        printf("%I64d\n",sum[r]-sum[l-1]);
    }

    return 0;
}


hunnu11546:Sum of f(x)

标签:name   ret   eof   queue   mil   init   printf   names   inf   

原文地址:http://www.cnblogs.com/yjbjingcha/p/6789624.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!