hdu 3864 D_num

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=3864

题意很简单,就是确定一个数n是不是除了1和它本身,还只有两个不等的因子。解题的思想就是首先用费马小随机测素数算法确定这个数是不是素数,如果不是素数,用pollard随机算法找出它的一个因子p,则另一个因子q为n/p;如果p,q同时为素数,则这个n一定符合要求,否则,假设p<q,如果 p^2 = q,即n为一个素数的三次方,也同样符合要求,从小到大输出p,q,n即可。参考代码如下:

#include<iostream>
#include<cstdio>
#include<cstdlib>
using namespace std;

#define LL __int64

LL multi(LL a,LL b,LL mod)
{
    LL temp = a,sum = 0;
    while(b)
    {
        if(b&1) sum = (sum +  temp) % mod;
        temp = (temp + temp) % mod;
        b = b>>1;
    }
    return sum;
}

LL powerful(LL a,LL x,LL mod)
{
    LL t = a % mod,r = 1;
    while(x)
    {
        if(x & 1) r = multi(r,t,mod);
        t = multi(t,t,mod); x = x>>1;
    }
    return r;
}
bool MR(LL time,LL n)
{
    if(n < 2) return false;
    if(n == 2) return true;
    int num[3] = {2,7,61};
    bool flag = false;
    for(LL k = 0;k < 3;++k)
    {
        if(num[k] == n) return true;
        flag = false;
        LL d = n - 1,r = 0,t,a = rand()%(n-2) + 2;
        while((d & 1) == 0)
        {
            d = d>>1;
            r++;
        }
        t = powerful(a,d,n);
        if(t == 1 || t == n-1) {flag = true;continue;}
        for(LL i = 1;i < r;i++)
        {
            t = multi(t,t,n);
            if(t == 1) {flag = false;return  flag;}
            if(t == n-1) {flag = true;break;}
        }
        if(!flag) break;
    }
    return flag;
}

LL gcd(LL a,LL b) { return !b ? a : gcd(b,a % b); }

LL pollard(LL c,LL n)
{
    LL x = rand() % n;
    LL d,i = 1,k = 2,y = x;
    do{
        i++;
        if((x = multi(x,x,n) - c) < 0) x += n;
        if(x == y) break;
        d = gcd( (y - x + n) % n,n);
        if(d >1 && d < n) return d;
        if(i == k) y = x,k <<= 1;
    }while(true);
    return n;
}
LL Find(LL num)
{
    LL res = 0;
    do{
        res = pollard(rand()%(num-1)+1,num);
    }while(res >= num);
    return  res;
}
int main()
{
    //freopen("1.txt","r",stdin);
    //freopen("5.txt","w",stdout);
    LL n;
    while(~scanf("%I64d",&n))
    {
        bool flag = false;
        if(n == 1 || MR(5,n)) flag = true;
        else 
        {
            LL num1 = Find(n);
            LL num2 = n / num1;
            if(num1 > num2) swap(num1,num2);
            if(MR(5,num1) && MR(5,num2) && num1 != num2 || num1 * num1 == num2)
                printf("%I64d %I64d %I64d\n",num1,num2,n);
            else flag = true;
        }
        if(flag) printf("is not a D_num\n");
    }
    return 0;
}        

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值