結果

問題 No.1226 I hate Robot Arms
ユーザー chielochielo
提出日時 2020-09-11 23:34:41
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 325 ms / 2,000 ms
コード長 3,707 bytes
コンパイル時間 2,320 ms
コンパイル使用メモリ 206,876 KB
実行使用メモリ 33,476 KB
最終ジャッジ日時 2024-06-10 10:53:09
合計ジャッジ時間 10,208 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 11 ms
32,692 KB
testcase_01 AC 11 ms
32,596 KB
testcase_02 AC 98 ms
32,864 KB
testcase_03 AC 114 ms
32,756 KB
testcase_04 AC 127 ms
33,224 KB
testcase_05 AC 117 ms
33,476 KB
testcase_06 AC 297 ms
33,324 KB
testcase_07 AC 94 ms
32,648 KB
testcase_08 AC 37 ms
33,188 KB
testcase_09 AC 210 ms
32,776 KB
testcase_10 AC 34 ms
32,684 KB
testcase_11 AC 145 ms
33,208 KB
testcase_12 AC 83 ms
32,804 KB
testcase_13 AC 63 ms
33,348 KB
testcase_14 AC 254 ms
32,964 KB
testcase_15 AC 25 ms
33,156 KB
testcase_16 AC 296 ms
33,404 KB
testcase_17 AC 86 ms
32,756 KB
testcase_18 AC 61 ms
32,728 KB
testcase_19 AC 151 ms
33,196 KB
testcase_20 AC 140 ms
33,224 KB
testcase_21 AC 187 ms
33,084 KB
testcase_22 AC 306 ms
33,356 KB
testcase_23 AC 309 ms
33,464 KB
testcase_24 AC 307 ms
33,336 KB
testcase_25 AC 311 ms
33,396 KB
testcase_26 AC 303 ms
33,400 KB
testcase_27 AC 322 ms
33,464 KB
testcase_28 AC 325 ms
33,344 KB
testcase_29 AC 325 ms
33,432 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using db = long double;
const int maxn = 1e5 + 3, lo = 0x2'0000;
const db pi = acos(-1);
using std::cos, std::sin, std::sqrt;

struct SEG {
    struct node {
        db a[2][2], b[2];
        bool m;
        void tag(db u[2][2], db v[2]) {
            db c[2][2];
            for(int i = 0; i < 2; ++i) {
                for(int j = 0; j < 2; ++j) {
                    c[i][j] = 0.0;
                    for(int k = 0; k < 2; ++k) {
                        c[i][j] += a[i][k] * u[k][j];
                    }
                }
            }
            memcpy(a, c, sizeof(a));
            db d[2];
            for(int i = 0; i < 2; ++i) {
                d[i] = 0;
                for(int k = 0; k < 2; ++k)
                    d[i] += b[k] * u[k][i];
                d[i] += v[i];
            }
            memcpy(b, d, sizeof(b));
            m = true;
        }
        void down(node &u, node &v) {
            if(m) {
                u.tag(a, b);
                v.tag(a, b);
                *this = {{{1.0, 0.0}, {0.0, 1.0}}, {0.0, 0.0}, false};
            }
        }
    } seg[lo * 2];
    void init() {
        for(int i = 0; i < lo * 2; ++i) {
            seg[i] = {{{1.0, 0.0}, {0.0, 1.0}}, {0.0, 0.0}, false};
        }
    }
    void modify(int l, int r, db u[2][2], db v[2]) {
        down(l - 1); down(r);
        for(int i = l + lo, j = r + lo; i < j; i /= 2, j /= 2) {
            if(i & 1) {
                seg[i++].tag(u, v);
            }
            if(j & 1) {
                seg[j++].tag(u, v);
            }
        }
    }
    void down(int k) {
        if(k < 0 || k >= lo)
            return;
        for(int i = 1, j = lo; i < lo; j /= 2, i = i * 2 + bool(k & j)) {
            seg[i].down(seg[i * 2], seg[i * 2 + 1]);
        }
    }
    std::pair<db, db> query(int k, db v[2]) {
        down(k);
        db d[2];
        auto &&u = seg[k + lo];
        for(int i = 0; i < 2; ++i) {
            d[i] = 0;
            for(int j = 0; j < 2; ++j)
                d[i] += v[j] * u.a[j][i];
            d[i] += u.b[i];
        }
        return {d[0], d[1]};
    }
} s;

std::pair<db, db> get(int i) {
    db v[2] = {db(i), 0.0};
    return s.query(i, v);
}

int ct[maxn], cl[maxn];

void rotate(int k, int v) {
    db t = db(v) / 360. * 2. * pi;
    auto uu = get(k - 1);
    SEG::node p = {{{1.0, 0.0}, {0.0, 1.0}}, {0.0, 0.0}, false};
    {
        db u[2][2] = {{1.0, 0.0}, {0.0, 1.0}};
        db w[2] = {-uu.first, -uu.second};
        p.tag(u, w);
    }
    {
        db u[2][2] = {{cos(t), sin(t)}, {-sin(t), cos(t)}};
        db w[2] = {0.0, 0.0};
        p.tag(u, w);
    }
    {
        db u[2][2] = {{1.0, 0.0}, {0.0, 1.0}};
        db w[2] = {+uu.first, +uu.second};
        p.tag(u, w);
    }
    s.modify(k, lo, p.a, p.b);
}

void extend(int k, int v) {
    auto uu = get(k - 1);
    auto vv = get(k);
    db dx = vv.first - uu.first;
    db dy = vv.second - uu.second;
    db uvl = sqrt(dx * dx + dy * dy);
    {
        db u[2][2] = {{1.0, 0.0}, {0.0, 1.0}};
        db w[2] = {dx / uvl * v, dy / uvl * v};
        s.modify(k, lo, u, w);
    }
}

int main() {
    s.init();
    int n, q;
    scanf("%d%d", &n, &q);
    for(int i = 1; i <= n; ++i) {
        ct[i] = 0;
        cl[i] = 1;
    }
    for(int i = 0; i < q; ++i) {
        int op, k, v;
        scanf("%d%d", &op, &k);
        if(op == 0) {
            scanf("%d", &v);
            rotate(k, v - ct[k]);
            ct[k] = v;
        } else if(op == 1) {
            scanf("%d", &v);
            extend(k, v - cl[k]);
            cl[k] = v;
        } else {
            auto p = get(k);
            printf("%.9Lf %.9Lf\n", p.first, p.second);
        }
    }
    return 0;
}
0