#138692 (D) No.461 三角形はいくつ？

提出ソース

結果

問題	No.461 三角形はいくつ？
ユーザー	izuru_matsuura
提出日時	2016-12-15 04:25:08
言語	D (dmd 2.106.1)
結果	WA
実行時間	-
コード長	4,244 bytes
コンパイル時間	1,726 ms
コンパイル使用メモリ	154,928 KB
実行使用メモリ	6,948 KB
最終ジャッジ日時	2024-06-12 05:36:43
合計ジャッジ時間	16,522 ms
ジャッジサーバーID （参考情報）	judge3 / judge2

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

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	1 ms 5,248 KB
testcase_01	AC	1 ms 5,376 KB
testcase_02	AC	1 ms 5,376 KB
testcase_03	AC	1 ms 5,376 KB
testcase_04	AC	3 ms 5,376 KB
testcase_05	AC	528 ms 5,376 KB
testcase_06	AC	301 ms 5,376 KB
testcase_07	AC	354 ms 5,376 KB
testcase_08	AC	235 ms 5,376 KB
testcase_09	AC	2 ms 5,376 KB
testcase_10	AC	71 ms 5,376 KB
testcase_11	AC	423 ms 5,376 KB
testcase_12	AC	182 ms 5,376 KB
testcase_13	AC	628 ms 5,376 KB
testcase_14	AC	636 ms 5,376 KB
testcase_15	AC	620 ms 5,376 KB
testcase_16	AC	129 ms 5,376 KB
testcase_17	AC	126 ms 5,376 KB
testcase_18	WA	-
testcase_19	WA	-
testcase_20	AC	587 ms 5,376 KB
testcase_21	AC	622 ms 5,376 KB
testcase_22	AC	616 ms 5,376 KB
testcase_23	AC	181 ms 5,376 KB
testcase_24	AC	174 ms 5,376 KB
testcase_25	AC	6 ms 5,376 KB
testcase_26	AC	7 ms 5,376 KB
testcase_27	AC	1,231 ms 5,376 KB
testcase_28	AC	6 ms 5,376 KB
testcase_29	AC	170 ms 5,376 KB
testcase_30	AC	81 ms 5,376 KB
testcase_31	AC	331 ms 5,376 KB
testcase_32	AC	429 ms 5,376 KB
testcase_33	AC	5 ms 5,376 KB
testcase_34	AC	6 ms 5,376 KB
testcase_35	AC	6 ms 5,376 KB
testcase_36	AC	240 ms 5,376 KB
testcase_37	AC	254 ms 5,376 KB
testcase_38	AC	256 ms 5,376 KB
testcase_39	AC	250 ms 5,376 KB
testcase_40	WA	-
testcase_41	AC	258 ms 5,376 KB
testcase_42	WA	-
testcase_43	WA	-
testcase_44	WA	-

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

ソースコード

raw source code

import std.algorithm;
import std.array;
import std.ascii;
import std.container;
import std.conv;
import std.math;
import std.numeric;
import std.range;
import std.stdio;
import std.string;
import std.typecons;
void log(A...)(A arg) { stderr.writeln(arg); }
int size(T)(in T s) { return cast(int)s.length; }

immutable real EPS = 1e-15;

struct Point {
    real x, y;
    Point opBinary(string op)(in Point p) const if (op == "+" || op == "-") {
        return Point(mixin("x" ~ op ~ "p.x"), mixin("y" ~ op ~ "p.y"));
    }
    Point opBinary(string op)(real k) const if (op == "*" || op == "/") {
        return Point(mixin("x" ~ op ~ "k"), mixin("y" ~ op ~ "k"));
    }
};
real dot(in Point a, in Point b) { return a.x * b.x + a.y * b.y; }
real cross(in Point a, in Point b) { return a.x * b.y - a.y * b.x; }
real norm(in Point a) { return sqrt(dot(a, a)); }
Point rot90(in Point p) { return Point(-p.y, p.x); }
real angle(in Point a) { return atan2(a.y, a.x); }

int ccw(Point a, Point b, Point c){
    b = b - a; c = c - a;                                          
    if (cross(b, c) > EPS) return +1;      // a,b,cの順に反時計周り
    if (cross(b, c) < -EPS) return -1;     // a,b,cの順に時計周り                       
    if (dot(b, c) < 0) return +2;          // c--a--b 直線                       
    if (norm(b) < norm(c)) return -2;      // a--b--c 直線                       
    return 0;                              // a--c--b 直線
}

struct Line {
    Point a, b;
};
bool contains(in Line l, in Point p) { return ccw(l.a, l.b, p) % 2 == 0; }
bool parallel(in Line s, in Line t) { return abs(cross(s.b - s.a, t.b - t.a)) < EPS; }
bool orthogonal(in Line s, in Line t) { return abs(dot(s.b - s.a, t.b - t.a)) < EPS; }
bool equals(in Line s, in Line t) { return parallel(s, t) && contains(s, t.a); }
Point crosspoint(in Line s, in Line t) {
    real d = cross(t.b - t.a, s.b - s.a);
    assert(abs(d) >= EPS);
    return s.a + (s.b - s.a) * cross(t.b - t.a, t.b - s.a) / d;
}

Point projection(in Line l, in Point p) {
    Point u = (p - l.a), v = (l.b - l.a);
    return l.a + (v / norm(v)) * (dot(u, v) / norm(v));
}
Point reflection(in Line l, in Point p) {
    Point h = projection(l, p);
    return p + (h - p) * 2;
}

struct Segment {
    Point a, b;
};
bool intersects(in Segment s, in Segment t) {
    return ccw(s.a, s.b, t.a) * ccw(s.a, s.b, t.b) <= 0 && ccw(t.a, t.b, s.a) * ccw(t.a, t.b, s.b) <= 0;
}
bool contains(in Segment s, in Point p) {
    return ccw(s.a, s.b, p) == 0;
}
real dist(in Segment s, in Point p) {
    Point q = projection(Line(s.a, s.b), p);
    auto t = Segment(p, q);
    if (intersects(s, t)) return norm(t.b - t.a);
    return min(norm(s.a - p), norm(s.b - p));
}
real dist(in Segment s, in Segment t) {
    if (intersects(s, t)) return 0;
    return min( min(dist(s, t.a), dist(s, t.b)),
                min(dist(t, s.a), dist(t, s.b)) );
}

void main() {
    auto N = readln.chomp.to!int;
    real[] Y;
    auto S = new Segment[][2];
    foreach (i; 0 .. N) {
        int p; real a, b; readf("%s %s %s\n",&p, &a, &b);
        if (p == 0) {
            Y ~= b / (a + b);
        } else {
            p--;
            Point s, t;
            if (p == 0) {
                s.x = (a / 2) / (a + b);
                s.y = a / (a + b);
                t.x = a / (a + b);
                t.y = 0;
            } else {
                s.x = (a / 2 + b) / (a + b);
                s.y = a / (a + b);
                t.x = b / (a + b);
                t.y = 0;
            }
            S[p] ~= Segment(s, t);
        }
    }
    Y ~= 0;
    S[0] ~= Segment(Point(0.5, 1), Point(1, 0));
    S[1] ~= Segment(Point(0.5, 1), Point(0, 0));

    sort(Y);
    int count(real sy, real ty) {
        return cast(int)Y.assumeSorted.lowerBound(ty + EPS).upperBound(sy - EPS).length;
    }

    //log(S);
    int ans = 0;
    foreach (s; S[0]) {
        foreach (t; S[1]) {
            if (! s.intersects(t)) continue;
            Point c = crosspoint(cast(Line)s, cast(Line)t);
            //log(c);
            real sy = max(s.b.y, t.b.y);
            real ty = min(s.a.y, t.a.y);
            ans += count(sy, c.y - EPS * 5) + count(c.y + EPS * 5, ty);
        }
    }
    writeln(ans);
}

yukicoder

結果

テストケース

ソースコード