結果

問題 No.2443 特殊線形群の標準表現
コンテスト
ユーザー nonon
提出日時 2023-08-25 22:56:03
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,425 bytes
コンパイル時間 4,021 ms
コンパイル使用メモリ 269,692 KB
最終ジャッジ日時 2025-02-16 14:17:14
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 6 WA * 15
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
#include<atcoder/all>
using namespace std;
using namespace atcoder;
using ll = long long;
using db = long double;
using ch = char;
using bl = bool;
using st = string;
using pll = pair<ll,ll>;
using psl = pair<st,ll>;
using vst = vector<st>;
using vch = vector<ch>;
using vvch = vector<vch>;
using vbl = vector<bl>;
using vvbl = vector<vbl>;
using vdb = vector<db>;
using vpll = vector<pll>;
using vvpll = vector<vpll>;
using vpsl = vector<psl>;
using vi = vector<int>;
using vvi = vector<vi>;
using vvvi = vector<vvi>;
using vvvvi = vector<vvvi>;
using vll = vector<ll>;
using vvll = vector<vll>;
using vvvll = vector<vvll>;
using vvvvll = vector<vvvll>;
using vvvvvll = vector<vvvvll>;
template <class T> using pq = priority_queue<T>;
template <class T> using pqg = priority_queue<T,vector<T>,greater<T>>;
#define all(A) A.begin(),A.end()
#define sz(A) (ll)A.size()
#define pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define rep(i, n) for (ll i = 0; i < (ll)(n); i++)
#define rrep(i,a,b) for(ll i=(ll)(a);i<=(ll)(b);i++)
#define drep(i,n) for(ll i=(ll)n-1; i>=0; i--)
#define drrep(i,a,b) for(ll i=(ll)a; i>=(ll)b; i--)

// using mint = modint998244353;
// using mint = modint1000000007;
using mint = modint;
using vm = vector<mint>;
using vvm = vector<vm>; 
using vvvm = vector<vvm>;

// const ll mod = 998244353;
// const ll mod = 1000000007;
const ll INF = 3e18;

template<class T> inline bool chmin(T &a, T b)
{
    if (a>b)
    {
        a=b;
        return true;
    }
    return false;
}

template<class T> inline bool chmax(T &a, T b)
{
    if (a<b)
    {
        a=b;
        return true;
    }
    return false;
}

ll gcd(ll a, ll b)
{
    if (a==0)return b;
    if (b==0)return a;
    ll c=a;
    while(a%b!=0)
    {
        c=a%b;
        a=b;
        b=c;
    }
    return b;
}

ll lcm(ll a,ll b){
    return a/gcd(a,b)*b;
}

ll Pow(ll a, ll n)
{
    ll res=1;
    while(n>0)
    {
        if(n%2==1)res*=a;
        a*=a;
        n/=2;
    }
    return res;
}

ll Pow(ll a, ll n, ll p)
{
    ll res=1;
    while(n>0)
    {
        a%=p;
        if(n%2==1)res*=a;
        a*=a;
        res%=p;
        n/=2;
    }
    return res;
}

void yn(bl ok)
{
    cout << (ok?"Yes\n":"No\n");
    return;
}

ll N,P,Q;
vvm A[100011],A_inv[100011];
vvm S[100011],S_inv[100011];

vvm mtxmul(vvm A, vvm B)
{
    assert(A[0].size()==B.size());
    ll l=A.size(),m=B.size(), n=B[0].size();
    vvm C(l,vm(n,0));
    rep(i,l)rep(k,m)rep(j,n)
    {
        C[i][j]+=A[i][k]*B[k][j];
    }
    return C;
}

vvm mtxpow(vvm A, ll N)
{
    assert(A.size()==A[0].size());
    ll n=A.size();
    vvm B(n,vm(n,0));
    rep(i,n)B[i][i]=1;
    while(N>0)
    {
        if(N&1)B=mtxmul(A,B);
        A=mtxmul(A,A);
        N/=2;
    }
    return B;
}

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cin>>N>>P>>Q;
    mint::set_mod(P);
    vvm E={{mint(1),mint(0)},{mint(0),mint(1)}};
    rrep(i,1,N)
    {
      ll a,b,c,d;
      cin>>a>>b>>c>>d;
      A[i]={{a,b},{c,d}};
      A_inv[i]={{d,-b},{-c,a}};
    }
    S[0]=S_inv[N+1]=S_inv[0]=E;
    rrep(i,1,N)
    {
      S[i]=mtxmul(A[i],S[i-1]);
    }
    drrep(i,N,1)
    {
      S_inv[i]=mtxmul(A_inv[i],S_inv[i+1]);
    }
    while(Q--)
    {
      ll L,R,x,y;
      cin>>L>>R>>x>>y;
      vvm B={{x},{y}};
      B=mtxmul(mtxmul(S[R],S_inv[L]),B);
      cout << B[0][0].val() << ' ' << B[1][0].val() << endl;
    }
    return 0;
}
0