結果
| 問題 |
No.129 お年玉(2)
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2015-01-16 23:03:41 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 3 ms / 5,000 ms |
| コード長 | 1,233 bytes |
| コンパイル時間 | 1,388 ms |
| コンパイル使用メモリ | 161,676 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-11-27 22:52:20 |
| 合計ジャッジ時間 | 2,944 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 46 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef tuple<ll, int, int> Tp;
const ll MOD = 1e9;
ll mod_pow(ll x, ll e){
ll v = 1;
for(;e;e>>=1){
if(e & 1)
v = (v * x) % MOD;
x = (x * x) % MOD;
}
return v;
}
ll inverse(ll a, ll m){
ll b = m, u = 1, v = 0;
while (b) {
ll t = a / b;
swap(a -= t * b, b);
swap(u -= t * v, v);
}
return (u % m + m) % m;
}
Tp calc(int a, int b){
ll rem = 1;
int two = 0, five = 0;
for(int i=a;i<=b;i++){
int j = i;
while(j % 2 == 0){
++two;
j /= 2;
}
while(j % 5 == 0){
++five;
j /= 5;
}
rem = (rem * j) % MOD;
}
return Tp(rem, two, five);
}
int main(){
ll N, M;
cin >> N >> M;
N /= 1000;
N %= M;
auto num = calc(M - N + 1, M);
auto den = calc(1, N);
int two = get<1>(num) - get<1>(den);
int five = get<2>(num) - get<2>(den);
ll res = get<0>(num) * inverse(get<0>(den), MOD) % MOD;
res = (res * mod_pow(2, two)) % MOD;
res = (res * mod_pow(5, five)) % MOD;
cout << res << endl;
return 0;
}