結果

問題 No.502 階乗を計算するだけ
ユーザー あんこあんこ
提出日時 2021-04-08 22:47:40
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 3,623 bytes
コンパイル時間 1,556 ms
コンパイル使用メモリ 169,580 KB
実行使用メモリ 4,384 KB
最終ジャッジ日時 2023-09-06 06:23:23
合計ジャッジ時間 9,648 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 RE -
testcase_01 RE -
testcase_02 RE -
testcase_03 RE -
testcase_04 RE -
testcase_05 RE -
testcase_06 RE -
testcase_07 RE -
testcase_08 RE -
testcase_09 RE -
testcase_10 RE -
testcase_11 RE -
testcase_12 RE -
testcase_13 RE -
testcase_14 RE -
testcase_15 RE -
testcase_16 RE -
testcase_17 RE -
testcase_18 RE -
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 RE -
testcase_31 RE -
testcase_32 RE -
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
testcase_39 RE -
testcase_40 RE -
testcase_41 RE -
testcase_42 RE -
testcase_43 RE -
testcase_44 RE -
testcase_45 RE -
testcase_46 RE -
testcase_47 RE -
testcase_48 RE -
testcase_49 RE -
testcase_50 RE -
testcase_51 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll,ll> l_l;
typedef pair<int, int> i_i;

#define fi first
#define se second
// #define rep(i,n) for(int i = 0; i < (n); ++i)
// #define rrep(i,n) for(int i = 1; i <= (n); ++i)
// #define drep(i,n) for(int i = (n)-1; i >= 0; --i)
// #define srep(i,s,t) for (int i = s; i < t; ++i)
#define rng(a) a.begin(),a.end()
#define chmax(x,y) (x = max(x,y))
#define chmin(x,y) (x = min(x,y))
// #define limit(x,l,r) max(l,min(x,r))
// #define lims(x,l,r) (x = max(l,min(x,r)))
// #define isin(x,l,r) ((l) <= (x) && (x) < (r))
// #define pb push_back
// #define sz(x) (int)(x).size()
// #define pcnt __builtin_popcountll
// #define uni(x) x.erase(unique(rng(x)),x.end())
// #define snuke srand((unsigned)clock()+(unsigned)time(NULL));
// #define show(x) cout<<#x<<" = "<<x<<endl;
// #define PQ(T) priority_queue<T,v(T),greater<T> >
// #define bn(x) ((1<<x)-1)
// #define dup(x,y) (((x)+(y)-1)/(y))
// #define newline puts("")
// #define v(T) vector<T>
// #define vv(T) v(v(T))
#define v(T) vector<T>

const double PI = 3.14159265358979323846;
const double EPS = 1e-12;
const ll INF = 0x1fffffffffffffff; // or 0x7fffffffffff
const ll MOD = 1e9+7;
const ll MOD2 = 998244353;

#define rep(i, a, b) for(ll (i) = (a); (i) < (b);++i)

// 逆順になる reverse(rng(ans));
// vector<ll> G[100000];
// cout << fixed << setprecision(15) << y << endl;
// cout << bitset<20>(i) << endl;
// fill(distance[0], distance[row], INF);
// __builtin_popcount()
// stoi(S[j], 0, 2)
// for (ll i = 0;i < (ll)s.size() - k;i++)
// ll dp[桁数][未満フラグ][条件を満たしているか]
// stoi stoll
// bit全探索 n進数全探索
// while (t > 0) cc += t % 9, t /= 9;
// priority_queue<ll> que; que.push(3); que.top(); que.pop();
// rep(i,C)rep(j,C) if(i != j) rep(k,C) if(i != k && j != k)
// 始点終点をデータに含ませる
// memory restriction RE -> CE
// forr 逆順ループsnippet
// 座標 上下の回数さえ分かれば組合せ問題

long long gcd(long long a, long long b) {
	if (b) return gcd(b, a % b);
	return a;
}
long long lcm(long long a, long long b) {
	return a * b / gcd(a, b);
}
long long extgcd(long long a, long long b, long long &x, long long &y) {
	long long g = a; x = 1; y = 0;
	if (b) {
		g = extgcd(b, a % b, y, x);
		y -= a / b * x;
	}
	return g;
}
long long invmod(long long a, long long mod) {
	long long x, y;
	extgcd(a, mod, x, y);
	x %= mod;
	if (x < 0) x += mod;
	return x;
}
long long powmod(long long e, long long x, long long mod) {
	long long prod = 1;
	long long cur = e;
	while(x > 0) {
		if (x & 1) prod = prod * cur % mod;
		cur = cur * cur % mod;
		x <<= 1;
	}
	return prod;
}
vector<int> primes(MOD);
vector<int> eratosthenes(long long n) {
	for (long long i = 2;i < n;++i) primes[i] = i;
	for (long long i = 2;i*i < n;++i) {
		if (primes[i])
			for (long long j = i*i; j < n; j+=i)
				primes[j] = 0;
	}
	auto p = remove(primes.begin(), primes.end(), 0);
	primes.erase(p, primes.end());
	return primes;
}
long long fac(long long n, long long mod) {
	if (n >= mod) return 0;
	if (n) return (n * fac(n-1, mod)) % mod;
	return 1;
}
long long combination(long long n, long long k, long long mod) {
	// 階乗やその逆元をメモ化すると速い
	// でかいやつはnCk/2^Nをパスカルの三角形で作る
	long long res = fac(n, mod);
	res = (res * invmod(fac(k, mod), mod)) % mod;
	res = (res * invmod(fac(n-k, mod), mod)) % mod;
	return res;
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(nullptr);

	ll N;
	cin >> N;
	primes = eratosthenes(N);
	cout << fac(N, MOD) << endl;
	return 0;
}
0