ABC370-E Avoid K Partition

問題へのリンク

C#のソース

using System;
using System.Collections.Generic;
using System.Linq;

class Program
{
    static string InputPattern = "InputX";

    static List<string> GetInputList()
    {
        var WillReturn = new List<string>();

        if (InputPattern == "Input1") {
            WillReturn.Add("3 3");
            WillReturn.Add("1 2 3");
            //2
        }
        else if (InputPattern == "Input2") {
            WillReturn.Add("5 0");
            WillReturn.Add("0 0 0 0 0");
            //0
        }
        else if (InputPattern == "Input3") {
            WillReturn.Add("10 5");
            WillReturn.Add("-5 -1 -7 6 -6 -2 -5 10 2 -10");
            //428
        }
        else {
            string wkStr;
            while ((wkStr = Console.ReadLine()) != null) WillReturn.Add(wkStr);
        }
        return WillReturn;
    }

    const long Hou = 998244353;

    static void Main()
    {
        List<string> InputList = GetInputList();

        long[] wkArr = InputList[0].Split(' ').Select(pX => long.Parse(pX)).ToArray();
        long K = wkArr[1];
        long[] AArr = InputList[1].Split(' ').Select(pX => long.Parse(pX)).ToArray();

        long Result = Solve(K, AArr);
        Console.WriteLine(Result);
    }

    static long Solve(long pK, long[] pAArr)
    {
        long UB = pAArr.GetUpperBound(0);

        // RunSumの配列を作成
        long[] RunSumArr = pAArr.ToArray();
        for (long I = 1; I <= UB; I++) {
            RunSumArr[I] += RunSumArr[I - 1];
        }

        // 余計に配った数[累積和]なDict
        var MinusDict = new Dictionary<long, long>();

        // 場合の数[区間の開始位置]なDP表
        var InsDualSegmentTree = new DualSegmentTree(UB, Hou);
        InsDualSegmentTree.RangeAdd(0, 0, 1);

        long Answer = 0;
        for (long I = 0; I <= UB; I++) {
            long SendVal = InsDualSegmentTree.GetVal(I);
            if (I > 0) {
                if (MinusDict.ContainsKey(RunSumArr[I - 1])) {
                    SendVal -= MinusDict[RunSumArr[I - 1]];
                    SendVal %= Hou;
                    if (SendVal < 0) SendVal += Hou;
                }
            }

            long PrevRunSum = 0;
            if (I > 0) {
                PrevRunSum = RunSumArr[I - 1];
            }
            long NGRunSum = pK + PrevRunSum;
            if (MinusDict.ContainsKey(NGRunSum) == false) {
                MinusDict[NGRunSum] = 0;
            }
            MinusDict[NGRunSum] += SendVal;
            MinusDict[NGRunSum] %= Hou;

            // 解に計上
            if (RunSumArr[UB] != NGRunSum) {
                Answer += SendVal;
                Answer %= Hou;
            }

            long RangeSta = I + 1;
            long RangeEnd = UB;
            if (RangeSta > UB) continue;
            InsDualSegmentTree.RangeAdd(RangeSta, RangeEnd, SendVal);
        }

        return Answer;
    }
}

// 区間加算、1点取得な双対セグ木(フェニック木使用) (法を指定可能) 
#region DualSegmentTree
internal class DualSegmentTree
{
    private long[] mBitArr; // 内部配列(1オリジンなため、添字0は未使用)
    private long mExternalArrUB;
    private long mHou;

    // ノードのIndexの列挙を返す
    internal IEnumerable<long> GetNodeIndEnum()
    {
        for (long I = 0; I <= GetUB(); I++) {
            yield return I;
        }
    }

    // ノードのUBを返す
    internal long GetUB()
    {
        return mExternalArrUB;
    }

    // コンストラクタ
    // フェニック木の外部配列(0オリジン)のUBを指定
    internal DualSegmentTree(long pExternalArrUB, long pHou)
    {
        mExternalArrUB = pExternalArrUB;

        // フェニック木の外部配列は0オリジンで、
        // フェニック木の内部配列は1オリジンなため、2を足す
        mBitArr = new long[pExternalArrUB + 2];

        mHou = pHou;
    }

    // 双対セグメント木の機能
    // 区間加算
    internal void RangeAdd(long pSta, long pEnd, long AddVal)
    {
        pSta++; // 1オリジンに変更
        pEnd++; // 1オリジンに変更

        long ImosSta = pSta;
        long ImosEnd = pEnd + 1;

        // いもす法
        FenwickTree_Add(ImosSta, AddVal);
        if (ImosEnd <= mBitArr.GetUpperBound(0)) {
            FenwickTree_Add(ImosEnd, -AddVal);
        }
    }

    // 双対セグメント木の機能
    // 1点取得
    internal long GetVal(long pInd)
    {
        pInd++; // 1オリジンに変更
        return FenwickTree_GetSum(1, pInd);
    }

    // フェニック木の機能
    // [pSta,pEnd] のSumを返す
    private long FenwickTree_GetSum(long pSta, long pEnd)
    {
        long Result = FenwickTree_GetSum(pEnd) - FenwickTree_GetSum(pSta - 1);
        Result %= mHou;
        if (Result < 0) Result += mHou;
        return Result;
    }

    // フェニック木の機能
    // [0,pEnd] のSumを返す
    private long FenwickTree_GetSum(long pEnd)
    {
        long Sum = 0;
        while (pEnd >= 1) {
            Sum += mBitArr[pEnd];
            Sum %= mHou;
            pEnd -= pEnd & -pEnd;
        }
        if (Sum < 0) Sum += mHou;
        return Sum;
    }

    // フェニック木の機能
    // [I] に Xを加算
    private void FenwickTree_Add(long pI, long pX)
    {
        pX %= mHou;
        while (pI <= mBitArr.GetUpperBound(0)) {
            mBitArr[pI] += pX;
            mBitArr[pI] %= mHou;
            pI += pI & -pI;
        }
    }
}
#endregion

解説

場合の数[区間の開始位置]なDP表を双対セグ木に持ち、
累積和を使って遷移を高速化してます。

余計に配った数[累積和]をDictで管理し、必要なタイミングで引き算してます。