結果

問題 No.1143 面積Nの三角形
ユーザー emthrm
提出日時 2020-08-01 00:57:20
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 263 ms / 800 ms
コード長 2,065 bytes
コンパイル時間 2,771 ms
コンパイル使用メモリ 204,520 KB
最終ジャッジ日時 2025-01-12 12:05:46
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _USE_MATH_DEFINES
#include <bits/stdc++.h>
using namespace std;
#define FOR(i,m,n) for(int i=(m);i<(n);++i)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()
using ll = long long;
const int INF = 0x3f3f3f3f;
const ll LINF = 0x3f3f3f3f3f3f3f3fLL;
const double EPS = 1e-8;
const int MOD = 1000000007;
// const int MOD = 998244353;
const int dy[] = {1, 0, -1, 0}, dx[] = {0, -1, 0, 1};
const int dy8[] = {1, 1, 0, -1, -1, -1, 0, 1}, dx8[] = {0, -1, -1, -1, 0, 1, 1, 1};
template <typename T, typename U> inline bool chmax(T &a, U b) { return a < b ? (a = b, true) : false; }
template <typename T, typename U> inline bool chmin(T &a, U b) { return a > b ? (a = b, true) : false; }
struct IOSetup {
  IOSetup() {
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);
    cout << fixed << setprecision(20);
  }
} iosetup;

template <typename T>
vector<T> divisor(T val) {
  vector<T> res;
  for (T i = 1; i * i <= val; ++i) {
    if (val % i == 0) {
      res.emplace_back(i);
      if (i * i != val) res.emplace_back(val / i);
    }
  }
  sort(ALL(res));
  return res;
}

vector<double> quadratic_equation(double a, double b, double c) {
  if (abs(a) < EPS) {
    if (abs(b) < EPS) {
      assert(abs(c) >= EPS);
      return {};
    }
    return {-c / b};
  }
  double d = b * b - 4 * a * c;
  if (abs(d) < EPS) return {-0.5 * b / a};
  if (d <= -EPS) return {};
  d = sqrt(d);
  double x1 = -0.5 * (b >= 0 ? b + d : b - d) / a, x2 = c / (x1 * a);
  if (x1 > x2) swap(x1, x2);
  return {x1, x2};
}

int main() {
  ll n; cin >> n;
  vector<ll> d = divisor(n * n);
  int ans = 0;
  REP(i, d.size()) REP(j, d.size()) {
    ll x = d[i], y = d[j];
    if (n * n % (x * y) == 0) {
      for (double e : quadratic_equation(1, x + y, -n * n / (x * y))) {
        ll z = round(e);
        if (z > 0 && (n * n / (x * y)) % z == 0 && n * n / x / y / z == x + y + z) {
          ll s = x + y + z;
          ll a = s - x, b = s - y, c = s - z;
          ans += 0 < a && a <= b && b <= c;
        }
      }
    }
  }
  cout << ans << '\n';
  return 0;
}
0