結果

問題 No.502 階乗を計算するだけ
ユーザー kotamanegi
提出日時 2017-04-07 22:57:58
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
MLE  
実行時間 -
コード長 1,527 bytes
コンパイル時間 932 ms
コンパイル使用メモリ 92,604 KB
実行使用メモリ 814,592 KB
最終ジャッジ日時 2024-07-16 02:51:05
合計ジャッジ時間 4,268 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other MLE * 1 -- * 51
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define  _CRT_SECURE_NO_WARNINGS
#include <stdio.h>
#include <algorithm>
#include <utility>
#include <functional>
#include <cstring>
#include <queue>
#include <stack>
#include <math.h>
#include <iterator>
#include <vector>
#include <string>
#include <set>
#include <math.h>
#include <iostream> 
#include<map>
#include <iomanip>
#include <time.h>
#include <random>
#include <stdlib.h>
#include <list>
#include <typeinfo>
#include <list>
#include <set>
using namespace std;
#define LONG_INF 10000000000000000
#define MAX_MOD 1000000007
#define REP(i,n) for(long long i = 0;i < n;++i)
long long rev(long long n) {
	long long hoge = MAX_MOD-2;
	long long ans = 1;
	while (hoge != 0){
		if (hoge % 2 == 1) {
			ans *= n;
			ans %= MAX_MOD;
		}
		n *= n;
		n %= MAX_MOD;
		hoge /= 2;
	}
	return ans;
}
long long compare(long long a, long long b) {
	long long ans = 1;
	for (long long i = 0;i < b;++i) {
		ans *= (a-i);
		ans %= MAX_MOD;
	}
	for (int i = 1;i <= b;++i) {
		ans *= rev(i);
		ans %= MAX_MOD;
	}
	return ans;
}
long long solve(long long n) {
	long long ans = 1;
	if (n == 1) return 1;
	if (n % 2 == 0) {
		ans *= solve(n / 2);
		ans *= ans;
		ans %= MAX_MOD;
		ans *= compare(n, n / 2);
		ans %= MAX_MOD;
		return ans;
	}
	else {
		ans *= solve(n / 2);
		ans *= ans;
		ans %= MAX_MOD;
		ans *= compare(n, n / 2);
		ans %= MAX_MOD;
		ans *= n - n/ 2;
		ans %= MAX_MOD;
		return ans;
	}
}
int main() {
	long long n;
	cin >> n;
	if (n >= MAX_MOD) {
		cout << 0 << endl;
		return 0;
	}
	cout << solve(n) << endl;
	return 0;
}
0