結果
| 問題 |
No.25 有限小数
|
| コンテスト | |
| ユーザー |
codershifth
|
| 提出日時 | 2016-01-27 01:39:29 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 2,608 bytes |
| コンパイル時間 | 1,368 ms |
| コンパイル使用メモリ | 159,896 KB |
| 実行使用メモリ | 6,944 KB |
| 最終ジャッジ日時 | 2024-09-21 17:39:06 |
| 合計ジャッジ時間 | 2,205 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 29 WA * 2 |
ソースコード
#include <bits/stdc++.h>
typedef long long ll;
typedef unsigned long long ull;
#define FOR(i,a,b) for(int (i)=(a);i<(b);i++)
#define REP(i,n) FOR(i,0,n)
#define RANGE(vec) (vec).begin(),(vec).end()
using namespace std;
template<typename T>
T gcd(T a, T b) {
if ( a < b )
std::swap(a,b);
if ( b == 0 )
return a;
return gcd(b, a%b);
}
class FiniteDecimal {
public:
void solve(void) {
ll N,M;
cin>>N>>M;
ll g = gcd(N,M);
// 互いに素にしておく
N /= g;
M /= g;
if (N % M == 0)
{
ll x = N/M;
while (x%10 == 0)
x /= 10;
cout<<x%10<<endl;
return;
}
// 有限少数なら
// K
// N/M = A + ∑ a[k] * 10^(-k) (A は整数) と書ける
// k=1
//
// 式変形すると
// 10^K * N = M * (10*A'+a[K]) = M * B
//
// N,M は互いに素なので M*B が 10 の冪乗であることが必要
// B は 10 の倍数でないので M は 5 or 2 の冪で割りきれる
//
// 逆に M が 5 or 2 で割り切れるなら
//
// N = ∑ c[k] * 10^k と書けて N/M の各項は
//
// c[k]*10^k
// ----------- となる
// 5^k * 2^l
//
// 分子分母に適当に数字w書けて b/10^m の形にできる。
//
int k1 = 0;
ll m = M;
while (m%5 == 0)
{
m /= 5;
++k1;
}
int k2 = 0;
while (m%2 == 0)
{
m /= 2;
++k2;
}
if (m == 1)
{
// (ll)(N*pow(10,K)/M) % 10 だと overflow するので
// 10^K/M のかわりに 2^k1 * 5^k2 をかける
while (k1>0)
{
(N *= 2) %= 10;
--k1;
}
while (k2>0)
{
(N *= 5) %= 10;
--k2;
}
cout<<N<<endl;
}
else
cout<<-1<<endl;
}
};
#if 1
int main(int argc, char *argv[])
{
ios::sync_with_stdio(false);
auto obj = new FiniteDecimal();
obj->solve();
delete obj;
return 0;
}
#endif
codershifth