3

我正在编码以找到两个数字 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;
}
4

1 回答 1

2

您遇到整数溢出 -只有(~2^15*sqrt(2)) 使用 32 位整数i*i才是安全的，例如：i=46340

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

这会导致未定义的行为。

于 2012-09-04T19:53:11.110 回答