結果
問題 | No.781 円周上の格子点の数え上げ |
ユーザー |
![]() |
提出日時 | 2019-01-11 22:38:47 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 1,083 ms / 2,000 ms |
コード長 | 1,469 bytes |
コンパイル時間 | 1,790 ms |
コンパイル使用メモリ | 193,992 KB |
最終ジャッジ日時 | 2025-01-06 20:22:52 |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 21 |
ソースコード
#include <bits/stdc++.h>#define rep(i, a, n) for(int i = a; i < n; i++)#define repr(i, a, b) for(int i = a; i >= b; i--)#define int long long#define all(a) a.begin(), a.end()#define chmax(x, y) x = max(x, y)#define chmin(x, y) x = min(x, y)using namespace std;typedef pair<int, int> P;const int mod = 1000000007;const int INF = 1e18;const int MAX_N = 10000010;int prime[MAX_N];bool is_prime[MAX_N + 1];int pf[MAX_N];int sieve(int n){int p = 0;for(int i = 0; i <= n; i++) is_prime[i] = true;is_prime[0] = is_prime[1] = false;for(int i = 2; i <= n; i++){if(is_prime[i]){prime[p++] = i;for(int j = 2*i; j <= n; j += i){is_prime[j] = false;pf[j] = i;}}}return p;}int prime_factor(int n){int res = 1;while(n != 1){int now = pf[n];if(now == 0){if(n % 4 == 1) res *= 2;else if(n != 2) res = 0;break;}int cnt = 0;while(n % now == 0){n /= now;cnt++;}if(now % 4 == 1) res *= cnt + 1;else if(now % 4 == 3 && cnt % 2) return 0;}return res * 4;}signed main(){ios::sync_with_stdio(false);cin.tie(0);int x, y;cin >> x >> y;sieve(y + 1);int ans = 0;for(int i = x; i <= y; i++){ans = max(ans, prime_factor(i));}cout << ans << endl;}