#964393 (C#) No.2697 Range LIS Query

• 問題: No.2697 Range LIS Query
• ユーザー: tobisatis
• 提出日時: 2024-03-22 23:18:30
• 言語: C# (.NET 8.0.404)
• 結果: AC
• 実行時間: 3,633 ms / 10,000 ms
• コード長: 8,399 bytes
• コンパイル時間: 10,999 ms
• コンパイル使用メモリ: 171,012 KB
• 実行使用メモリ: 287,140 KB
• 最終ジャッジ日時: 2024-09-30 12:36:05
• 合計ジャッジ時間: 46,400 ms
• ジャッジサーバーID （参考情報）: judge2 / judge5

ファイルパターン	結果
sample	AC * 3
other	AC * 15

コンパイルメッセージ

  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (99 ms)。
MSBuild のバージョン 17.9.6+a4ecab324 (.NET)
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.0/publish/

ソースコード

namespace AtCoder;

#nullable enable

using System.Numerics;

static class Extensions
{
    public static T[] Repeat<T>(this int time, Func<T> F) => Enumerable.Range(0, time).Select(_ => F()).ToArray();
}

abstract class SegmentTree<TMonoid>
{
    public interface IOperator
    {
        TMonoid IdentityElement { get; }
        TMonoid Operation(TMonoid x, TMonoid y);
    }

    public interface IBinarySearchOperator { bool F(TMonoid value); }

    public class Body<TOperator> where TOperator : struct, IOperator
    {
        public int BaseSize { get; init; }
        protected int Height { get; init; }
        protected TOperator Op { get; init; }
        protected TMonoid[] Values { get; init; }

        public Body(TMonoid[] initial, TOperator op = default)
        {
            Op = op;
            (BaseSize, Height) = (1, 0);
            while (BaseSize < initial.Length) (BaseSize, Height) = (BaseSize << 1, Height + 1);
            var size = BaseSize << 1;
            Values = new TMonoid[size];
            for (var i = 0; i < size; i++) Values[i] = op.IdentityElement;
            for (var i = 0; i < initial.Length; i++) Values[i + BaseSize] = initial[i];
            for (var i = BaseSize - 1; i > 0; i--) Update(i);
        }

        public TMonoid Product(int l, int r)
        {
            (l, r) = (l + BaseSize, r + BaseSize);
            var (resL, resR) = (Op.IdentityElement, Op.IdentityElement);
            while (l < r)
            {
                if ((l & 1) == 1) resL = Op.Operation(resL, Values[l++]);
                if ((r & 1) == 1) resR = Op.Operation(Values[--r], resR);
                (l, r) = (l >> 1, r >> 1);
            }
            return Op.Operation(resL, resR);
        }

        public void Set(int index, TMonoid value)
        {
            var j = index + BaseSize;
            Values[j] = value;
            for (var i = 1; i <= Height; i++) Update(j >> i);
        }

        // the result might be greater than n
        public int MinIndex<TMinIndexOperator>(TMinIndexOperator op = default) where TMinIndexOperator : struct, IBinarySearchOperator
        {
            if (!op.F(Values[1])) return BaseSize + 1;
            var product = Op.IdentityElement;
            var (res, w, i) = (0, BaseSize, 1);
            while (w > 0)
            {
                var p = Op.Operation(product, Values[i]);
                if (!op.F(p))
                {
                    res += w;
                    i++;
                    product = p;
                }
                w >>= 1;
                i <<= 1;
            }
            return res + 1;
        }

        protected void Update(int k) { Values[k] = Op.Operation(Values[k << 1], Values[(k << 1) + 1]); }
    }
}

abstract class LazySegmentTree<TMonoid, TFunction>
{
    public interface IOperator : SegmentTree<TMonoid>.IOperator
    {
        TFunction IdentityFunction { get; }
        TMonoid Mapping(TFunction f, TMonoid x);
        TFunction Composition(TFunction f, TFunction g);
    }

    public class Body<TOperator> : SegmentTree<TMonoid>.Body<TOperator> where TOperator : struct, IOperator
    {
        TFunction[] Lazy { get; init; }

        public Body(TMonoid[] initial, TOperator op = default) : base(initial, op)
        {
            var size = BaseSize << 1;
            Lazy = new TFunction[size];
            for (var i = 0; i < size; i++) Lazy[i] = op.IdentityFunction;
        }

        public new TMonoid Product(int l, int r)
        {
            if (l == r) return Op.IdentityElement;
            Ascend(l, r);
            return base.Product(l, r);
        }

        public new void Set(int index, TMonoid value)
        {
            var j = index + BaseSize;
            for (var i = Height; i >= 1; i--) Propagate(j >> i);
            Values[j] = value;
            for (var i = 1; i <= Height; i++) Update(j >> i);
        }

        public void Apply(int l, int r, TFunction f)
        {
            if (l == r) return;
            Ascend(l, r);
            (l, r) = (l + BaseSize, r + BaseSize);
            var (l2, r2) = (l, r);
            while (l2 < r2)
            {
                if ((l2 & 1) == 1) Apply(l2++, f);
                if ((r2 & 1) == 1) Apply(--r2, f);
                (l2, r2) = (l2 >> 1, r2 >> 1);
            }
            for (int i = 1; i <= Height; i++)
            {
                if (((l >> i) << i) != l) Update(l >> i);
                if (((r >> i) << i) != r) Update((r - 1) >> i);
            }
        }

        void Ascend(int l, int r)
        {
            (l, r) = (l + BaseSize, r + BaseSize);
            for (int i = Height; i >= 1; i--)
            {
                if (((l >> i) << i) != l) Propagate(l >> i);
                if (((r >> i) << i) != r) Propagate((r - 1) >> i);
            }
        }
        void Propagate(int k)
        {
            Apply(k << 1, Lazy[k]);
            Apply((k << 1) + 1, Lazy[k]);
            Lazy[k] = Op.IdentityFunction;
        }
        void Apply(int k, TFunction f)
        {
            Values[k] = Op.Mapping(f, Values[k]);
            if (k < BaseSize) Lazy[k] = Op.Composition(Lazy[k], f);
        }
    }
}

readonly struct Op : LazySegmentTree<(int[,], int), (int, int)>.IOperator
{
    public (int, int) IdentityFunction => (0, -1);
    public (int[,], int) IdentityElement => (new int[4, 4], 0);
    public (int, int) Composition((int, int) f, (int, int) g)
    {
        if (g.Item2 < 0) return f;
        return g;
    }
    public (int[,], int) Mapping((int, int) f, (int[,], int) x)
    {
        var (fv, fc) = f;
        if (fc < 0) return x;
        var a = new int[4, 4];
        var xc = x.Item2;
        a[fc, fc] = xc;
        return (a, xc);
    }
    public (int[,], int) Operation((int[,], int) x, (int[,], int) y)
    {
        var (xa, xc) = x;
        var (ya, yc) = y;
        var a = new int[4, 4];
        for (var fi = 0; fi < 4; fi++) for (var fj = fi; fj < 4; fj++)
        {
            for (var gi = fj; gi < 4; gi++) for (var gj = gi; gj < 4; gj++)
            {
                a[fi, gj] = Math.Max(a[fi, gj], xa[fi, fj] + ya[gi, gj]);
            }
        }
        return (a, xc + yc);
    }
}

class AtCoder
{
    object? Solve()
    {
        var n = Int();
        var init = new (int[,], int)[n];
        for (var i = 0; i < n; i++)
        {
            var v = Int() - 1;
            var a = new int[4, 4];
            a[v, v] = 1;
            init[i] = (a, 1);
        }
        var tree = new LazySegmentTree<(int[,], int), (int, int)>.Body<Op>(init);
        var ans = new List<int>();
        var q = Int();
        for (var i = 0; i < q; i++)
        {
            var t = Int();
            var l = Int() - 1;
            var r = Int();
            if (t == 1)
            {
                var tns = 0;
                var (a, _) = tree.Product(l, r);
                for (var x = 0; x < 4; x++) for (var y = x; y < 4; y++) tns = Math.Max(tns, a[x, y]);
                ans.Add(tns);
            }
            else
            {
                var x = Int() - 1;
                tree.Apply(l, r, (0, x));
            }
        }
        Out(ans);
        return null;
    }

    public static void Main() => new AtCoder().Run();
    public void Run()
    {
        var res = Solve();
        if (res != null)
        {
            if (res is bool yes) res = yes ? "Yes" : "No";
            sw.WriteLine(res);
        }
        sw.Flush();
    }

    string[] input = Array.Empty<string>();
    int iter = 0;
    readonly StreamWriter sw = new(Console.OpenStandardOutput()) { AutoFlush = false };

    #pragma warning disable IDE0051
    string String()
    {
        while (iter >= input.Length) (input, iter) = (Console.ReadLine()!.Split(' '), 0);
        return input[iter++];
    }
    T Input<T>() where T : IParsable<T> => T.Parse(String(), null);
    int Int() => Input<int>();
    void Out(object? x, string? separator = null)
    {
        separator ??= Environment.NewLine;
        if (x is System.Collections.IEnumerable obj and not string)
        {
            var firstLine = true;
            foreach (var item in obj)
            {
                if (!firstLine) sw.Write(separator);
                firstLine = false;
                sw.Write(item);
            }
        }
        else sw.Write(x);
        sw.WriteLine();
    }
    #pragma warning restore IDE0051
}