結果
問題 | No.136 Yet Another GCD Problem |
ユーザー | kaffelun |
提出日時 | 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; }