poj1423

公式直接算的，取对数的会超时。

#include<iostream>

#include<cmath>
using namespace std;


int input;
int log(int n){
double sum=0;
int i;
for(int i=2;i<=n;i++){
sum+=log10(1.0*i); 
}
return (int)sum+1;
}
int stirling(int n)
{
const double p=acos(-1.0);
double c1=log10(2*p);
double c2=0.434294481903;
double c3=log10(1.0*n);
int s=1;
if(n>3)
s=(c3+c1)/2+n*(c3-c2)+1;
return s;
}
int main(){
cin>>input;
int size=input;
for(int i=0;i<size;i++){
cin>>input;
cout<<stirling(input)<<endl;
}
return 0;

}

还写了一种算n!的，不过这个是不能用来做这道题的

#include<iostream>
using namespace std;
int bigsize;
unsigned __int64 big[10000000];
unsigned __int64 input;
const unsigned __int64 mask=10000000;
unsigned __int64 value,remain;
void multiple(unsigned __int64 x){
value = big[0]*x;
for(int i=0;i<bigsize;i++){
big[i]=value%mask;
remain=value/mask;
value=big[i+1]*x+remain;
}
if(value!=0){
big[bigsize++]=value;
}


}
int countDigits(){
value=big[0];
int i=0;
int digits=0;
while(i<bigsize){
if(big[i+1]==0){
while(value>0){
digits+=1;
value=value/10;
}
}else{
digits+=7;
}
i++;
value=big[i];
}
return digits;
}
int main(){
cin>>input;
int size=input;
for(int i=0;i<size;i++){
cin>>input;
memset(big,0,sizeof(big));
big[0]=1;
bigsize=1;
for(int i=2;i<=input;i++){
multiple(i);
}
cout<<countDigits()<<endl;
}
return 0;
}