結果
問題 |
No.1006 Share an Integer
|
ユーザー |
![]() |
提出日時 | 2020-03-06 23:03:21 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,973 bytes |
コンパイル時間 | 1,696 ms |
コンパイル使用メモリ | 178,280 KB |
実行使用メモリ | 8,884 KB |
最終ジャッジ日時 | 2024-10-14 09:33:23 |
合計ジャッジ時間 | 2,731 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 12 WA * 7 |
ソースコード
#include<bits/stdc++.h> using namespace std; #define int long long const int MAX_N = 2e6+10; const int inf = 1e9; bool arr[MAX_N]; int division[MAX_N]; set<int> se; void eratosthenes(){ for(int i = 0; i < MAX_N; i++){ arr[i] = true;} arr[0] = false; arr[1] = false; for(int i = 2; i < MAX_N; i++){ //素因数列挙 if(arr[i] ){ for(int j = 2; i*j < MAX_N; j++){ arr[i*j] = false; } } } for(int i = 0; i < (int)sqrt(MAX_N); i++){ if( arr[i] ){ se.insert(i); } } return ; } signed main(void){ int N; cin >> N; if( N % 2 == 0 ){ cout << 0 << endl; return 0;} set<pair<int,int>> ans; eratosthenes(); division[1] = 1; int left = max(N/2 - (int)sqrt(N) - 100,1LL); int right = min(N/2 + (int)sqrt(N) + 100,N ); for(int i = left; i <= right; i++){ // iの約数の個数を計算しておく // N/2の√N周りらへんを計算しておけば十分そう(予想) int res = i; int s = 1LL; if( arr[res] ){ division[i] = 2; continue;} // p と 1 for(auto d : se){ if(res % d == 0 ){ int cnt = 1LL; while( res % d == 0 ){ res /= d; cnt++; } s *= cnt; if( res == 1 ){ break;} //for-脱出 }// if res%d }//for-d if( arr[res] ){ s *= 2;} division[i] = s; } int rem = inf; for(int A = left; A <= right; A++){ int B = N-A; if( !(A&&B) ){ continue;} int FA = A-division[A]; int FB = B-division[B]; if( rem > abs(FA-FB) ){ rem = abs(FA-FB); ans.clear(); ans.insert({A,B}); } else if( rem == abs(FA-FB) ){ ans.insert({A,B}); } } for(auto ne : ans ){ cout << ne.first << " " << ne.second << endl;} return 0; }