ybt1172 万以内阶乘
时空限制 1000ms/64M
题目描述
求10000以内n的阶乘。
输入
只有一行输入,整数n(0≤n≤100000)。
输出
一行,即n!的值。
样例输入
4
样例输出
24
代码
法一:数组模拟
#include<iostream>
#include<cstring>
#include<iomanip>
using namespace std;
const int N = 10005, MOD = 10000;
int n,a[N];
void fact(int n,int a[]){
memset(a,0,sizeof(int)*N);
a[0] = a[1] = 1;
for (int i=1; i<=n; ++i){
for (int j=1; j<=a[0]; ++j) a[j]*=i;
for (int j=1; j<=a[0]; ++j) //万进制
a[j+1]+=a[j]/MOD,a[j]%=MOD;
if (a[a[0]+1]) a[0]++;
}
}
void output(int a[]){
cout<<a[a[0]]; //单独处理最高位
for (int i=a[0]-1; i>=1; --i) cout<<setw(4)<<setfill('0')<<a[i];
cout<<'\n';
}
int main(){
while (cin>>n){
fact(n,a);
output(a);
}
return 0;
}
法二:重载运算符
#include<iostream>
#include<cstring>
#include<iomanip>
using namespace std;
const int N = 10005, MOD = 10000;
struct bign{
int len,d[N];
bign(){ memset(d,0,sizeof(d)); len=1; } //构造函数,默认初始化
bign(int x){ //构造函数,int初始化
memset(d,0,sizeof(d));
len = 0;
do {
d[++len] = x%MOD;
x /= MOD;
} while (x);
}
bign operator * (const int &b){ //重载* 高*低
bign c; //c=结构体本身*b
c.len = len;
for (int i=1; i<=len; ++i) c.d[i]=d[i]*b;
for (int i=1; i<=len; ++i)
c.d[i+1]+=c.d[i]/MOD,c.d[i]%=MOD;
while (c.d[c.len+1]){
c.len++;
c.d[c.len+1] += c.d[c.len]/MOD;
c.d[c.len] %= MOD;
}
return c;
}
bign operator *= (const int &b){ //重载*=
*this = *this*b;
return *this;
}
friend ostream& operator << (ostream& output,const bign &x){//重载输出流
cout<<x.d[x.len]; //单独处理最高位
for (int i=x.len-1; i>=1; --i) cout<<setw(4)<<setfill('0')<<x.d[i];
return output;
}
};
int main(){
int n;
while (cin>>n){
bign f=1; //int初始化
for (int i=1; i<=n; ++i) f*=i; //需重载高*低、*=
cout<<f<<'\n'; //需重载输出流
}
return 0;
}