結果

問題 No.502 階乗を計算するだけ
コンテスト
ユーザー Kmcode1
提出日時 2017-04-07 23:30:53
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
RE  
実行時間 -
コード長 1,207 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,083 ms
コンパイル使用メモリ 184,344 KB
実行使用メモリ 19,584 KB
最終ジャッジ日時 2026-03-30 17:58:40
合計ジャッジ時間 13,705 ms
ジャッジサーバーID
（参考情報） 		judge3_1 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 42 RE * 10
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
#include<unordered_set>
#include<unordered_map>
using namespace std;




#define MOD 1000000007

map<pair<long long int,long long int>, long long int> mp;

#define M 1000002
long long int F[M];
long long int ppow(long long int i, long long int j){
	long long int res = 1;
	while (j){
		if(j&1)res *= i;
		if (res >= MOD)res %= MOD;
		i *= i;
		i %= MOD;
		j >>= 1;
	}
	return res;
}
long long int inv[M];

long long int recfact(long long int start, long long int n) {
	long long int i;
	if (mp.count(make_pair(start, n) )){
		return mp[make_pair(start, n)];
	}
	if (n <= M) {
		long long int r = inv[start]*F[start+n-1];
		r %= MOD;
		//for (i = start + 1; i < start + n; i++) r *= i,r%=MOD;
		return r%MOD;
	}
	i = n / 2;
	mp[make_pair(start,n)]=(recfact(start, i) * recfact(start + i, n - i))%MOD;
	return mp[{start, n}];
}
long long int factorial(long  n) { return recfact(1, n); }

int main(){
	F[0] = 1;
	inv[0] = 1;
	for (int i = 1; i < M; i++){
		F[i] = F[i - 1];
		F[i] *= i;
		F[i] %= MOD;
		inv[i] = ppow(F[i], MOD - 1);
	}
	long long int n;
	scanf("%lld", &n);
	if (n >= MOD){
		puts("0");
		return 0;
	}
	long long int c = factorial(n);
	printf("%lld\n", c);
	return 0;
}
0