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

最新推荐文章于 2021-09-20 12:22:00 发布

原创最新推荐文章于 2021-09-20 12:22:00 发布 · 1.2k 阅读

0 ·

CC 4.0 BY-SA版权

CF 同时被 2 个专栏收录

197 篇文章

订阅专栏

暴力

102 篇文章

订阅专栏

本文介绍了一种算法，用于找出不超过10^12的大数的所有因子中，不能整除任何平方数的最大因子。通过将该数分解为素数因子，并确保因子不包含重复的素数因子来实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：给一个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;
 
}