結果

問題 No.1144 Triangles
ユーザー chielochielo
提出日時 2020-08-01 00:13:59
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 234 ms / 3,000 ms
コード長 1,947 bytes
コンパイル時間 1,602 ms
コンパイル使用メモリ 171,088 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-06 22:34:50
合計ジャッジ時間 5,080 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using ll = long long;
const int maxn = 2e3 + 3;
const ll mo = 1e9 + 7, tp = 166666668;

int n, m;
struct vec {
    ll x, y;
    double a;
    vec operator-(const vec &k) const {
        return {x - k.x, y - k.y};
    }
    ll cross(const vec &k) const {
        return x * k.y - y * k.x;
    }
    ll dot(const vec &k) const {
        return x * k.x + y * k.y;
    }
    void arc() {
        a = atan2(y, x);
    }
} a[maxn], b[maxn];

ll ans = 0, uxs, vxs, uys, vys;

void addU(int i, ll d) {
    uxs += b[i].x * d;
    uys += b[i].y * d;
    uxs %= mo;
    uys %= mo;
}
void addV(int i, ll d) {
    vxs += b[i].x * d;
    vys += b[i].y * d;
    vxs %= mo;
    vys %= mo;
}

bool check(int i, int j) {
    ll c = b[i].cross(b[j]);
    ll d = b[i].dot(b[j]);
    return c > 0 || c == 0 && d > 0;
}

void subsolve() {
    std::sort(b, b + m, [](const vec &u, const vec &v) {
        return u.a < v.a;
    });
    uxs = vxs = uys = vys = 0;
    for(int i = 0; i < m; ++i)
        addV(i, 1);
    for(int i = 0, j = 0; i < m;) {
        while(j < m && j - i < m && check(i, j % m)) {
            addV(j % m, -1);
            addU(j % m, +1);
            j += 1;
        }
        vec us{uxs, uys};
        vec vs{vxs, vys};
        ans = (ans + b[i].cross(us) % mo - b[i].cross(vs) % mo) % mo;
        if(j > i) {
            addU(i, -1);
            i += 1;
        } else {
            addV(i, -1);
            i += 1;
            j += 1;
        }
    }
}

void solve(int s) {
    m = 0;
    for(int j = s + 1; j < n; ++j) {
        if(a[s].x != a[j].x || a[s].y != a[j].y) {
            b[m] = a[j] - a[s];
            b[m].arc();
            m += 1;
        }
    }
    subsolve();
}

int main() {
    scanf("%d", &n);
    for(int i = 0; i < n; ++i) {
        scanf("%lld%lld", &a[i].x, &a[i].y);
    }
    for(int i = 0; i < n; ++i) {
        solve(i);
    }
    printf("%lld\n", (ans % mo + mo) % mo);
    return 0;
}
0