結果
問題 | No.2606 Mirror Relay |
ユーザー | KumaTachiRen |
提出日時 | 2024-01-12 23:37:44 |
言語 | C# (.NET 8.0.203) |
結果 |
TLE
|
実行時間 | - |
コード長 | 18,394 bytes |
コンパイル時間 | 10,351 ms |
コンパイル使用メモリ | 170,852 KB |
実行使用メモリ | 51,744 KB |
最終ジャッジ日時 | 2024-09-28 00:21:16 |
合計ジャッジ時間 | 20,552 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 48 ms
28,928 KB |
testcase_01 | AC | 50 ms
25,472 KB |
testcase_02 | AC | 53 ms
25,600 KB |
testcase_03 | AC | 51 ms
25,592 KB |
testcase_04 | AC | 57 ms
26,624 KB |
testcase_05 | AC | 56 ms
26,748 KB |
testcase_06 | AC | 57 ms
26,752 KB |
testcase_07 | AC | 56 ms
26,744 KB |
testcase_08 | AC | 58 ms
26,752 KB |
testcase_09 | AC | 58 ms
26,880 KB |
testcase_10 | AC | 59 ms
27,008 KB |
testcase_11 | AC | 55 ms
27,008 KB |
testcase_12 | AC | 56 ms
26,880 KB |
testcase_13 | AC | 55 ms
26,880 KB |
testcase_14 | AC | 52 ms
25,728 KB |
testcase_15 | AC | 51 ms
25,596 KB |
testcase_16 | AC | 52 ms
25,600 KB |
testcase_17 | AC | 51 ms
25,856 KB |
testcase_18 | AC | 51 ms
25,596 KB |
testcase_19 | AC | 176 ms
27,772 KB |
testcase_20 | AC | 177 ms
27,760 KB |
testcase_21 | AC | 99 ms
27,904 KB |
testcase_22 | AC | 91 ms
27,624 KB |
testcase_23 | AC | 94 ms
27,648 KB |
testcase_24 | AC | 99 ms
27,776 KB |
testcase_25 | AC | 66 ms
27,648 KB |
testcase_26 | AC | 56 ms
27,136 KB |
testcase_27 | AC | 134 ms
27,864 KB |
testcase_28 | AC | 75 ms
27,776 KB |
testcase_29 | AC | 56 ms
26,880 KB |
testcase_30 | AC | 57 ms
27,008 KB |
testcase_31 | AC | 55 ms
26,748 KB |
testcase_32 | AC | 57 ms
26,880 KB |
testcase_33 | AC | 57 ms
26,752 KB |
testcase_34 | AC | 57 ms
26,880 KB |
testcase_35 | AC | 55 ms
26,624 KB |
testcase_36 | AC | 56 ms
27,136 KB |
testcase_37 | AC | 56 ms
26,880 KB |
testcase_38 | AC | 55 ms
26,752 KB |
testcase_39 | AC | 178 ms
44,024 KB |
testcase_40 | AC | 178 ms
43,904 KB |
testcase_41 | AC | 178 ms
44,028 KB |
testcase_42 | AC | 185 ms
43,900 KB |
testcase_43 | AC | 191 ms
43,900 KB |
testcase_44 | AC | 179 ms
43,776 KB |
testcase_45 | AC | 212 ms
43,772 KB |
testcase_46 | AC | 183 ms
43,776 KB |
testcase_47 | AC | 180 ms
44,032 KB |
testcase_48 | AC | 181 ms
43,648 KB |
testcase_49 | TLE | - |
testcase_50 | -- | - |
testcase_51 | -- | - |
testcase_52 | -- | - |
testcase_53 | -- | - |
testcase_54 | -- | - |
testcase_55 | -- | - |
testcase_56 | -- | - |
testcase_57 | -- | - |
testcase_58 | -- | - |
testcase_59 | -- | - |
testcase_60 | -- | - |
testcase_61 | -- | - |
testcase_62 | -- | - |
testcase_63 | -- | - |
testcase_64 | -- | - |
testcase_65 | -- | - |
testcase_66 | -- | - |
testcase_67 | -- | - |
testcase_68 | -- | - |
testcase_69 | -- | - |
testcase_70 | -- | - |
testcase_71 | -- | - |
testcase_72 | -- | - |
コンパイルメッセージ
復元対象のプロジェクトを決定しています... /home/judge/data/code/main.csproj を復元しました (100 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/
ソースコード
using System; using System.ComponentModel; using System.Collections; using System.Collections.Generic; using System.Diagnostics; using System.Diagnostics.Contracts; using System.IO; using System.Linq; using System.Numerics; using System.Runtime.CompilerServices; using System.Runtime.InteropServices; using System.Text; using static Functions; using static OutputLib; using ModInt = StaticModInt<Mod998244353>; static class Program { static public void Main(string[] args) { Console.SetOut(new StreamWriter(Console.OpenStandardOutput()) { AutoFlush = false }); new Solver(); Console.Out.Flush(); } } public class Solver { public Solver() { Solve(); } void Solve() { var s = rs; int n = s.Length; var rh = new RollingHash(s); var pa = StringLib.Palindrome(s); pa.Sort((p1, p2) => -((p1.r - p1.l) - (p2.r - p2.l))); var set = new HashSet<ulong>(); var lr_list = new List<(int l, int r)>(); for (int i = 0; i < pa.Count; i++) { var (l, r) = pa[i]; while (l < r) { var hash = rh[l, r]; if (set.Contains(hash)) break; set.Add(hash); lr_list.Add((l, r)); l++; r--; } } lr_list.Sort((p1, p2) => (p1.r - p1.l) - (p2.r - p2.l)); var lr = lr_list.ToArray(); int m = lr.Length; var dic = new Dictionary<ulong, int>(); for (int i = 0; i < m; i++) { var (l, r) = lr[i]; dic.Add(rh[l, r], i); } var child = new int[m]; Array.Fill(child, -1); for (int i = 0; i < m; i++) { var (l, r) = lr[i]; l++; r--; if (l < r) child[i] = dic[rh[l, r]]; } var cnt = new int[m]; for (int i = 0; i < pa.Count; i++) { var (l, r) = pa[i]; cnt[dic[rh[l, r]]]++; } for (int i = m - 1; i >= 0; i--) if (child[i] != -1) cnt[child[i]] += cnt[i]; var dp = new long[m]; for (int i = 0; i < m; i++) { var (l, r) = lr[i]; dp[i] = (long)cnt[i] * (r - l); for (r--; l < r; r--) { var hash = rh[l, r]; if (dic.ContainsKey(hash)) { dp[i] += dp[dic[hash]]; break; } } } Write(dp.Max()); } const long INF = 1L << 60; int ri { [MethodImpl(256)] get => (int)sc.Integer(); } (int, int) ri2 { [MethodImpl(256)] get => (ri, ri); } (int, int, int) ri3 { [MethodImpl(256)] get => (ri, ri, ri); } (int, int, int, int) ri4 { [MethodImpl(256)] get => (ri, ri, ri, ri); } (int, int, int, int, int) ri5 { [MethodImpl(256)] get => (ri, ri, ri, ri, ri); } long rl { [MethodImpl(256)] get => sc.Integer(); } (long, long) rl2 { [MethodImpl(256)] get => (rl, rl); } (long, long, long) rl3 { [MethodImpl(256)] get => (rl, rl, rl); } (long, long, long, long) rl4 { [MethodImpl(256)] get => (rl, rl, rl, rl); } (long, long, long, long, long) rl5 { [MethodImpl(256)] get => (rl, rl, rl, rl, rl); } uint rui { [MethodImpl(256)] get => (uint)sc.UInteger(); } ulong rul { [MethodImpl(256)] get => sc.UInteger(); } double rd { [MethodImpl(256)] get => sc.Double(); } string rs { [MethodImpl(256)] get => sc.Scan(); } string rline { [MethodImpl(256)] get => sc.Line(); } public StreamScanner sc = new StreamScanner(Console.OpenStandardInput()); } public static class Functions { [MethodImpl(256)] public static int Popcount(ulong x) { x = (x & 0x5555555555555555UL) + ((x >> 1) & 0x5555555555555555UL); x = (x & 0x3333333333333333UL) + ((x >> 2) & 0x3333333333333333UL); x = (x & 0x0f0f0f0f0f0f0f0fUL) + ((x >> 4) & 0x0f0f0f0f0f0f0f0fUL); x = (x & 0x00ff00ff00ff00ffUL) + ((x >> 8) & 0x00ff00ff00ff00ffUL); x = (x & 0x0000ffff0000ffffUL) + ((x >> 16) & 0x0000ffff0000ffffUL); x = (x & 0x00000000ffffffffUL) + ((x >> 32) & 0x00000000ffffffffUL); return (int)x; } [MethodImpl(256)] public static int Popcount(long x) { x = (x & 0x5555555555555555L) + ((x >> 1) & 0x5555555555555555L); x = (x & 0x3333333333333333L) + ((x >> 2) & 0x3333333333333333L); x = (x & 0x0f0f0f0f0f0f0f0fL) + ((x >> 4) & 0x0f0f0f0f0f0f0f0fL); x = (x & 0x00ff00ff00ff00ffL) + ((x >> 8) & 0x00ff00ff00ff00ffL); x = (x & 0x0000ffff0000ffffL) + ((x >> 16) & 0x0000ffff0000ffffL); x = (x & 0x00000000ffffffffL) + ((x >> 32) & 0x00000000ffffffffL); return (int)x; } [MethodImpl(256)] public static int Popcount(int x) { x = (x & 0x55555555) + ((x >> 1) & 0x55555555); x = (x & 0x33333333) + ((x >> 2) & 0x33333333); x = (x & 0x0f0f0f0f) + ((x >> 4) & 0x0f0f0f0f); x = (x & 0x00ff00ff) + ((x >> 8) & 0x00ff00ff); x = (x & 0x0000ffff) + ((x >> 16) & 0x0000ffff); return x; } [MethodImpl(256)] public static int Ctz(long x) { if (x == 0) return -1; return Popcount((ulong)((x & -x) - 1)); } [MethodImpl(256)] public static int CeilPow2(int n) { int x = 0; while ((1 << x) < n) x++; return x; } [MethodImpl(256)] public static int SafeMod(int x, int m) { int r = x % m; return r < 0 ? r + Math.Abs(m) : r; } [MethodImpl(256)] public static long SafeMod(long x, long m) { long r = x % m; return r < 0 ? r + Math.Abs(m) : r; } [MethodImpl(256)] public static int Sign(long x) => x == 0 ? 0 : (x < 0 ? -1 : 1); [MethodImpl(256)] public static int Sign(double x) => x == 0 ? 0 : (x < 0 ? -1 : 1); [MethodImpl(256)] public static int DigitSum(long n, int d = 10) { long s = 0; while (n > 0) { s += n % d; n /= d; } return (int)s; } [MethodImpl(256)] public static long Floor(long a, long b) => a >= 0 ? a / b : (a + 1) / b - 1; [MethodImpl(256)] public static long Ceil(long a, long b) => a > 0 ? (a - 1) / b + 1 : a / b; [MethodImpl(256)] public static int Gcd(int a, int b) => b == 0 ? a : Gcd(b, a % b); [MethodImpl(256)] public static uint Gcd(uint a, uint b) => b == 0 ? a : Gcd(b, a % b); [MethodImpl(256)] public static long Gcd(long a, long b) => b == 0 ? a : Gcd(b, a % b); [MethodImpl(256)] public static ulong Gcd(ulong a, ulong b) => b == 0 ? a : Gcd(b, a % b); [MethodImpl(256)] public static (long x, long y, long g) ExtGcd(long a, long b) { if (b == 0) return (Sign(a), 0, Math.Abs(a)); long c = SafeMod(a, b); var (x2, y2, g) = ExtGcd(b, c); long x = SafeMod(y2, b); long y = (g - a * x) / b; return (x, y, g); } [MethodImpl(256)] public static void Swap(ref int x, ref int y) { x ^= y; y ^= x; x ^= y; } [MethodImpl(256)] public static void Swap(ref long x, ref long y) { x ^= y; y ^= x; x ^= y; } [MethodImpl(256)] public static void Swap<T>(ref T x, ref T y) { T t = y; y = x; x = t; } [MethodImpl(256)] public static T Clamp<T>(T x, T l, T r) where T : IComparable<T> => x.CompareTo(l) <= 0 ? l : (x.CompareTo(r) <= 0 ? x : r); [MethodImpl(256)] public static T Clamp<T>(ref T x, T l, T r) where T : IComparable<T> => x = x.CompareTo(l) <= 0 ? l : (x.CompareTo(r) <= 0 ? x : r); [MethodImpl(256)] public static void Chmin<T>(ref T x, T y) where T : IComparable<T> { if (x.CompareTo(y) > 0) x = y; } [MethodImpl(256)] public static void Chmax<T>(ref T x, T y) where T : IComparable<T> { if (x.CompareTo(y) < 0) x = y; } [MethodImpl(256)] public static int LowerBound(long[] arr, long val, int l = -1, int r = -1) => LowerBound(arr.AsSpan(), t => Sign(t - val), l, r); [MethodImpl(256)] public static int LowerBound(int[] arr, int val, int l = -1, int r = -1) => LowerBound(arr.AsSpan(), t => t - val, l, r); [MethodImpl(256)] public static int LowerBound<T>(T[] arr, T val, int l = -1, int r = -1) where T : IComparable<T> => LowerBound(arr.AsSpan(), t => t.CompareTo(val), l, r); [MethodImpl(256)] public static int LowerBound<T>(T[] arr, Func<T, int> comp, int l = -1, int r = -1) => LowerBound(arr.AsSpan(), comp, l, r); [MethodImpl(256)] public static int LowerBound<T>(Span<T> data, Func<T, int> comp, int l = -1, int r = -1) { if (data.Length == 0) return -1; if (l == -1) l = 0; if (r == -1) r = data.Length; while (l < r) { int x = (l + r) / 2; if (comp(data[x]) < 0) l = x + 1; else r = x; } return l; } [MethodImpl(256)] public static string ToBase2(long v, int digit = -1) { if (digit == -1) { digit = 0; while ((v >> digit) > 0) digit++; } var c = new string[digit]; for (int i = 0; i < digit; i++) c[digit - 1 - i] = ((v >> i) & 1) == 0 ? "0" : "1"; return string.Join("", c); } [MethodImpl(256)] public static string ToBaseN(long v, int n, int digit = -1) { if (digit == -1) { digit = 0; long pow = 1; while (v >= pow) { digit++; pow *= n; } } var c = new int[digit]; for (int i = 0; i < digit; i++, v /= n) c[digit - 1 - i] = (int)(v % n); return string.Join("", c); } } public class StreamScanner { public StreamScanner(Stream stream) { str = stream; } private readonly Stream str; private readonly byte[] buf = new byte[1024]; private int len, ptr; public bool isEof = false; public bool IsEndOfStream { get { return isEof; } } [MethodImpl(256)] private byte read() { if (isEof) throw new EndOfStreamException(); if (ptr >= len) { ptr = 0; if ((len = str.Read(buf, 0, 1024)) <= 0) { isEof = true; return 0; } } return buf[ptr++]; } [MethodImpl(256)] public char Char() { byte b = 0; do b = read(); while (b < 33 || 126 < b); return (char)b; } [MethodImpl(256)] public string Line() { var sb = new StringBuilder(); for (var b = Char(); b != 10 && !isEof; b = (char)read()) sb.Append(b); return sb.ToString(); } [MethodImpl(256)] public string Scan() { var sb = new StringBuilder(); for (var b = Char(); b >= 33 && b <= 126; b = (char)read()) sb.Append(b); return sb.ToString(); } [MethodImpl(256)] public long Integer() { long ret = 0; byte b = 0; var ng = false; do b = read(); while (b != '-' && (b < '0' || '9' < b)); if (b == '-') { ng = true; b = read(); } for (; true; b = read()) { if (b < '0' || '9' < b) return ng ? -ret : ret; else ret = ret * 10 + b - '0'; } } [MethodImpl(256)] public ulong UInteger() { ulong ret = 0; byte b = 0; do b = read(); while (b < '0' || '9' < b); for (; true; b = read()) { if (b < '0' || '9' < b) return ret; else ret = ret * 10 + b - '0'; } } [MethodImpl(256)] public double Double() => double.Parse(Scan()); } public static class OutputLib { [MethodImpl(256)] public static void WriteJoin<T>(string s, IEnumerable<T> t) => Console.WriteLine(string.Join(s, t)); [MethodImpl(256)] public static void WriteMat<T>(T[,] a) { int sz1 = a.GetLength(0); int sz2 = a.GetLength(1); for (int i = 0; i < sz1; i++) { var b = new T[sz2]; for (int j = 0; j < sz2; j++) b[j] = a[i, j]; WriteJoin(" ", b); } } [MethodImpl(256)] public static void WriteMat<T>(T[][] a) { for (int i = 0; i < a.Length; i++) WriteJoin(" ", a[i]); } [MethodImpl(256)] public static void Write(object t) => Console.WriteLine(t.ToString()); [MethodImpl(256)] public static void Write(string str) => Console.WriteLine(str); [MethodImpl(256)] public static void Write(string str, object arg1) => Console.WriteLine(str, arg1); [MethodImpl(256)] public static void Write(string str, object arg1, object arg2) => Console.WriteLine(str, arg1, arg2); [MethodImpl(256)] public static void Write(string str, object arg1, object arg2, object arg3) => Console.WriteLine(str, arg1, arg2, arg3); [MethodImpl(256)] public static void Write(string str, params object[] arg) => Console.WriteLine(str, arg); [MethodImpl(256)] public static void WriteFlush(object t) { Console.WriteLine(t.ToString()); Console.Out.Flush(); } [MethodImpl(256)] public static void WriteError(object t) => Console.Error.WriteLine(t.ToString()); [MethodImpl(256)] public static void Flush() => Console.Out.Flush(); [MethodImpl(256)] public static void YN(bool t) => Console.WriteLine(t ? "YES" : "NO"); [MethodImpl(256)] public static void Yn(bool t) => Console.WriteLine(t ? "Yes" : "No"); [MethodImpl(256)] public static void yn(bool t) => Console.WriteLine(t ? "yes" : "no"); } public interface IStaticMod { uint Mod { get; } bool IsPrime { get; } } public readonly struct Mod1000000007 : IStaticMod { public uint Mod => 1000000007; public bool IsPrime => true; } public readonly struct Mod998244353 : IStaticMod { public uint Mod => 998244353; public bool IsPrime => true; } public readonly struct StaticModInt<T> : IEquatable<StaticModInt<T>> where T : struct, IStaticMod { internal readonly uint _v; private static readonly T op = default; public int Value => (int)_v; public static int Mod => (int)op.Mod; public static StaticModInt<T> Zero => default; public static StaticModInt<T> One => new StaticModInt<T>(1u); [MethodImpl(256)] public static StaticModInt<T> Raw(int v) { var u = unchecked((uint)v); return new StaticModInt<T>(u); } [MethodImpl(256)] public StaticModInt(long v) : this(Round(v)) { } [MethodImpl(256)] public StaticModInt(ulong v) : this((uint)(v % op.Mod)) { } [MethodImpl(256)] private StaticModInt(uint v) => _v = v; [MethodImpl(256)] private static uint Round(long v) { var x = v % op.Mod; if (x < 0) x += op.Mod; return (uint)x; } [MethodImpl(256)] public static StaticModInt<T> operator ++(StaticModInt<T> v) { var x = v._v + 1; if (x == op.Mod) x = 0; return new StaticModInt<T>(x); } [MethodImpl(256)] public static StaticModInt<T> operator --(StaticModInt<T> v) { var x = v._v; if (x == 0) x = op.Mod; return new StaticModInt<T>(x - 1); } [MethodImpl(256)] public static StaticModInt<T> operator +(StaticModInt<T> lhs, StaticModInt<T> rhs) { var v = lhs._v + rhs._v; if (v >= op.Mod) v -= op.Mod; return new StaticModInt<T>(v); } [MethodImpl(256)] public static StaticModInt<T> operator -(StaticModInt<T> lhs, StaticModInt<T> rhs) { unchecked { var v = lhs._v - rhs._v; if (v >= op.Mod) v += op.Mod; return new StaticModInt<T>(v); } } [MethodImpl(256)] public static StaticModInt<T> operator *(StaticModInt<T> lhs, StaticModInt<T> rhs) => new StaticModInt<T>((uint)((ulong)lhs._v * rhs._v % op.Mod)); [MethodImpl(256)] public static StaticModInt<T> operator /(StaticModInt<T> lhs, StaticModInt<T> rhs) => new StaticModInt<T>((uint)((ulong)lhs._v * Inv(rhs._v) % op.Mod)); [MethodImpl(256)] public static StaticModInt<T> operator +(StaticModInt<T> v) => v; [MethodImpl(256)] public static StaticModInt<T> operator -(StaticModInt<T> v) => new StaticModInt<T>(v._v == 0 ? 0 : op.Mod - v._v); [MethodImpl(256)] public static bool operator ==(StaticModInt<T> lhs, StaticModInt<T> rhs) => lhs._v == rhs._v; [MethodImpl(256)] public static bool operator !=(StaticModInt<T> lhs, StaticModInt<T> rhs) => lhs._v != rhs._v; [MethodImpl(256)] public static implicit operator StaticModInt<T>(int v) => new StaticModInt<T>(v); [MethodImpl(256)] public static implicit operator StaticModInt<T>(uint v) => new StaticModInt<T>((long)v); [MethodImpl(256)] public static implicit operator StaticModInt<T>(long v) => new StaticModInt<T>(v); [MethodImpl(256)] public static implicit operator StaticModInt<T>(ulong v) => new StaticModInt<T>(v); [MethodImpl(256)] public static implicit operator long(StaticModInt<T> v) => v._v; [MethodImpl(256)] public static implicit operator ulong(StaticModInt<T> v) => v._v; [MethodImpl(256)] public StaticModInt<T> Pow(long n) { var x = this; var r = new StaticModInt<T>(1U); while (n > 0) { if ((n & 1) > 0) r *= x; x *= x; n >>= 1; } return r; } [MethodImpl(256)] public StaticModInt<T> Inv() => new StaticModInt<T>(Inv(_v)); [MethodImpl(256)] static ulong Inv(ulong x) { long u = op.Mod, xu = 1, yu = 0, v = (long)x, xv = 0, yv = 1; while (v != 0) { long w = SafeMod(u, v); long q = (u - w) / v; long xw = xu - xv * q; long yw = yu - yv * q; u = v; xu = xv; yu = yv; v = w; xv = xw; yv = yw; } return (ulong)(yu < 0 ? yu + op.Mod : yu); } [MethodImpl(256)] static long SafeMod(long x, long m) { long r = x % m; if (r < 0) r += m; return r; } [MethodImpl(256)] public override string ToString() => _v.ToString(); [MethodImpl(256)] public string ToString(string format, IFormatProvider formatProvider) => _v.ToString(format, formatProvider); [MethodImpl(256)] public override bool Equals(object obj) => obj is StaticModInt<T> m && Equals(m); [MethodImpl(256)] public bool Equals(StaticModInt<T> other) => _v == other._v; [MethodImpl(256)] public override int GetHashCode() => _v.GetHashCode(); } public partial class StringLib { public static List<(int l, int r)> Palindrome(string s) { int n = s.Length; var rh = new RollingHash(s); var irh = new RollingHash(string.Concat(s.Reverse())); ulong Hash(int l, int r) => rh[l, r]; ulong IHash(int l, int r) => irh[n - r, n - l]; bool Check(int l, int r) => Hash(l, r) == IHash(l, r); var ret = new List<(int l, int r)>(); for (int i = 0; i < n; i++) { int ok = 0, ng = Math.Min(i, n - 1 - i); if (Check(i - ng, i + ng + 1)) { ret.Add((i - ng, i + ng + 1)); continue; } while (ng - ok > 1) { int x = (ok + ng) / 2; if (Check(i - x, i + x + 1)) ok = x; else ng = x; } ret.Add((i - ok, i + ok + 1)); } for (int i = 0; i < n - 1; i++) { if (s[i] != s[i + 1]) continue; int ok = 0, ng = Math.Min(i, n - 2 - i); if (Check(i - ng, i + ng + 2)) { ret.Add((i - ng, i + ng + 2)); continue; } while (ng - ok > 1) { int x = (ok + ng) / 2; if (Check(i - x, i + x + 2)) ok = x; else ng = x; } ret.Add((i - ok, i + ok + 2)); } return ret; } } class RollingHash { const ulong MOD = (1UL << 61) - 1; const ulong MASK30 = (1UL << 30) - 1; const ulong MASK31 = (1UL << 31) - 1; const ulong MASK61 = (1UL << 61) - 1; const ulong POSITIVIZER = MOD * 4; int n; public static ulong b { get; private set; } public ulong[] hash, pow; public RollingHash(string s) { n = s.Length; if (b < (1 << 8)) b = (ulong)(new Random().Next(1 << 16) + (1 << 8)); hash = new ulong[n + 1]; pow = new ulong[n + 1]; pow[0] = 1; for (int i = 0; i < n; i++) { hash[i + 1] = Mod(Mul(hash[i], b) + s[i]); pow[i + 1] = Mod(Mul(pow[i], b)); } } public ulong this[int l, int r] { [MethodImpl(256)] get => Get(l, r); } [MethodImpl(256)] public ulong Get(int l, int r) => Mod(hash[r] + POSITIVIZER - Mul(hash[l], pow[r - l])); [MethodImpl(256)] public static ulong Mul(ulong a, ulong b) { ulong au = a >> 31; ulong ad = a & MASK31; ulong bu = b >> 31; ulong bd = b & MASK31; ulong mid = ad * bu + au * bd; ulong midu = mid >> 30; ulong midd = mid & MASK30; return au * bu * 2 + midu + (midd << 31) + ad * bd; } [MethodImpl(256)] public static ulong Mod(ulong x) { ulong xu = x >> 61; ulong xd = x & MASK61; ulong res = xu + xd; if (res >= MOD) res -= MOD; return res; } }