結果
| 問題 | No.1068 #いろいろな色 / Red and Blue and more various colors (Hard) |
| コンテスト | |
| ユーザー |
eSeF
|
| 提出日時 | 2020-07-13 01:38:30 |
| 言語 | C#(csc) (csc 3.9.0) |
| 結果 |
TLE
|
| 実行時間 | - |
| コード長 | 12,668 bytes |
| 記録 | |
| コンパイル時間 | 1,348 ms |
| コンパイル使用メモリ | 122,184 KB |
| 実行使用メモリ | 117,680 KB |
| 最終ジャッジ日時 | 2024-11-06 12:36:04 |
| 合計ジャッジ時間 | 9,428 ms |
|
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 10 TLE * 1 -- * 18 |
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc) Copyright (C) Microsoft Corporation. All rights reserved.
ソースコード
using System;
using System.Collections.Generic;
using System.Linq;
using System.IO;
using System.Text;
using System.Numerics;
using System.Threading;
using System.Runtime.CompilerServices;
using System.Diagnostics;
using static System.Math;
using static System.Array;
using static AtCoder.Sc_out;
using static AtCoder.Tool;
using static AtCoder.ModInt;
namespace AtCoder
{
class AC
{
//const int MOD = 1000000007;
const int MOD = 998244353;
const int INF = int.MaxValue / 2;
const long SINF = long.MaxValue / 2;
const double EPS = 1e-8;
static readonly int[] dI = { 0, 1, 0, -1, 1, -1, -1, 1 };
static readonly int[] dJ = { 1, 0, -1, 0, 1, 1, -1, -1 };
//static List<List<int>> G = new List<List<int>>();
//static List<List<Edge>> G = new List<List<Edge>>();
//static List<Edge> E = new List<Edge>();
static void Main(string[] args)
{
var sw = new StreamWriter(Console.OpenStandardOutput()) { AutoFlush = false }; Console.SetOut(sw);
var th = new Thread(Run, 1 << 26);
th.Start();
th.Join();
//Run();
Console.Out.Flush();
}
static void Run()
{
int Testcase = 1;
for (var _ = 0; _ < Testcase; _++) Solve();
}
static void Solve()
{
Cin.Input(out int N, out int Q);
var A = Cin.ReadSplitLong;
var B = Cin.ReadSplitInt;
var calc = new Math_NTT(MOD, 3);
List<long> f(List<long> x, List<long> y)
{
return calc.Convolution_NTT(x.ToArray(), y.ToArray()).ToList();
}
var seg = new SegmentTree<List<long>>(N + 1, f, new List<long>() { 1 });
for (var i = 0; i < N; i++) seg.Assign(i, new List<long>() { (A[i] - 1 + MOD) % MOD, 1 });
seg.Set_All();
var ans = seg.Query_Answer(0, N);
for (var i = 0; i < Q; i++) OutL(ans[B[i]]);
}
public struct Edge
{
public int from;
public int to;
public long dist;
public Edge(int t, long c)
{
from = -1;
to = t;
dist = c;
}
public Edge(int f, int t, long c)
{
from = f;
to = t;
dist = c;
}
}
}
public class Math_NTT
{
private long mod;
private long root;
public Math_NTT(long prime_mod, long premitive_root) { mod = prime_mod; root = premitive_root; }
[MethodImpl(MethodImplOptions.AggressiveInlining)]
long modpow(long x, long nn) => nn == 0 ? 1 : (nn % 2 == 0 ? modpow((x * x) % mod, nn / 2) : (x * modpow(x, nn - 1)) % mod);
[MethodImpl(MethodImplOptions.AggressiveInlining)]
long inverse(long x) => modpow(x, mod - 2);
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public void NTT(ref long[] a, bool inv = false)
{
int n = a.Length;
int h = 0;
for (var i = 0; (1 << i) < n; i++) h++;
//バタフライ演算準備
for (var i = 0; i < n; i++)
{
int j = 0;
for (var k = 0; k < h; k++) j |= ((i >> k) & 1) << (h - 1 - k);
//2進法表記が逆の数字とスワップする
if (i < j)
{
var kep = a[i];
a[i] = a[j];
a[j] = kep;
}
}
//バタフライ演算パート
for (var b = 1; b < n; b <<= 1)
{
long W = modpow(root, (inv ? (mod - 1) - (mod - 1) / (b << 1) : (mod - 1) / (b << 1)));
long z = 1;
for (var j = 0; j < b; j++)
{
for (var k = 0; k < n; k += (b << 1))
{
var s = a[k + j];
var t = (a[k + j + b] * z) % mod;
a[k + j] = (s + t) % mod;
a[k + j + b] = (s - t + mod) % mod;
}
z = (z * W) % mod;
}
}
if (inv)
{
for (var i = 0; i < n; i++) a[i] = (a[i] * modpow(n, mod - 2)) % mod;
}
}
public long[] Convolution_NTT(long[] a, long[] b)
{
//memo: mod 998244353 => primitive_root = 3
// mod 924844033 => root = 5
int N = a.Length + b.Length - 1;
int t = 1;
while (t < N) t <<= 1;
var nxa = new long[t];
var nxb = new long[t];
for (var i = 0; i < t; i++)
{
nxa[i] = i < a.Length ? a[i] : 0;
nxb[i] = i < b.Length ? b[i] : 0;
}
NTT(ref nxa);
NTT(ref nxb);
for (var i = 0; i < t; i++) nxa[i] = (nxa[i] * nxb[i]) % mod;
NTT(ref nxa, true);
var ret = new long[N];
for (var i = 0; i < N; i++) ret[i] = nxa[i];
return ret;
}
}
public class SegmentTree<T>
{
int n;
T[] Tree;
Func<T, T, T> f;
T ex;
public SegmentTree(int size, Func<T, T, T> fun, T exvalue)
{
f = fun;
ex = exvalue;
n = 1;
while (n < size) n <<= 1;
Tree = new T[(n << 1) - 1];
for (var i = 0; i < Tree.Length; i++) Tree[i] = ex;
}
public SegmentTree(int size, Func<T, T, T> f, T exvalue, T initial) : this(size, f, exvalue)
{
for (var i = 0; i < size; i++) Tree[i + n - 1] = initial;
Set_All();
}
public void Set_All()
{
for (var i = n - 2; i >= 0; i--)
{
Tree[i] = f(Tree[(i << 1) + 1], Tree[(i << 1) + 2]);
}
}
public void Assign(int ind, T val)
{
Tree[ind + n - 1] = val;
}
public void Update(int ind)
{
int i = ind + n - 1;
while (i > 0)
{
i = (i - 1) >> 1;
Tree[i] = f(Tree[(i << 1) + 1], Tree[(i << 1) + 2]);
}
}
public void Query_Update(int ind, T val)
{
Assign(ind, val);
Update(ind);
}
public T Query_Answer(int l, int r) //[l,r)
{
return Scanning(l, r, 0, 0, n);
}
private T Scanning(int l, int r, int ind, int ind_l, int ind_r)
{
if (ind_r <= l || r <= ind_l) { return ex; }
if (l <= ind_l && ind_r <= r) { return Tree[ind]; }
return f(Scanning(l, r, (ind << 1) + 1, ind_l, (ind_l + ind_r) >> 1), Scanning(l, r, (ind << 1) + 2, (ind_l + ind_r) >> 1, ind_r));
}
public T Look(int ind)
{
return Tree[ind + n - 1];
}
public void Debug_Display_All()
{
for (var i = 0; i < n; i++) Console.Write($"{Tree[i + n - 1]} ");
Console.WriteLine();
}
}
struct ModInt
{
public long value;
private const int MOD = 1000000007;
//private const int MOD = 998244353;
public ModInt(long value) { this.value = value; }
public static implicit operator ModInt(long a)
{
var ret = a % MOD;
return new ModInt(ret < 0 ? (ret + MOD) : ret);
}
public static ModInt operator +(ModInt a, ModInt b) => (a.value + b.value);
public static ModInt operator -(ModInt a, ModInt b) => (a.value - b.value);
public static ModInt operator *(ModInt a, ModInt b) => (a.value * b.value);
public static ModInt operator /(ModInt a, ModInt b) => a * Modpow(b, MOD - 2);
public static ModInt operator <<(ModInt a, int n) => (a.value << n);
public static ModInt operator >>(ModInt a, int n) => (a.value >> n);
public static ModInt operator ++(ModInt a) => a.value + 1;
public static ModInt operator --(ModInt a) => a.value - 1;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static ModInt Modpow(ModInt a, long n)
{
var k = a;
ModInt ret = 1;
while (n > 0)
{
if ((n & 1) != 0) ret *= k;
k *= k;
n >>= 1;
}
return ret;
}
private static readonly List<long> Factorials = new List<long>() { 1 };
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static ModInt Fac(long n)
{
for (var i = Factorials.Count(); i <= n; i++)
{
Factorials.Add((Factorials[i - 1] * i) % MOD);
}
return Factorials[(int)n];
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public static ModInt nCr(long n, long r)
{
return n < r ? 0 : Fac(n) / (Fac(r) * Fac(n - r));
}
public static explicit operator int(ModInt a) => (int)a.value;
}
static class Cin
{
public static string[] ReadSplit => Console.ReadLine().Split();
public static int[] ReadSplitInt => ConvertAll(Console.ReadLine().Split(), int.Parse);
public static long[] ReadSplitLong => ConvertAll(Console.ReadLine().Split(), long.Parse);
public static double[] ReadSplit_Double => ConvertAll(Console.ReadLine().Split(), double.Parse);
public static string Str => Console.ReadLine();
public static int Int => int.Parse(Console.ReadLine());
public static long Long => long.Parse(Console.ReadLine());
public static double Double => double.Parse(Console.ReadLine());
public static T Conv<T>(string input)
{
if (typeof(T).Equals(typeof(ModInt)))
{
return (T)(dynamic)(long.Parse(input));
}
return (T)Convert.ChangeType(input, typeof(T));
}
public static void Input<T>(out T a) => a = Conv<T>(Console.ReadLine());
public static void Input<T, U>(out T a, out U b)
{ var q = ReadSplit; a = Conv<T>(q[0]); b = Conv<U>(q[1]); }
public static void Input<T, U, V>(out T a, out U b, out V c)
{ var q = ReadSplit; a = Conv<T>(q[0]); b = Conv<U>(q[1]); c = Conv<V>(q[2]); }
public static void Input<T, U, V, W>(out T a, out U b, out V c, out W d)
{ var q = ReadSplit; a = Conv<T>(q[0]); b = Conv<U>(q[1]); c = Conv<V>(q[2]); d = Conv<W>(q[3]); }
public static void Input<T, U, V, W, X>(out T a, out U b, out V c, out W d, out X e)
{ var q = ReadSplit; a = Conv<T>(q[0]); b = Conv<U>(q[1]); c = Conv<V>(q[2]); d = Conv<W>(q[3]); e = Conv<X>(q[4]); }
}
static class Sc_out
{
public static void OutL(object s) => Console.WriteLine(s);
public static void Out_Sep<T>(IEnumerable<T> s) => Console.WriteLine(string.Join(" ", s));
public static void Out_Sep<T>(IEnumerable<T> s, string sep) => Console.WriteLine(string.Join($"{sep}", s));
public static void Out_Sep(params object[] s) => Console.WriteLine(string.Join(" ", s));
public static void Out_One(object s) => Console.Write($"{s} ");
public static void Out_One(object s, string sep) => Console.Write($"{s}{sep}");
public static void Endl() => Console.WriteLine();
}
public static class Tool
{
static public void Initialize<T>(ref T[] array, T initialvalue)
{
array = ConvertAll(array, x => initialvalue);
}
static public void Swap<T>(ref T a, ref T b)
{
T keep = a;
a = b;
b = keep;
}
static public void Display<T>(T[,] array2d, int n, int m)
{
for (var i = 0; i < n; i++)
{
for (var j = 0; j < m; j++)
{
Console.Write($"{array2d[i, j]} ");
}
Console.WriteLine();
}
}
static public long Gcd(long a, long b)
{
if (a == 0 || b == 0) return Max(a, b);
return a % b == 0 ? b : Gcd(b, a % b);
}
static public long LPow(int a, int b) => (long)Pow(a, b);
static public bool Bit(long x, int dig) => ((1L << dig) & x) != 0;
static public int Sig(long a) => a == 0 ? 0 : (int)(a / Abs(a));
}
}
eSeF