結果

問題 No.502 階乗を計算するだけ
ユーザー srup٩(๑`н´๑)۶
提出日時 2017-04-08 00:15:42
言語 C++11(廃止可能性あり)
(gcc 13.3.0)
結果
TLE  
実行時間 -
コード長 1,340 bytes
コンパイル時間 1,132 ms
コンパイル使用メモリ 159,008 KB
実行使用メモリ 13,884 KB
最終ジャッジ日時 2024-07-16 03:30:48
合計ジャッジ時間 6,183 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 33 TLE * 2 -- * 17
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef vector<int> vint;
typedef pair<int,int> pint;
typedef vector<pint> vpint;
#define rep(i,n) for(int i=0;i<(n);i++)
#define REP(i,n) for(int i=n-1;i>=(0);i--)
#define reps(i,f,n) for(int i=(f);i<(n);i++)
#define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)
#define all(v) (v).begin(),(v).end()
#define eall(v) unique(all(v), v.end())
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define chmax(a, b) a = (((a)<(b)) ? (b) : (a))
#define chmin(a, b) a = (((a)>(b)) ? (b) : (a))
const int MOD = 1e9 + 7;
const int INF = 1e9;
const ll INFF = 1e18;

ll N;
//x^k mod
long long powmod(long long x, long long k, long long m){ // 繰り返し二乗法(logk)
	if(k == 0) return 1;
	if(k % 2 == 0) return powmod(x * x % m, k / 2, m);
	else return x * powmod(x, k - 1, m) % m;
}


int main(void){
	cin >> N;
	if(N == 0 || N == 1){
		printf("1\n");
		return 0;
	}

	if(N >= MOD){
		printf("0\n");
	}else{
		ll ans = 1;
		if(N <= MOD / 2){
			reps(i, 1, N + 1) ans *= i, ans %= MOD;
			printf("%lld\n", ans);
		}else{
			reps(i, N + 1, MOD){
				// printf("i\n");
				ans *= i, ans %= MOD;
			}
			// printf("k\n");
			ans = powmod(ans, MOD - 2, MOD);
			ans = -1 * ans + MOD;
			printf("%lld\n", ans);
		}
	}
	return 0;
}
0