找到除数的数量不知道出了什么问题

时间:2012-09-04 18:46:01

标签: c++ primes factorization

我正在编码以找到两个数字n和k的公约数。

我正在使用寻找GCD g的方法来实现它,然后找到GCD的除数。

然而,代码编译但在运行时给出了无响应的消息:(

我已经对此感到震惊..任何人都可以帮忙调试.. 提前致谢

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<vector>

using namespace std;

vector<bool>p(1000002,true);

long int getgcd(long int a, long int b)
{  if(b == 0)
          return a;
   else
          return getgcd(b, a % b);
}

void prime()
{   int i,j;
    for(i=2;i<1000001;i++)
      if(p[i])
        for(j=i*i;j<1000001;j=j+i)
             p[j]=false;
    p[0]=false;
    p[1]=false;
}

void divfind(long int a, long int b)
{
     long int g;
     g = getgcd(a,b);
     int i,s,j=0,ans=1,num=0;
     //short int fo[1000000];
     s=(int) sqrt(g);
     for(i=2;i<s+1;i++)
       if(p[i])
        {
          while(g%i==0)
          {g=g/i;
           num++;
          }
          if(num)
            ans*=++num;
          num=0;
        }
      printf("%d\n",ans);
}

int main()
{
    prime();
    long int n;
    int q;
    scanf("%ld %d",&n,&q);
    while(q--)
    {
        int t;
        long int k;
        scanf("%d%ld",&t,&k);
        if(t==1)
        divfind(n,k);
//        else if(t==2)
//       divi(n,k);
//        else
//        nodivi(n,k);
    }
    return 0;
}

1 个答案:

答案 0 :(得分:2)

你遇到了整数溢出 - i*i只有i=46340(~2 ^ 15 * sqrt(2))才能保证32位整数,例如:

i     i*i
46339 2147302921
46340 2147395600
46341 -2147479015
46342 -2147386332

这会导致未定义的行为。