結果
| 問題 |
No.136 Yet Another GCD Problem
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2018-11-04 17:09:10 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 5 ms / 5,000 ms |
| コード長 | 1,212 bytes |
| コンパイル時間 | 1,313 ms |
| コンパイル使用メモリ | 162,396 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-11-20 19:59:53 |
| 合計ジャッジ時間 | 2,751 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 39 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
// 素数列を作る
// N: 上限
// v: 出力先
void GetPrimes(int N, vector<int> &v) {
v = vector<int>();
int c = 1;
while (c < N) {
c++;
bool flag = true;
int d = (int)sqrt(c);
for (int i = 0; i < v.size(); i++) {
if (c % v[i] == 0) {
flag = false;
break;
}
if (d < v[i])
break;
}
if (flag)
v.push_back(c);
}
}
// 素因数分解を行う
// N: 対象
// p: 素数セット(GetPrimesで取得)
// v: 導かれた指数部
// [return]: 割り切れない場合の残り, 失敗すると-1を返す
int64_t Factorize(int64_t N, const vector<int> &p, vector<int> &v) {
if (p.back() < sqrt(N))
return -1;
for (int i = 0; N != 1 && i < p.size(); i++) {
v.push_back(0);
while (N % p[i] == 0)
{
N /= p[i];
v[i]++;
}
}
return N;
}
int main() {
#ifdef DEBUG
std::ifstream in("/home/share/inputf.in");
std::cin.rdbuf(in.rdbuf());
#endif
int N, K;
cin >> N >> K;
vector<int> p, v;
GetPrimes(N, p);
Factorize(N, p, v);
for(int i = 0; i < v.size(); i++) {
if(v[i] > 0) {
cout << (N / p[i]) << endl;
break;
}
}
return 0;
}