結果
問題 |
No.1143 面積Nの三角形
|
ユーザー |
![]() |
提出日時 | 2020-08-01 00:40:53 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 199 ms / 800 ms |
コード長 | 2,555 bytes |
コンパイル時間 | 1,448 ms |
コンパイル使用メモリ | 134,696 KB |
最終ジャッジ日時 | 2025-01-12 12:03:34 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 2 |
other | AC * 18 |
ソースコード
#include <algorithm> #include <cassert> #include <cctype> #include <climits> #include <cmath> #include <complex> #include <cstdio> #include <cstring> #include <deque> #include <functional> #include <iomanip> #include <iostream> #include <map> #include <numeric> #include <queue> #include <random> #include <set> #include <stack> #include <string> #include <tuple> #include <vector> #include <list> #define rep(i, n) for (int i = 0; i < (int)(n); ++i) //#define cerr if(false) cerr #ifdef DEBUG #define show(...) cerr << #__VA_ARGS__ << " = ", debug(__VA_ARGS__); #else #define show(...) 42 #endif using namespace std; using ll = long long; using pii = pair<int, int>; template <typename T, typename S> ostream& operator<<(ostream& os, pair<T, S> a) { os << '(' << a.first << ',' << a.second << ')'; return os; } template <typename T> ostream& operator<<(ostream& os, vector<T> v) { for (auto x : v) os << x << ' '; return os; } void debug() { cerr << '\n'; } template <typename H, typename... T> void debug(H a, T... b) { cerr << a; if (sizeof...(b)) cerr << ", "; debug(b...); } template <typename T> vector<T> divisor(T n) { vector<T> res; for (T i = 1; i * i <= n; i++) { if (n % i == 0) { res.push_back(i); if (i * i != n) { res.push_back(n / i); } } } sort(res.begin(), res.end()); return res; } template <typename T> map<T, int> prime_factor(T n) { map<T, int> res; for (T i = 2; i * i <= n; i++) { while (n % i == 0) { res[i]++; n /= i; } } if (n != 1) { res[n]++; } return res; } template <typename T> T gcd(T a, T b) { while (b) { a %= b; swap(a, b); } return a; } template <typename T> T lcm(T a, T b) { return a / gcd(a, b) * b; } int main(){ ll n; cin >> n; auto d = divisor(n * n); ll ans = 0; set<vector<ll>>st; for(auto x : d){ for(auto y : d){ if((n * n) % (x * y))continue; //z(z + x + y) = n * n / x / y ll b = x + y; ll c = - n * n / x / y; if(b*b-4*c<0)continue; ll r = (ll) sqrt(b*b-4*c); if(r * r < (b*b-4*c))r++; if(r * r > (b*b-4*c))r--; if(r * r != (b*b-4*c))continue; if(abs(-b+r)%2)continue; auto v = vector<ll>{x,y,(-b+r)/2}; sort(v.begin(), v.end()); st.insert(v); } } cout << st.size() << endl; }