主页 / 问题 / 1193474
Accepted
A_Hatake
Asked:2021-10-21 23:54:55 +0800 CST2021-10-21 23:54:55 +0800 CST

大数的阶乘

  • 772

您需要计算数字 N 的阶乘，其中 1<=N<=1000。您还需要在响应中提供最多 3000 个字符。我写了一个程序，它在编译器中运行良好：

#include <stdio.h>
#include <gmp.h>

static void
factorial (long n, mpz_t r)
{ 
  mpz_init_set_si (r, 1);
  for (; n > 1; n--) {
    mpz_mul_si (r, r, n);
  }
}

int
main (void)
{
  int n;
  mpz_t r;
  while (scanf ("%d", &n) == 1) {
    factorial (n, r);
    gmp_printf ("%Zd\n", r);
  }
  return 0;
}

但是测试系统不接受代码，因为 不知道 GMP 库。告诉我如何解决这个问题？

c

3 个回答

  1. 最有可能的是，测试系统要求您知道算术的工作时间，以及用代码编写它的能力（即对编程语言有足够的了解）。系统的目的不太可能是检查代码的简洁性或简单性。

    但是，如果您真的不想手动实现长算法，那么我看到了另一种方法：获取库代码并将其复制粘贴到您的文件中。这可能会导致问题：

    • 系统很可能对解决方案的大小有限制，您可以通过这种方式超过它；
    • 如果你的决定也会被人检查（例如，如果你以学生的身份完成任务，有时老师会这样做），那么你的伎俩就会被揭露；

    事实上，如果你的任务是专门用 C 编写代码，这种技巧几乎等同于作弊。如果 Java 或 Python 也可以从选项中获得（它们在标准库级别支持长数字），那么这更像是一个机会的平衡:)

    PS：这是可以内置到您的解决方案中的库的变体：. 测试甚至有一个使用它的阶乘计算示例。

    • 3
  2. Best Answer

    简而言之，这是系统最终接受的代码：

    #include<stdio.h>
 
#define MAX 3000

int multiply(int x, int res[], int res_size);
 
void factorial(int n)
{
    int res[MAX];
     
    res[0] = 1;
    int res_size = 1;
     
    for (int x=2; x<=n; x++)
        res_size = multiply(x, res, res_size);
 
    for (int i=res_size-1; i>=0; i--)
        printf ("%d",res[i]);
}
 
int multiply(int x, int res[], int res_size)
{
    int carry = 0;  // Инициализируем перенос
     
    for (int i=0; i<res_size; i++)
    {
        int prod = res[i] * x + carry;
         
        res[i] = prod % 10;  
         
        carry  = prod/10;    
    }
    
    while (carry)
    {
        res[res_size] = carry%10;
        carry = carry/10;
        res_size++;
    }
    return res_size;
}
 
int main()
{   int n;
    scanf ("%d", &n);
    factorial(n);
    return 0;
}
    • 3

  3. 这是另一种解决方案，但测试系统再次不想接受它......

    int main(int argc, char** argv)
{
    int i,j,d,sz,n=30;
    char *f=calloc(2,sizeof(char));
    *f='1';
    for(d=0,i=1; i<=n; i++)
    {
        sz=strlen(f);
        for(j=sz-1; j>=0; j--,d/=10)
        {
            d+=(f[j]-'0')*i;
            f[j]=d%10+'0';
        }
        for(sz++; d; d/=10,sz++)
        {
            f=realloc(f,(sz+1)*sizeof(char));
            memmove(f+1,f,sz*sizeof(char));
            *f=d%10+'0';
        }
        //printf("%s\n",f);
    }
    printf("%d!=%s\n",n,f);
    free(f);
    system("pause");
    return 0;
}
    • 0

相关问题