codeforces#326-B-Duff in Love-暴力

题意：给一个n 【1，10^12】

求出所有因子中，  不能整除一个平方数 的 最大因子；

直接把n 拆分成素数，答案就是不重复的因子的乘积； （只要重复了，也就是有平方数为因子，不合法 ）

当时直接把每个因子去试除。。。暴力判断。。最后 少写了一个控制边界的条件。。。systemtest时导致re了。。23333

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <iostream>
#include <queue>
#include <map>
#include <set>
#include <vector>
using namespace std;
const __int64 inf=2147483647;
__int64 tm[2000005]; 
__int64 prime[2000005];
__int64 f[2000005];
__int64 cmp(__int64 a,__int64 b)
{
    return a>b;
}
int main()
{
    __int64 n;
    __int64 i,j;
scanf("%I64d",&n);
        f[1]=true;
        for (i=2;i<=2000000;i++)    //1处
    {
        if (f[i]==false)           //优化
        {
            for (j=(__int64  )i*i;j<=2000000;j=j+i)   // 如果1处用1000必须___int6464强制转换
            {
                f[j]=true;  
            } 
        }
    } 
        __int64 ok=0; 
        for (i=1;i<=2000000;i++)
        {
            if (f[i]==false)
                prime[++ok]=i;
        }
__int64 tol_prime=ok;
__int64 len=sqrt(n);
ok=0;
for (i=1;i<=len;i++)
{
        if (n%i==0)
        {
            tm[++ok]=i;
            tm[++ok]=n/i;
        }
}
//n; 
sort(tm+1,tm+1+2*ok,cmp);

for (i=1;i<=ok*2;i++)
{
    __int64 tmp=tm[i];
    j=1;
    __int64 flag=0;
    while(tmp&&prime[j]<tmp&&j<=tol_prime)<span style="white-space:pre">	</span>//少写了j<=tot_prime导致re
    {
			 
        if (tmp%prime[j]==0)
        { 
            if ((tmp/prime[j])%prime[j]==0)
            {flag=1;break;}
            else
                j++;
        }
        else
            j++;
    }
    if (flag)
    continue;
    else
    break ;
} 
    printf("%I64d\n",tm[i]);
 








    return 0;
 
}