結果
問題 | No.864 四方演算 |
ユーザー | hitonanode |
提出日時 | 2019-08-16 21:26:06 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 25 ms / 1,000 ms |
コード長 | 4,991 bytes |
コンパイル時間 | 1,686 ms |
コンパイル使用メモリ | 173,884 KB |
実行使用メモリ | 7,700 KB |
最終ジャッジ日時 | 2024-09-25 10:42:02 |
合計ジャッジ時間 | 3,081 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 23 ms
7,520 KB |
testcase_01 | AC | 24 ms
7,576 KB |
testcase_02 | AC | 24 ms
7,612 KB |
testcase_03 | AC | 25 ms
7,520 KB |
testcase_04 | AC | 23 ms
7,532 KB |
testcase_05 | AC | 24 ms
7,588 KB |
testcase_06 | AC | 25 ms
7,680 KB |
testcase_07 | AC | 24 ms
7,424 KB |
testcase_08 | AC | 25 ms
7,680 KB |
testcase_09 | AC | 24 ms
7,468 KB |
testcase_10 | AC | 23 ms
7,580 KB |
testcase_11 | AC | 23 ms
7,532 KB |
testcase_12 | AC | 24 ms
7,528 KB |
testcase_13 | AC | 24 ms
7,592 KB |
testcase_14 | AC | 23 ms
7,576 KB |
testcase_15 | AC | 24 ms
7,692 KB |
testcase_16 | AC | 23 ms
7,700 KB |
testcase_17 | AC | 23 ms
7,628 KB |
testcase_18 | AC | 23 ms
7,572 KB |
testcase_19 | AC | 23 ms
7,520 KB |
testcase_20 | AC | 24 ms
7,512 KB |
testcase_21 | AC | 23 ms
7,576 KB |
testcase_22 | AC | 24 ms
7,680 KB |
testcase_23 | AC | 23 ms
7,648 KB |
testcase_24 | AC | 23 ms
7,512 KB |
testcase_25 | AC | 24 ms
7,648 KB |
testcase_26 | AC | 23 ms
7,676 KB |
testcase_27 | AC | 23 ms
7,584 KB |
testcase_28 | AC | 23 ms
7,604 KB |
testcase_29 | AC | 23 ms
7,672 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; using lint = long long int; using pint = pair<int, int>; using plint = pair<lint, lint>; struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_; #define ALL(x) (x).begin(), (x).end() #define SZ(x) ((lint)(x).size()) #define POW2(n) (1LL << (n)) #define FOR(i, begin, end) for(int i=(begin),i##_end_=(end);i<i##_end_;i++) #define IFOR(i, begin, end) for(int i=(end)-1,i##_begin_=(begin);i>=i##_begin_;i--) #define REP(i, n) FOR(i,0,n) #define IREP(i, n) IFOR(i,0,n) template<typename T> void ndarray(vector<T> &vec, int len) { vec.resize(len); } template<typename T, typename... Args> void ndarray(vector<T> &vec, int len, Args... args) { vec.resize(len); for (auto &v : vec) ndarray(v, args...); } template<typename T> bool mmax(T &m, const T q) { if (m < q) {m = q; return true;} else return false; } template<typename T> bool mmin(T &m, const T q) { if (m > q) {m = q; return true;} else return false; } template<typename T1, typename T2> pair<T1, T2> operator+(const pair<T1, T2> &l, const pair<T1, T2> &r) { return make_pair(l.first + r.first, l.second + r.second); } template<typename T1, typename T2> pair<T1, T2> operator-(const pair<T1, T2> &l, const pair<T1, T2> &r) { return make_pair(l.first - r.first, l.second - r.second); } template<typename T> istream &operator>>(istream &is, vector<T> &vec){ for (auto &v : vec) is >> v; return is; } ///// This part below is only for debug, not used ///// template<typename T> ostream &operator<<(ostream &os, const vector<T> &vec){ os << "["; for (auto v : vec) os << v << ","; os << "]"; return os; } template<typename T> ostream &operator<<(ostream &os, const deque<T> &vec){ os << "deq["; for (auto v : vec) os << v << ","; os << "]"; return os; } template<typename T> ostream &operator<<(ostream &os, const set<T> &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template<typename T> ostream &operator<<(ostream &os, const unordered_set<T> &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template<typename T> ostream &operator<<(ostream &os, const multiset<T> &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template<typename T> ostream &operator<<(ostream &os, const unordered_multiset<T> &vec){ os << "{"; for (auto v : vec) os << v << ","; os << "}"; return os; } template<typename T1, typename T2> ostream &operator<<(ostream &os, const pair<T1, T2> &pa){ os << "(" << pa.first << "," << pa.second << ")"; return os; } template<typename TK, typename TV> ostream &operator<<(ostream &os, const map<TK, TV> &mp){ os << "{"; for (auto v : mp) os << v.first << "=>" << v.second << ","; os << "}"; return os; } template<typename TK, typename TV> ostream &operator<<(ostream &os, const unordered_map<TK, TV> &mp){ os << "{"; for (auto v : mp) os << v.first << "=>" << v.second << ","; os << "}"; return os; } #define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl; ///// END ///// /* #include <ext/pb_ds/assoc_container.hpp> #include <ext/pb_ds/tree_policy.hpp> #include <ext/pb_ds/tag_and_trait.hpp> using namespace __gnu_pbds; // find_by_order(), order_of_key() template<typename TK> using pbds_set = tree<TK, null_type, less<TK>, rb_tree_tag, tree_order_statistics_node_update>; template<typename TK, typename TV> using pbds_map = tree<TK, TV, less<TK>, rb_tree_tag, tree_order_statistics_node_update>; */ vector<int> makePrimeLst(int N) { vector<int> ans; vector<int> alive(N+1, 1); lint now = 2; for (; now * now <= N; now++) { if (alive[now]) ans.push_back(now); for (int t = now; t <= N; t += now) alive[t] = 0; } for (; now <= N; now++) if (alive[now]) ans.push_back(now); return ans; } // 素因数分解 map<lint, int> primeFactorize(lint N, const vector<int> &primeLst) { map<lint, int> ans; for (auto v : primeLst) { while (!(N % v)) { N /= v; ans[v]++; } if (N == 1) return ans; } ans[N]++; return ans; // exit(1); } // 約数 unsorted vector<lint> divisor(lint n, const vector<int> &primeLst) { map<lint, int> factor = primeFactorize(n, primeLst); vector<lint> now{1}, nxt; for (auto pa : factor) { nxt.clear(); for (auto v : now) { REP(i, pa.second + 1) { nxt.push_back(v); v *= pa.first; } } swap(now, nxt); } return now; } lint solve(lint x, lint N) { if (x > N * 2 or x < 2) return 0; if (x <= N + 1) return x - 1; return 2 * N + 1 - x; } int main() { auto primes = makePrimeLst(1000000); lint N, K; cin >> N >> K; auto v = divisor(K, primes); lint ret = 0; for (auto x : v) { lint y = K / x; ret += solve(x, N) * solve(y, N); } cout << ret << endl; }