CodeForces 172 D.Calendar Reform（数论）

Description
第一年a天，第二年a+1天，以此类推，第n年a+n-1天，一年可以被分成几个月，一个月的天数必须是完全平方数，一年中每个月的天数相同，问这n天最多有几个月
Input
两个整数a和n(1<=a,n<=1e7,a+n-1<=1e7)
Output
输出最少月数
Sample Input
25 3
Sample Output
30
Solution
问题转化为把a~a+n-1中的每个数拆成a^2*b形式，且b最小，所以把每个数的平方因子全部除掉即可
Code

#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
const int maxn=10000001;
int a,n,x[maxn];
int main()
{
    scanf("%d%d",&a,&n);
    for(int i=a;i<=a+n-1;i++)x[i]=i;
    for(int i=2;i*i<=a+n-1;i++)
    {
        int p=i*i;
        int t=((a+p-1)/p)*p;
        for(int j=t;j<=a+n-1;j+=p)
            while(x[j]%p==0)x[j]/=p;
    }
    ll ans=0;
    for(int i=a;i<=a+n-1;i++)ans+=x[i];
    printf("%I64d\n",ans);
    return 0;
}