結果

問題 No.3227 Matrix Query
ユーザー kakel-san
提出日時 2025-08-13 23:22:18
言語 C#
(.NET 8.0.404)
結果
AC  
実行時間 765 ms / 8,000 ms
コード長 5,879 bytes
コンパイル時間 8,444 ms
コンパイル使用メモリ 170,012 KB
実行使用メモリ 234,320 KB
最終ジャッジ日時 2025-08-13 23:22:56
合計ジャッジ時間 26,587 ms
ジャッジサーバーID
(参考情報)
judge1 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 28
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (106 ミリ秒)。
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.0/publish/

ソースコード

diff #

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;

class Program
{
    static int NN => int.Parse(ReadLine());
    static long[] NList => ReadLine().Split().Select(long.Parse).ToArray();
    static long[][] NArr(long n) => Enumerable.Repeat(0, (int)n).Select(_ => NList).ToArray();
    public static void Main()
    {
        Solve();
    }
    static void Solve()
    {
        var c = NList;
        var (k, n) = (c[0], c[1]);
        mod = (int)k;
        var init = new Matrix[n];
        for (var i = 0; i < n; ++i) init[i] = new Matrix(NArr(2));
        var seg = new SegmentTree<Matrix>(init, new SegOp());
        var q = NN;
        var ans = new string[q];
        for (var i = 0; i < q; ++i)
        {
            c = NList;
            var (pos, l, r) = ((int)c[0] - 1, (int)c[1] - 1, (int)c[2]);
            var y = NArr(2);
            seg[pos] = new Matrix(y);
            var ai = seg.Prod(l, r);
            ans[i] = $"{ai.A} {ai.B}\n{ai.C} {ai.D}";
        }
        WriteLine(string.Join("\n", ans));
    }
    struct Matrix
    {
        public long A;
        public long B;
        public long C;
        public long D;
        public Matrix(long a, long b, long c, long d)
        {
            A = a;
            B = b;
            C = c;
            D = d;
        }
        public Matrix(long[][] a)
        {
            A = (a[0][0] % mod + mod) % mod;
            B = (a[0][1] % mod + mod) % mod;
            C = (a[1][0] % mod + mod) % mod;
            D = (a[1][1] % mod + mod) % mod;
        }
    }
    static int mod = 1;
    struct SegOp : ISegmentTreeOperator<Matrix>
    {
        public Matrix Identity => new Matrix(1, 0, 0, 1);
        public Matrix Operate(Matrix x, Matrix y)
        {
            return new Matrix(
                (x.A * y.A + x.B * y.C) % mod,
                (x.A * y.B + x.B * y.D) % mod,
                (x.C * y.A + x.D * y.C) % mod,
                (x.C * y.B + x.D * y.D) % mod);
        }
    }

    interface ISegmentTreeOperator<T>
    {
        T Identity { get; }
        T Operate(T x, T y);
    }
    class SegmentTree<T>
    {
        int _n;
        int size;
        int log;
        T[] d;
        ISegmentTreeOperator<T> op;
        void Update(int k)
        {
            d[k] = op.Operate(d[2 * k], d[2 * k + 1]);
        }
        public SegmentTree(int n, ISegmentTreeOperator<T> op)
        {
            this.op = op;
            _n = n;
            size = 1;
            while (size < n) size <<= 1;
            log = CountRZero(size);
            d = new T[2 * size];
            for (var i = 0; i < d.Length; ++i) d[i] = op.Identity;
        }
        public SegmentTree(T[] v, ISegmentTreeOperator<T> op)
        {
            this.op = op;
            _n = v.Length;
            size = 1;
            while (size < v.Length) size <<= 1;
            log = CountRZero(size);
            d = new T[2 * size];
            for (var i = 0; i < size; ++i) d[i] = op.Identity;
            for (var i = 0; i < v.Length; ++i) d[size + i] = v[i];
            for (var i = size - 1; i >= 1; --i) Update(i);
        }
        int CountRZero(int n)
        {
            var ans = 0;
            while (n % 2 == 0)
            {
                ++ans;
                n >>= 1;
            }
            return ans;
        }
        public T this[int p]
        {
            get { return d[p + size]; }
            set
            {
                p += size;
                d[p] = value;
                for (var i = 1; i <= log; ++i) Update(p >> i);
            }
        }
        public T Prod(int l, int r)
        {
            var sml = op.Identity;
            var smr = op.Identity;
            l += size;
            r += size;
            while (l < r)
            {
                if ((l & 1) != 0) sml = op.Operate(sml, d[l++]);
                if ((r & 1) != 0) smr = op.Operate(d[--r], smr);
                l >>= 1;
                r >>= 1;
            }
            return op.Operate(sml, smr);
        }
        T AllProd() => d[1];
        int MinLeft(int r, Predicate<T> f)
        {
            if (r == 0) return 0;
            r += size;
            T sm = op.Identity;
            do
            {
                r--;
                while (r > 1 && (r % 2) != 0) r >>= 1;
                if (!f(op.Operate(d[r], sm)))
                {
                    while (r < size)
                    {
                        r = 2 * r + 1;
                        if (f(op.Operate(d[r], sm)))
                        {
                            sm = op.Operate(d[r], sm);
                            r--;
                        }
                    }
                    return r + 1 - size;
                }
                sm = op.Operate(d[r], sm);
            }
            while ((r & -r) != r);
            return 0;
        }
        public int MaxRight(int l, Predicate<T> f)
        {
            if (l == _n) return _n;
            l += size;
            T sm = op.Identity;
            do
            {
                while (l % 2 == 0) l >>= 1;
                if (!f(op.Operate(sm, d[l])))
                {
                    while (l < size)
                    {
                        l = 2 * l;
                        if (f(op.Operate(sm, d[l])))
                        {
                            sm = op.Operate(sm, d[l]);
                            ++l;
                        }
                    }
                    return l - size;
                }
                sm = op.Operate(sm, d[l]);
                ++l;
            } while ((l & -l) != l);
            return _n;
        }
    }
}
0