結果
問題 | No.2718 Best Consonance |
ユーザー |
👑 |
提出日時 | 2024-04-06 06:01:42 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 3,027 bytes |
コンパイル時間 | 4,719 ms |
コンパイル使用メモリ | 264,924 KB |
最終ジャッジ日時 | 2025-02-20 22:26:46 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 14 WA * 21 TLE * 1 |
ソースコード
#include<bits/stdc++.h>#include<atcoder/all>#define rep(i,n) for(int i=0;i<n;i++)using namespace std;using namespace atcoder;typedef long long ll;typedef vector<int> vi;typedef vector<long long> vl;typedef vector<vector<int>> vvi;typedef vector<vector<long long>> vvl;typedef long double ld;typedef pair<int, int> P;ostream& operator<<(ostream& os, const modint& a) {os << a.val(); return os;}template <int m> ostream& operator<<(ostream& os, const static_modint<m>& a) {os << a.val(); return os;}template <int m> ostream& operator<<(ostream& os, const dynamic_modint<m>& a) {os << a.val(); return os;}template<typename T> istream& operator>>(istream& is, vector<T>& v){int n = v.size(); assert(n > 0); rep(i, n) is >> v[i]; return is;}template<typename U, typename T> ostream& operator<<(ostream& os, const pair<U, T>& p){os << p.first << ' ' << p.second; return os;}template<typename T> ostream& operator<<(ostream& os, const vector<T>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : " "); returnos;}template<typename T> ostream& operator<<(ostream& os, const vector<vector<T>>& v){int n = v.size(); rep(i, n) os << v[i] << (i == n - 1 ? "\n" : "");return os;}template<typename T> void chmin(T& a, T b){a = min(a, b);}template<typename T> void chmax(T& a, T b){a = max(a, b);}// Sieve of Eratosthenes// https://youtu.be/UTVg7wzMWQc?t=2774struct Sieve {int n;vector<int> f, primes;Sieve(int n=1):n(n), f(n+1) {f[0] = f[1] = -1;for (ll i = 2; i <= n; ++i) {if (f[i]) continue;primes.push_back(i);f[i] = i;for (ll j = i*i; j <= n; j += i) {if (!f[j]) f[j] = i;}}}bool isPrime(int x) { return f[x] == x;}vector<int> factorList(int x) {vector<int> res;while (x != 1) {res.push_back(f[x]);x /= f[x];}return res;}vector<pair<ll,int>> factor(ll x) {vector<pair<ll,int>> res;for (int p : primes) {int y = 0;while (x%p == 0) x /= p, ++y;if (y != 0) res.emplace_back(p,y);}if (x != 1) res.emplace_back(x,1);return res;}int num_of_divisors(long long x){auto f = factor(x);long long res = 1;for(auto [p, val] : f){res *= val + 1;}return res;}set<long long> divisors(long long x){set<long long> se = {1};auto f = factor(x);for(auto [key, value] : f){set<long long> old = se;for(long long y : old){rep(i, value + 1){se.insert(y);y *= key;}}}return se;}};int main(){int n;cin >> n;long long ans = 0;vector<P> v;rep(i, n){int a, b;cin >> a >> b;v.emplace_back(a, b);}sort(v.begin(), v.end(), [](P& left, P& right){return left.first * left.second > right.first * right.second;});Sieve sieve(450);vector<int> c(200001, 1001001001);rep(i, n){auto f = sieve.divisors(v[i].first);for(int g : f){chmax(ans, (long long)v[i].second * g / c[g]);chmin(c[g], v[i].first);}}cout << ans << "\n";return 0;}