結果
問題 | No.741 AscNumber(Easy) |
ユーザー |
![]() |
提出日時 | 2018-10-05 22:46:03 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 42 ms / 2,000 ms |
コード長 | 1,653 bytes |
コンパイル時間 | 1,607 ms |
コンパイル使用メモリ | 196,280 KB |
最終ジャッジ日時 | 2025-01-06 14:19:51 |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 55 |
ソースコード
#include<bits/stdc++.h>#include<random>using namespace std;using ll=long long;using ld=long double;using P=pair<ll,ll>;#define MOD 1000000007LL#define INF 1000000000LL#define EPS 1e-10#define FOR(i,n,m) for(ll i=n;i<(ll)m;i++)#define REP(i,n) FOR(i,0,n)#define DUMP(a) REP(d,a.size()){cout<<a[d];if(d!=a.size()-1)cout<<" ";else cout<<endl;}#define ALL(v) v.begin(),v.end()#define UNIQUE(v) sort(ALL(v));v.erase(unique(ALL(v)),v.end());#define pb push_back#define MAX_N 1000010long long extgcd(long long a, long long b, long long& x, long long& y) {long long d = a;if (b != 0) {d = extgcd(b, a % b, y, x);y -= (a / b) * x;} else {x = 1; y = 0;}return d;}long long mod_inverse(long long a, long long m) {long long x, y;if(extgcd(a, m, x, y) == 1) return (m + x % m) % m;else return -1;}vector<long long> fact(MAX_N+1, INF);long long mod_fact(long long n, long long& e) {if(fact[0] == INF) {fact[0]=1;if(MAX_N != 0) fact[1]=1;for(ll i = 2; i <= MAX_N; ++i) {fact[i] = (fact[i-1] * i) % MOD;}}e = 0;if(n == 0) return 1;long long res = mod_fact(n / MOD, e);e += n / MOD;if((n / MOD) % 2 != 0) return (res * (MOD - fact[n % MOD])) % MOD;return (res * fact[n % MOD]) % MOD;}// return nCklong long mod_comb(long long n, long long k) {if(n < 0 || k < 0 || n < k) return 0;long long e1, e2, e3;long long a1 = mod_fact(n, e1), a2 = mod_fact(k, e2), a3 = mod_fact(n - k, e3);if(e1 > e2 + e3) return 0;return (a1 * mod_inverse((a2 * a3) % MOD, MOD)) % MOD;}int main() {cin.tie(0);ios::sync_with_stdio(false);int n;cin>>n;cout<<mod_comb(n+9,9)<<endl;return 0;}