結果

問題 No.461 三角形はいくつ?
ユーザー te-sh
提出日時 2017-12-18 17:14:11
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 179 ms / 5,000 ms
コード長 1,118 bytes
コンパイル時間 820 ms
コンパイル使用メモリ 113,280 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-12 23:11:22
合計ジャッジ時間 4,915 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
/home/linuxbrew/.linuxbrew/opt/dmd/include/dlang/dmd/std/numeric.d(2999): Warning: cannot inline function `std.numeric.gcdImpl!uint.gcdImpl`

ソースコード

diff #

import std.algorithm, std.conv, std.range, std.stdio, std.string;
import std.numeric;   // gcd, fft

void main()
{
  auto n = readln.chomp.to!int;
  auto f = new Frac[][](3);
  foreach (i; 0..n) {
    auto rd = readln.splitter;
    auto p = rd.front.to!int; rd.popFront();
    auto a = rd.front.to!int; rd.popFront();
    auto b = rd.front.to!int;
    auto g = gcd(b, a+b);
    f[p] ~= Frac(b/g, (a+b)/g);
  }

  foreach (i; 0..3) {
    f[i] ~= Frac(0, 1);
    f[i] = f[i].sort().uniq.array;
  }

  auto ans = 0L, u = Frac(1, 1);
  foreach (f0; f[0])
    foreach (f1; f[1]) {
      if (f0+f1 > u) continue;
      ans += f[2].length;
      ans -= f[2].assumeSorted.equalRange(u-f0-f1).length;
      ans -= f[2].assumeSorted.upperBound(u-max(f0, f1)).length;
    }

  writeln(ans);
}

struct Frac
{
  long num, den;

  auto opCmp(Frac rhs)
  {
    return num * rhs.den - rhs.num * den;
  }

  auto opBinary(string op: "+")(Frac rhs)
  {
    return Frac(num * rhs.den + rhs.num * den, den * rhs.den);
  }

  auto opBinary(string op: "-")(Frac rhs)
  {
    return Frac(num * rhs.den - rhs.num * den, den * rhs.den);
  }
}
0