結果

問題 No.1932 動く点 P / Moving Point P
ユーザー hir355hir355
提出日時 2022-05-06 22:19:10
言語 C++23(draft)
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 3,491 ms / 6,000 ms
コード長 6,466 bytes
コンパイル時間 5,963 ms
コンパイル使用メモリ 286,560 KB
実行使用メモリ 73,988 KB
最終ジャッジ日時 2023-09-20 03:16:12
合計ジャッジ時間 45,645 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,376 KB
testcase_01 AC 1,563 ms
69,092 KB
testcase_02 AC 1,215 ms
68,040 KB
testcase_03 AC 1,092 ms
7,740 KB
testcase_04 AC 1,558 ms
20,936 KB
testcase_05 AC 1,407 ms
20,848 KB
testcase_06 AC 623 ms
36,424 KB
testcase_07 AC 3,491 ms
73,932 KB
testcase_08 AC 3,476 ms
73,792 KB
testcase_09 AC 2,576 ms
73,660 KB
testcase_10 AC 2,592 ms
73,656 KB
testcase_11 AC 2,574 ms
73,988 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// #include <atcoder/modint>
// #include <atcoder/fenwicktree>
// #include <atcoder/lazysegtree>
#include <atcoder/segtree>
// #include <atcoder/maxflow>
// #include <atcoder/scc>
// #include <atcoder/mincostflow>
using namespace atcoder;

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#include <bits/stdc++.h>
using namespace std;

constexpr long long INF_LL = 2000000000000000000LL;
constexpr int INF = 2000000000;
constexpr long long MOD = 998244353;
// constexpr long long MOD = 1000000007;

const double PI = acos(-1);

#define all(x) x.begin(), x.end()
#define REP(i, a, b) for(int i = a; i < b; i++)
#define rep(i, n) REP(i, 0, n)

typedef long long ll;
typedef pair<ll, ll> P;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<P> vp;
typedef vector<ll> vl;

int dx[4] = {0, -1, 0, 1};
int dy[4] = {1, 0, -1, 0};
int sign[2] = {1, -1};

template <class T> bool chmax(T &a, T b) {
    if(a < b) {
        a = b;
        return 1;
    }
    return 0;
}
template <class T> bool chmin(T &a, T b) {
    if(a > b) {
        a = b;
        return 1;
    }
    return 0;
}

ll modpow(ll a, ll b, ll m) {
    if(b == 0)
        return 1;
    ll t = modpow(a, b / 2, m);
    if(b & 1) {
        return (t * t % m) * a % m;
    } else {
        return t * t % m;
    }
}

template <class T>
T gcd(T m, T n) {
    if(n == 0) return m;
    return gcd(n, m % n);
}

struct edge {
    int to;
    ll cost;
    edge(int t, ll c) { to = t, cost = c; }
};

typedef vector<vector<edge>> graph;

// using mint = modint998244353;
// // using mint = modint1000000007;
// constexpr int MAX_COM = 1000001;
// mint fac[MAX_COM], ifac[MAX_COM];
// void initfac() {
//     fac[0] = ifac[0] = 1;
//     REP(i, 1, MAX_COM) fac[i] = i * fac[i - 1];
//     REP(i, 1, MAX_COM) ifac[i] = 1 / fac[i];
// }
// mint nCr(int n, int r){
//     if(r < 0 || n < r) return 0;
//     return fac[n] * ifac[n - r] * ifac[r];
// }

template< class T >
struct Matrix {
    vector< vector< T > > A;

    Matrix() {}

    Matrix(size_t n, size_t m) : A(n, vector< T >(m, 0)) {}

    Matrix(size_t n) : A(n, vector< T >(n, 0)) {};

    size_t height() const {
        return (A.size());
    }

    size_t width() const {
        return (A[0].size());
    }

    inline const vector< T > &operator[](int k) const {
        return (A.at(k));
    }

    inline vector< T > &operator[](int k) {
        return (A.at(k));
    }

    static Matrix I(size_t n) {
        Matrix mat(n);
        for(int i = 0; i < n; i++) mat[i][i] = 1;
        return (mat);
    }

    Matrix &operator+=(const Matrix &B) {
        size_t n = height(), m = width();
        assert(n == B.height() && m == B.width());
        for(int i = 0; i < n; i++)
        for(int j = 0; j < m; j++)
            (*this)[i][j] += B[i][j];
        return (*this);
    }

    Matrix &operator-=(const Matrix &B) {
        size_t n = height(), m = width();
        assert(n == B.height() && m == B.width());
        for(int i = 0; i < n; i++)
        for(int j = 0; j < m; j++)
            (*this)[i][j] -= B[i][j];
        return (*this);
    }

    Matrix &operator*=(const Matrix &B) {
        size_t n = height(), m = B.width(), p = width();
        assert(p == B.height());
        vector< vector< T > > C(n, vector< T >(m, 0));
        for(int i = 0; i < n; i++)
        for(int j = 0; j < m; j++)
            for(int k = 0; k < p; k++)
            C[i][j] = (C[i][j] + (*this)[i][k] * B[k][j]);
        A.swap(C);
        return (*this);
    }

    Matrix &operator^=(long long k) {
        Matrix B = Matrix::I(height());
        while(k > 0) {
        if(k & 1) B *= *this;
        *this *= *this;
        k >>= 1LL;
        }
        A.swap(B.A);
        return (*this);
    }

    Matrix operator+(const Matrix &B) const {
        return (Matrix(*this) += B);
    }

    Matrix operator-(const Matrix &B) const {
        return (Matrix(*this) -= B);
    }

    Matrix operator*(const Matrix &B) const {
        return (Matrix(*this) *= B);
    }

    Matrix operator^(const long long k) const {
        return (Matrix(*this) ^= k);
    }

    friend ostream &operator<<(ostream &os, Matrix &p) {
        size_t n = p.height(), m = p.width();
        for(int i = 0; i < n; i++) {
        os << "[";
        for(int j = 0; j < m; j++) {
            os << p[i][j] << (j + 1 == m ? "]\n" : ",");
        }
        }
        return (os);
    }


    T determinant() {
        Matrix B(*this);
        assert(width() == height());
        T ret = 1;
        for(int i = 0; i < width(); i++) {
        int idx = -1;
        for(int j = i; j < width(); j++) {
            if(B[j][i] != 0) idx = j;
        }
        if(idx == -1) return (0);
        if(i != idx) {
            ret *= -1;
            swap(B[i], B[idx]);
        }
        ret *= B[i][i];
        T vv = B[i][i];
        for(int j = 0; j < width(); j++) {
            B[i][j] /= vv;
        }
        for(int j = i + 1; j < width(); j++) {
            T a = B[j][i];
            for(int k = 0; k < width(); k++) {
            B[j][k] -= B[i][k] * a;
            }
        }
        }
        return (ret);
    }
};

typedef Matrix<double> S;
S op(S x, S y){ return y * x; }
S e(){ return Matrix<double>::I(3); }
// typedef int F;
// S mapping(F f, S x){ 
//     return f + x;
// }
// F composition(F f, F g){ 
//     return f + g;
// }
// F id(){ return 0; }

void solve(){
    int n;
    cin >> n;
    segtree<S, op, e> seg(n);
    rep(i, n) {
        double p, q, r;
        cin >> p >> q >> r;
        double rad = r / 180 * PI;
        Matrix<double> a(3), b(3), c(3);
        a = Matrix<double>::I(3);
        c = Matrix<double>::I(3);
        a[0][2] = -p;
        a[1][2] = -q;
        c[0][2] = p;
        c[1][2] = q;
        b[0][0] = cos(rad);
        b[0][1] = -sin(rad);
        b[1][0] = sin(rad);
        b[1][1] = cos(rad);
        b[2][2] = 1;
        auto mat = c * b * a;
        seg.set(i, mat);
    }
    int q;
    cin >> q;
    rep(i, q){
        int s, t;
        double x, y;
        cin >> s >> t >> x >> y;
        auto a = seg.prod(s - 1, t);
        cout << a[0][0] * x + a[0][1] * y + a[0][2] << " " << a[1][0] * x + a[1][1] * y + a[1][2] << endl;
    }
}

int main(){
    cin.tie(nullptr);
    std::ios_base::sync_with_stdio(false);
    std::cout << std::fixed << std::setprecision(16);
    // initfac();
    int t;
    t = 1;
    // cin >> t;
    while (t--)
    {
        solve();
    }
}
0