n!的位数-CSDN博客

求一个数阶乘的位数有公式

#include<stdio.h>
#include<memory.h>
#define MAX 100000
long long  h[MAX];
int main()
{
    long long n,num;
    scanf_s("%lld", &n);
    while (n--)
    {
        long long i,j,k,carry;
        long long cnt = 0;
        scanf_s("%lld", &num);
        memset(h, 0, MAX*(sizeof(long long)));
        for(i=1,h[MAX-1]=1;i<=num;i++)
        { 
            for (k = MAX - 1, carry = 0; k >= 0; k--)
            {
                carry += i*h[k];
                h[k] = carry % 100000;
                carry /= 100000;
            }
        }
        for (j = 0; j < MAX&&h[j] == 0; ++j);
        for (; j < MAX; ++j)
        {
            if (h[j] / 10000 > 0)
            {
                cnt += 5;
            }
            else
                while (h[j] > 0)
                {
                    cnt++;
                    h[j] /= 10;
                }
        }
        printf("%lld\n", cnt);
    }
    return 0;
}

原先这样写会超时，输入1000的时候运算时间就很久了，所以需要优化
求一个数n的位数：lg(n!)+1=lgn+lg(n-1)+…+lg1+1;

    for(int i=1;i<=n;i++)//
    {
        sum=sum+log10(i);
    }
    printf("%d\n",(int)sum+1);

可以用斯特林公式，再取对数
原文：https://blog.csdn.net/liangzhaoyang1/article/details/51145807
这里写图片描述

#include<stdio.h>
#include<math.h>
#define PI 3.1415926
#define e 2.71828182845
double a = log10(2.0*PI) / 2;
double b = log10(e);
int main()
{
    long long n,num;
    scanf("%lld", &n);
    while (n--)
    {
        long long cnt = 0;
        scanf("%lld", &num);
        if (num == 1)
        {
            printf("1\n");
        }
        else
        {
            cnt = ceil((a + log10(num) / 2.0 + num*(log10(num) - b)));
            printf("%lld\n", cnt);
        }
    }
    return 0;
}