lightoj 1045 - Digits of Factorial

1045 - Digits of Factorial
Time Limit: 2 second(s)Memory Limit: 32 MB

Factorial of an integer is defined by the following function

f(0) = 1

f(n) = f(n - 1) * n, if(n > 0)

So, factorial of 5 is 120. But in different bases, the factorial may be different. For example, factorial of 5 in base 8 is 170.

In this problem, you have to find the number of digit(s) of the factorial of an integer in a certain base.

Input

Input starts with an integer T (≤ 50000), denoting the number of test cases.

Each case begins with two integers n (0 ≤ n ≤ 106) and base (2 ≤ base ≤ 1000). Both of these integers will be given in decimal.

Output

For each case of input you have to print the case number and the digit(s) of factorial n in the given base.

Sample Input

Output for Sample Input

5

5 10

8 10

22 3

1000000 2

0 100

Case 1: 3

Case 2: 5

Case 3: 45

Case 4: 18488885

Case 5: 1


算N!在K进制下的位数,即计算 [ log(1)+log(2)+...+log(N) ]+1  其中log的底数都是K。在此要先了解计算机是怎么表示对数的。计算机的log默认为自然对数,即以e为底。   或者log10(a), 就是以10为底,其他的都得通过换底公式来表示。

loga(b)=logc(b)/logc(a)




 N!在十进制下的位数就是  log10 ( N ! )  +1  ( 自己找个数验证 ) ,N!在K进制下的位数就是  logK( N ! )  +1  ( K 为底数 ),而计算机不能直接以K为底求对数。所以运用换底公式, logK( N ! )  = log( K ) / log ( N!)   ,注意,等号右边的 log 都是默认以e为底。


正式进入解题了,题目给出N,K,如果每一次输入N,K 都要计算 N!的话,那就有很大的开销,开销为O(N)。

所以可以采取一种办法,先预处理,用double数组 sum[ 1000000 ] 把 log ( N ! )  存起来。数组 应该够的,需要说的是定义大数组尽量放到外面,这样比较妥当,要不运行时会出错。

用sum [ i ] 表示  log(1 )+log(2)+。。。+log(i) , 这样时间开始时花费O(N),之后每次花费O(1),等输入的时候用换底公式处理一下就得到答案了。

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
using namespace std;
double sum[1000010];
void fas()
{
memset(sum,0,sizeof(sum));
for(int i=1;i<=1000000;i++)
{
sum[i]=sum[i-1]+log(i);
}
}
int main()
{
int t,n,base;
scanf("%d",&t);
fas();
int l=1;
while(t--)
{
scanf("%d%d",&n,&base);
printf("Case %d: ",l++);
if(n==0)
printf("1\n");
else
{
printf("%d\n",(int)(sum[n]/log(base)+1));
}

}
return 0;
}


  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值