#369930 (D) No.869 ふたつの距離

提出ソース

結果

問題	No.869 ふたつの距離
ユーザー	👑 hos.lyric
提出日時	2019-08-17 00:19:01
言語	D (dmd 2.107.1)
結果	WA
実行時間	-
コード長	5,517 bytes
コンパイル時間	2,551 ms
コンパイル使用メモリ	169,620 KB
実行使用メモリ	8,016 KB
最終ジャッジ日時	2023-09-04 02:30:13
合計ジャッジ時間	14,459 ms
ジャッジサーバーID （参考情報）	judge15 / judge14

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 4,380 KB
testcase_01	AC	2 ms 4,380 KB
testcase_02	AC	1 ms 4,376 KB
testcase_03	AC	1 ms 4,376 KB
testcase_04	WA	-
testcase_05	AC	2 ms 4,380 KB
testcase_06	AC	1 ms 4,380 KB
testcase_07	AC	1 ms 4,376 KB
testcase_08	AC	1 ms 4,376 KB
testcase_09	AC	1 ms 4,376 KB
testcase_10	AC	2 ms 4,380 KB
testcase_11	AC	2 ms 4,380 KB
testcase_12	AC	1 ms 4,376 KB
testcase_13	AC	1 ms 4,376 KB
testcase_14	RE	-
testcase_15	WA	-
testcase_16	WA	-
testcase_17	AC	2 ms 4,380 KB
testcase_18	AC	1 ms 4,376 KB
testcase_19	AC	2 ms 4,376 KB
testcase_20	AC	2 ms 4,380 KB
testcase_21	WA	-
testcase_22	WA	-
testcase_23	AC	2 ms 4,376 KB
testcase_24	AC	2 ms 4,376 KB
testcase_25	AC	1 ms 4,380 KB
testcase_26	AC	1 ms 4,380 KB
testcase_27	WA	-
testcase_28	AC	2 ms 4,376 KB
testcase_29	AC	1 ms 4,376 KB
testcase_30	AC	1 ms 4,380 KB
testcase_31	AC	2 ms 4,380 KB
testcase_32	AC	2 ms 4,376 KB
testcase_33	AC	1 ms 4,376 KB
testcase_34	AC	2 ms 4,376 KB
testcase_35	AC	2 ms 4,380 KB
testcase_36	AC	2 ms 4,380 KB
testcase_37	AC	1 ms 4,380 KB
testcase_38	AC	1 ms 4,376 KB
testcase_39	AC	2 ms 4,376 KB
testcase_40	AC	2 ms 4,376 KB
testcase_41	AC	1 ms 4,380 KB
testcase_42	AC	48 ms 4,376 KB
testcase_43	AC	21 ms 4,380 KB
testcase_44	WA	-
testcase_45	WA	-
testcase_46	AC	5 ms 4,380 KB
testcase_47	WA	-
testcase_48	AC	30 ms 4,376 KB
testcase_49	WA	-
testcase_50	WA	-
testcase_51	WA	-
testcase_52	WA	-
testcase_53	AC	12 ms 4,376 KB
testcase_54	AC	48 ms 4,376 KB
testcase_55	WA	-
testcase_56	AC	11 ms 4,380 KB
testcase_57	WA	-
testcase_58	WA	-
testcase_59	WA	-
testcase_60	AC	46 ms 4,376 KB
testcase_61	WA	-
testcase_62	WA	-
testcase_63	WA	-
testcase_64	WA	-
testcase_65	AC	38 ms 4,376 KB
testcase_66	WA	-
testcase_67	WA	-
testcase_68	WA	-
testcase_69	AC	61 ms 4,380 KB
testcase_70	WA	-
testcase_71	WA	-
testcase_72	AC	63 ms 4,376 KB
testcase_73	WA	-
testcase_74	AC	61 ms 4,384 KB
testcase_75	WA	-
testcase_76	WA	-
testcase_77	WA	-
testcase_78	WA	-
testcase_79	WA	-
testcase_80	WA	-
testcase_81	WA	-
testcase_82	AC	60 ms 4,376 KB
testcase_83	AC	59 ms 4,380 KB
testcase_84	AC	62 ms 4,380 KB
testcase_85	AC	61 ms 4,376 KB
testcase_86	AC	59 ms 4,376 KB
testcase_87	AC	61 ms 4,376 KB
testcase_88	RE	-
testcase_89	RE	-
testcase_90	WA	-
testcase_91	WA	-
testcase_92	AC	64 ms 4,376 KB
testcase_93	AC	1 ms 4,380 KB
testcase_94	AC	2 ms 4,380 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

import std.conv, std.functional, std.range, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.complex, std.container, std.math, std.numeric, std.regex, std.typecons;
import core.bitop;

class EOFException : Throwable { this() { super("EOF"); } }
string[] tokens;
string readToken() { for (; tokens.empty; ) { if (stdin.eof) { throw new EOFException; } tokens = readln.split; } auto token = tokens.front; tokens.popFront; return token; }
int readInt() { return readToken.to!int; }
long readLong() { return readToken.to!long; }
real readReal() { return readToken.to!real; }

bool chmin(T)(ref T t, in T f) { if (t > f) { t = f; return true; } else { return false; } }
bool chmax(T)(ref T t, in T f) { if (t < f) { t = f; return true; } else { return false; } }

int binarySearch(alias pred, T)(in T[] as) { int lo = -1, hi = cast(int)(as.length); for (; lo + 1 < hi; ) { const mid = (lo + hi) >> 1; (unaryFun!pred(as[mid]) ? hi : lo) = mid; } return hi; }
int lowerBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a >= val)); }
int upperBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a > val)); }


int N, A, B;
real[] X, Y;

int calc(real d) {
  int cnt;
  foreach (i; 0 .. N) foreach (j; i + 1 .. N) {
    if ((X[i] - X[j])^^2 + (Y[i] - Y[j])^^2 <= d^^2) {
      ++cnt;
    }
  }
  return cnt;
}

void doIt(int l, int x, int y, int r, int s0, int s1, int dir) {
  debug {
    writefln("N = %d, A = %d, B = %d; B' = %d", N, A, B, N * (N - 1) / 2 - B);
    writefln("  l = %d, x = %d, y = %d, r = %d; s0 = %d, s1 = %d, dir = %d", l, x, y, r, s0, s1, dir);
  }
  
  /*
    |--------20----|
        |----20--------|
      l          x   1   y
  */
  enum real C = 1e-4;
  enum real D = 1e-4 * tan(1.0);
  debug {
    writefln("  C = %.10f", C);
    writefln("  D = %.10f", D);
  }
  
  X = new real[N];
  Y = new real[N];
  foreach (i; 0 .. l) {
    X[i] = C * i;
    Y[i] = 0.0;
  }
  foreach (j; 0 .. x) {
    X[l + j] = 18.0 - D * j;
    Y[l + j] = 0.0;
  }
  X[l + x] = 20.0 + C * (r - 0.5);
  Y[l + x] = 0.0;
  foreach (j; 0 .. y) {
    X[l + x + 1 + j] = 25.0;
    Y[l + x + 1 + j] = 0.0;
  }
  
  debug {
    writeln("  X = ", X);
    writeln("  Y = ", Y);
    writeln("  10: ", calc(10.0));
    writeln("  20: ", calc(20.0));
  }
  assert(calc(10.0) == s0);
  assert(calc(20.0) == B);
  
  int now = calc(10.0);
  foreach (j; 0 .. x) {
    real t0 = 18.0 + C * j;
    real t1 = 8.0 + C * j;
    Tuple!(real, int)[] events;
    events ~= tuple(t0, 0);
    events ~= tuple(t1, 0);
    foreach (i; 0 .. j + x) {
      real t = X[i] + 10.0;
      if (t1 <= t && t <= t0) {
        events ~= tuple(t, +1);
      }
    }
    foreach (i; j + x + 1 .. N) {
      real t = X[i] - 10.0;
      if (t1 <= t && t <= t0) {
        events ~= tuple(t, -1);
      }
    }
    sort(events);
    debug {
      writeln("  events = ", events);
    }
    const eventsLen = cast(int)(events.length);
    int mx = now, cur;
    int hm;
    cur = now;
    foreach_reverse (h; 1 .. eventsLen) {
      cur += events[h][1];
      if (chmax(mx, cur)) {
        hm = h;
      }
    }
    bool ok;
    cur = now;
    foreach_reverse (h; 1 .. eventsLen) {
      cur += events[h][1];
      if (cur == A) {
        ok = true;
        X[l + j] = (events[h][0] + events[h - 1][0]) / 2.0;
        break;
      }
    }
    if (ok) {
      break;
    }
    X[l + j] = (events[hm][0] + events[hm - 1][0]) / 2.0;
    now = mx;
  }
  assert(calc(10.0) == A);
  assert(calc(20.0) == B);
}

void main() {
  debug {
    foreach (n; 1 .. 50) {
      foreach (b; (n * (n - 1) + 3) / 4 + n .. n * (n - 1) / 2 + 1) {
        // writefln("n = %d, b = %d", n, b);
        int mn = n * n, mx = 0;
        foreach (l; 0 .. n + 1) {
          foreach (x; 0 .. (n - l - 1) + 1) {
            const y = n - l - 1 - x;
            foreach (r; 0 .. l + 1) {
              if (n * (n - 1) / 2 - b == 0 * x + r + l * y) {
                const scores = [
                    l * (l - 1) / 2 + (n - l) * (n - l - 1) / 2,
                    (l + x) * (l + x - 1) / 2 + (n - l - x) * (n - l - x - 1) / 2
                ];
                if (scores[0] <= scores[1]) {
                  foreach (score; scores) {
                    chmin(mn, score);
                    chmax(mx, score);
                  }
                }
              }
            }
          }
        }
        // writefln("  [%d, %d]", mn, mx);
        foreach (a; (n * (n - 1) + 3) / 4 .. b - n + 1) {
          assert(mn <= a && a <= mx);
        }
      }
    }
  }
  
  try {
    for (; ; ) {
      N = readInt();
      A = readInt();
      B = readInt();
      
      foreach (l; 0 .. N + 1) {
        foreach (x; 0 .. (N - l - 1) + 1) {
          const y = N - l - 1 - x;
          int r = N * (N - 1) / 2 - B;
          if (y > 0) {
            r %= y;
          }
          if (0 <= r && r <= l && N * (N - 1) / 2 - B == 0 * x + r + l * y) {
            int s0 = l * (l - 1) / 2 + (N - l) * (N - l - 1) / 2;
            int s1 = (l + x) * (l + x - 1) / 2 + (N - l - x) * (N - l - x - 1) / 2;
            if (s0 <= A && A <= s1) {
              doIt(l, x, y, r, s0, s1, 0);
              goto done;
            }
            /*
            if (s1 <= A && A <= s0) {
              doIt(l, x, y, r, s0, s1, 1);
              goto done;
            }
            */
          }
        }
      }
      assert(false);
     done:
      foreach (i; 0 .. N) {
        writefln("%.10f %.10f", X[i], Y[i]);
      }
    }
  } catch (EOFException e) {
  }
}

yukicoder

結果

テストケース

ソースコード