AOJ 758.大数的位数

大数的位数

Time Limit: 1000 ms   Case Time Limit: 1000 ms   Memory Limit: 64 MB
Total Submission: 30   Submission Accepted: 9

 

Description

给出一个整数，确定这个数的阶乘的位数。

Input

多组数据，每行为一个大于等于1且小于等于107的整数n.

Output

对应每个输入数据输出一个结果。

Sample Input

Original

Transformed

10
20

Sample Output

Original

Transformed

7
19

两种方法可能可以实现：

1. 抛去一定的精度来计算
2. 采用数学公式计算

对于第一种方法，我们可以得知，对于两个大数相乘，其较低位对结果的影响较小，因此我们可以舍弃一些数据

采用类似于科学记数法的形式来计算(a*10^b)

这种做法，当n较大时，会由于误差的累积导致错误。

而第二种方法可以采用一些数学公式计算

log₁₀n! = log₁₀(1*2*3*......*n) = log₁₀1+log₁₀2+log₁₀3+...+log₁₀n

显然，等式右侧我们可以轻易得出。

而等式左侧，log₁₀n的意义是n对10的对数，即10^log₁₀n=n

对于一个整数k，若10^k<n<10^k+1，我们可以得出结论，n是k位数

所以，大数n!的位数就是log₁₀n!，由于int的向下取整性质，我们还要将这个数再加上1

虽然题目上说的最大值是107，然而，最大值至少是197(方法1WA)

AC代码:GitHub

/*
By:OhYee
Github:OhYee
HomePage:http://www.oyohyee.com
Email:oyohyee@oyohyee.com
Blog:http://www.cnblogs.com/ohyee/
 
かしこいかわいい？
エリーチカ！
要写出来Хорошо的代码哦~
*/
 
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
using namespace std;
 
//DEBUG MODE
#define debug 0
 
//循环
#define REP(n) for(int o=0;o<n;o++)
 
const int maxn = 10005;
bool tree[maxn];
 
int bit(int n) {
    int a = 1;
    int e = 0;
 
    for(int i = 1;i <= n;i++) {
        a *= i;
        while(a / 100000) {//保证精度
            a /= 10;
            e++;
        }
    }
 
    while(a) {//转化为位数
        a /= 10;
        e++;
    }
 
    return e;
}
 
//n! = 1*2*3*4*5*......*(n-1)*n
//ln n! = ln 1+len2+......ln(n-1)+ln n
int mathbit(int n) {
    double ans = 0;
    for(int i = 1;i <= n;i++)
        ans += log10(i);
    return (int)ans+1;
}
 
 
bool Do() {
    int n;
    if(scanf("%d",&n) == EOF)
        return false;
 
 
    printf("%d
",mathbit(n));
    return true;
}
 
int main() {
    while(Do());
    return 0;
}