結果
問題 |
No.375 立方体のN等分 (1)
|
ユーザー |
|
提出日時 | 2020-08-21 01:36:20 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 224 ms / 5,000 ms |
コード長 | 1,718 bytes |
コンパイル時間 | 15,966 ms |
コンパイル使用メモリ | 275,120 KB |
最終ジャッジ日時 | 2025-01-13 04:35:17 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 32 |
ソースコード
#pragma GCC target("avx2") #pragma GCC optimize("03") #pragma GCC optimize("unroll-loops") #include <bits/stdc++.h> using namespace std; typedef long double ld; typedef long long ll; typedef unsigned long long ull; #define endl "\n" #define FOR(i,a,b) for(int i=(a);i<=(b);i++) #define rep(i,n) for(int i=0;i<(n);i++) #define PII pair<int, int> #define PLL pair<ll, ll> #define ALL(x) (x).begin(), (x).end() constexpr int INF=1<<30; constexpr ll LINF=1LL<<60; constexpr ll mod=1e9+7; constexpr int NIL = -1; template<class T>vector<T> vec(int len, T elem) { return vector<T>(len, elem); } // auto dp = vec(52, vec(103, vec(103, INF))); template<class T>inline bool chmax(T &a, const T &b) { if (a<b) { a = b; return 1; } return 0; } template<class T>inline bool chmin(T &a, const T &b) { if (b<a) { a = b; return 1; } return 0; } template<class T>inline int popcount(T a) {return __builtin_popcount(a);} template<class T>inline T emod(T a, T p) { return (a%p + p) % p;} //------------------- vector<ll> fact(ll n) { vector<ll> ans; for(ll i=1; i*i<=n; i++) { if(n%i==0) { ans.push_back(i); if(n!=i*i) ans.push_back(n/i); } } return ans; } int main() { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(15); ll n; cin >> n; ll Tmax = n-1; ll Tmin = n-1; auto factors = fact(n); // for(auto e:factors) cout << e << endl; for(ll e: factors) { auto f3 = fact(n/e); ll mn3 = LINF; for(ll e3: f3) { chmin(mn3, (e3-1)+(n/e/e3-1)); } // cout << mn2 << " " << mn3 << endl; chmin(Tmin, e-1+mn3); } cout << Tmin << " " << Tmax << endl; return 0; }