結果
問題 | No.666 1000000007で割るだけ |
ユーザー | polylogK |
提出日時 | 2018-03-23 22:21:36 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,685 bytes |
コンパイル時間 | 1,655 ms |
コンパイル使用メモリ | 170,280 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-06-24 21:38:04 |
合計ジャッジ時間 | 2,407 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 24 |
ソースコード
#include <bits/stdc++.h> #define pq priority_queue using namespace std; typedef pair<int, int> P; typedef pair<int, P> P2; typedef pair<int, P2> P3; constexpr long long gcd(long long a, long long b){return b ? gcd(b, a % b) : a;} constexpr long long lcm(long long a, long long b){return a / gcd(a, b) * b;} constexpr int INF = 1e9, MOD = INF + 7, around[] = {0, 1, 1, -1, -1, 0, -1, 1, 0, 0}; constexpr int mod_pow(long long x, long long n, const int mod){long long ret=1;while(n){if(n&1)(ret*=x)%=mod;(x*=x)%=mod;n>>=1;}return ret;} template<int n> struct Prime{bool arr[n+1];constexpr bool operator[](int k){return arr[k];}constexpr Prime():arr(){for(int i=2;i<n;i++){arr[i]=true;for(int j=2;j*j<=i;j++){if(!(i%j))arr[i]=false;}}}}; template<int n> struct Factorial{long long arr[n+1],ary[n+1];constexpr Factorial():arr(),ary(){arr[0]=1;ary[0]=1;for(int i=0;i<n;i++){arr[i+1]=arr[i]*(i+1)%MOD;ary[i+1]=mod_pow(arr[i+1],MOD-2,MOD);}}}; constexpr Factorial<1010> fact; constexpr Prime<10> prime; constexpr long long comb(int a, int b){if(a < b) return 0LL;if(!a or !b) return 1LL; long long pos = fact.arr[a], pot = fact.ary[a - b], por = fact.ary[b]; return pos * pot % MOD * por % MOD;} template<int n> struct Bernoulli{long long arr[n+1];constexpr Bernoulli():arr(){arr[0]=1;for(int i=1;i<=n;i++){long long sum=0;for(int j=0;j<i;j++){(sum+=comb(i+1,j)*arr[j]%MOD)%=MOD;}arr[i]=(MOD-mod_pow(i+1,MOD-2,MOD))%MOD*sum%MOD;}}}; constexpr int vx[] = {1, 0, -1, 0}, vy[] = {0, 1, 0, -1}; constexpr long double PI = abs(acos(-1)); constexpr int sqrtN = 512, logN = 32; constexpr long long LINF=1e18; int main(){ long long a, b; cin >> a >> b; cout << a * b % MOD << endl; return 0; }