B - 完美的数字

Bob是个很喜欢数字的孩子，现在他正在研究一个与数字相关的题目，我们知道一个数字的完美度是 把这个数字分解成三个整数相乘A*A*B（0<A<=B）的方法数，例如数字80可以分解成1*1*80，2*2*20 ，4*4*5，所以80的完美度是3；数字5只有一种分解方法1*1*5，所以完美度是1，假设数字x的完美度为d(x)，现在给定a，b（a<=b），请你帮Bob求出

S，S表示的是从a到b的所有数字的流行度之和，即S=d(a)+d(a+1)+…+d(b)。

Input

输入两个整数a，b（1<=a<=b<=10^15）

Output

输出一个整数，表示从a到b的所有数字流行度之和。

Sample Input

1 80

Sample Output

107  

#include<cstdio>
#include<math.h>
using namespace std;
int main()
{
    long long a,b,sum1,sum2;
    while(~scanf("%I64d%I64d",&a,&b))
    {
        sum1=0,sum2=0;
        for(long long i=1;i<=(long long)pow(b,1.0/3)+1;i++)//之前i是定义为int，而pow前也没有加long long，结果一直报错，我搞了好久都
        {
            if(i<=b/(i*i))sum1+=b/(i*i)-i+1;//没弄明白咋回事，但在自己的电脑上运行又很正常。以后的注意不然会帆船；oj上比较死板；
        }
        for(long long i=1;i<=(long long)pow(a-1,1.0/3)+1;i++)//而且在用pow开三次方时要加1才正确，我也不明白原因；
        {
            if(i<=(a-1)/(i*i))sum2+=(a-1)/(i*i)-i+1;
        }
        printf("%I64d\n",sum1-sum2);
    }
    return 0;
}
之前觉的数据范围很大，一个一个求肯定会超时，只有找规律，通过数学关系大范围计算，但我当时只局限于单个元素之间的关系，结果找了半下午都每找到。后来经学长点醒，从区间上下手，然后发现了快速计算的方法；