求n!首位数
使用斯特林公式得到近似结果:
limn→+∞n!2πn−−−√(ne)n=1
之后用10的幂计算取首位即可,注意数字较小时特判,代码如下:
#include <stdio.h>
#include <math.h>
const double PI =3.14159265358979;
const double E =2.718281828;
int main()
{
int n,fn;
double log_n_fac;
while (scanf("%d", &n) != EOF){
log_n_fac = 0.5 * log10(2 * PI *(double)n) + (double)n * log10((double)n / E);
log_n_fac -=(int)log_n_fac;
fn = pow(10, log_n_fac);//Stirling's approximation
switch(n){
case 0:printf("1\n");break;
case 1:printf("1\n");break;
case 2:printf("2\n");break;
case 3:printf("6\n");break;
case 7:printf("5\n");break;
case 8:printf("4\n");break;
default:printf("%d\n", fn);
}
}
return 0;
}
求n^n首位
方法类似,代码:
#include <stdio.h>
#include <math.h>
int main()
{
int n;
scanf("%d",&n);
while(n != 0){
printf("%d\n",(int)pow(10,n*log10(n)-(int)(n*log10(n))));
scanf("%d",&n);
}
return 0;
}
整数质因子分解
从小到大枚举因数,如果这个因数不是素数(如15),那么之前枚举过的素数已经消耗掉了这个因数(如3和5),因此n%15不会为0。所以这题并不需要筛选素数。
#include <stdio.h>
int n;
void solve(){
int i;
int m = n;
for (i = 2; i <= n; i++){
int cnt = 0;
if (m % i) continue;
while (m % i == 0){
m /= i;
cnt++;
}
printf("(%d,%d)", i, cnt);
if (m == 1) break;
}
printf("\n");
}
int main(void)
{
int t,i;
scanf("%d", &t);
for (i = 0; i < t; i++){
scanf("%d", &n);
solve();
}
return 0;
}
n!右端的0的个数
令 n!= 5m2px , m,p,x为非负整数。同时由于n!中2的倍数远多于5的倍数,即p>m,所以m即为答案。要得到m,首先找从1到n中5的倍数,即n/5,而25= 52 ,一个数就含有两个5在内,因此第二轮我们找从1到n/5中5的倍数。由于第一轮中已经被拿掉了一个5,第二轮25就只需再拿出一个5即可,n/25。第三轮则n/125,第i论 n5i ,以此类推:
#include <stdio.h>
int main()
{
int t,i,n,m,z;
scanf("%d", &t);
for (i = 0; i < t; i++){
scanf("%d", &n);
m = 5;z = 0;
while (n >= m){
z += n / m;
m *= 5;
}
printf("case #%d:\n%d\n", i, z);
}
return 0;
}