2021-01-31

最新推荐文章于 2024-06-01 21:43:54 发布

提莫队长丶

最新推荐文章于 2024-06-01 21:43:54 发布

阅读量50

点赞数

本文链接：https://blog.csdn.net/qq_51920503/article/details/113481607

版权

素数溢出平方根算法优化数论

关键词由CSDN通过智能技术生成

hd3826

#include <stdio.h>
#include <math.h>
#define N 1000000
typedef long long ll;
int vit[N+1]= {1,1},prime[1000000],count=0;
void creat_prime()
{
    int m=sqrt(N+0.5);
    for(int i=2; i<=m; i++)
        if(vit[i]==0)
            for(int j=i*i; j<=N; j+=i)
                vit[j]=1;
    for(int i=2; i<=N; i++)
        if(vit[i]==0)
            prime[count++]=i;
}
int main()
{
    ll t, n;
    int l=1;
    creat_prime();
    scanf("%lld", &t);
    for (int i = 0; i < t; i++)
    {
        int flag;
        scanf("%lld", &n);
        int m=sqrt(n+0.5);
        for (int j = 0; j<count&&prime[j]<=m; j++)
        {
            int k = 0;
            flag=0;
                while (n % prime[j] == 0)
                {
                    n = n / prime[j];
                    k++;
                    if (k >= 2)
                    {
                        flag = 1;
                        break;
                    }
                }
             if (flag == 1)
                break;
        }
        if (flag ==1 )
            printf("Case %d: No\n",l++);
        else if(flag==0&&n>1)
        {
            double a=sqrt(n*1.0);
            if(a-floor(a)<0.000001)
                printf("Case %d: No\n",l++);
            else
                printf("Case %d: Yes\n",l++);
        }
        else
            printf("Case %d: Yes\n",l++);
            flag=0;
    }
    return 0;
}

问题一；下面a*a=n会发生溢出

#include <stdio.h>
#include <math.h>
#define N 1000000
typedef long long ll;
int vit[N+1]= {1,1},prime[1000000],count=0;
void creat_prime()
{
    int m=sqrt(N+0.5);
    for(int i=2; i<=m; i++)
        if(vit[i]==0)
            for(int j=i*i; j<=N; j+=i)
                vit[j]=1;
    for(int i=2; i<=N; i++)
        if(vit[i]==0)
            prime[count++]=i;
}
int main()
{
    ll t, n;
    int l=1;
    creat_prime();
    scanf("%lld", &t);
    for (int i = 0; i < t; i++)
    {
        int flag;
        scanf("%lld", &n);
        int m=sqrt(n+0.5);
        for (int j = 0; j<count&&prime[j]<=m; j++)
        {
            int k = 0;
            flag=0;
                while (n % prime[j] == 0)
                {
                    n = n / prime[j];
                    k++;
                    if (k >= 2)
                    {
                        flag = 1;
                        break;
                    }
                }
             if (flag == 1)
                break;
        }
        if (flag ==1 )
            printf("Case %d: No\n",l++);
        else if(flag==0&&n>1)
        {
            int a=sqrt(n+0.5);
            if(a*a==n)
                printf("Case %d: No\n",l++);
            else
                printf("Case %d: Yes\n",l++);
        }
        else
            printf("Case %d: Yes\n",l++);
            flag=0;
    }
    return 0;
}