#572005 (C#(csc)) No.854 公平なりんご分配

提出ソース

結果

問題	No.854 公平なりんご分配
ユーザー	fairy_lettuce
提出日時	2020-10-28 08:59:29
言語	C#(csc) (csc 3.9.0)
結果	TLE
実行時間	-
コード長	14,866 bytes
コンパイル時間	1,283 ms
コンパイル使用メモリ	123,508 KB
実行使用メモリ	175,780 KB
最終ジャッジ日時	2024-07-21 22:11:22
合計ジャッジ時間	17,518 ms
ジャッジサーバーID （参考情報）	judge1 / judge2

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	42 ms 26,368 KB
testcase_01	AC	43 ms 20,992 KB
testcase_02	AC	46 ms 21,248 KB
testcase_03	AC	46 ms 21,504 KB
testcase_04	AC	47 ms 21,376 KB
testcase_05	AC	46 ms 20,992 KB
testcase_06	AC	41 ms 20,864 KB
testcase_07	AC	45 ms 21,376 KB
testcase_08	AC	44 ms 21,120 KB
testcase_09	AC	43 ms 21,248 KB
testcase_10	AC	46 ms 21,376 KB
testcase_11	AC	45 ms 21,120 KB
testcase_12	AC	48 ms 21,632 KB
testcase_13	AC	49 ms 22,016 KB
testcase_14	AC	47 ms 21,248 KB
testcase_15	AC	45 ms 21,376 KB
testcase_16	AC	47 ms 21,760 KB
testcase_17	AC	48 ms 21,888 KB
testcase_18	AC	45 ms 21,632 KB
testcase_19	AC	47 ms 21,632 KB
testcase_20	AC	44 ms 21,376 KB
testcase_21	AC	47 ms 21,632 KB
testcase_22	AC	60 ms 23,296 KB
testcase_23	AC	64 ms 25,984 KB
testcase_24	AC	70 ms 27,648 KB
testcase_25	AC	59 ms 24,192 KB
testcase_26	AC	71 ms 27,904 KB
testcase_27	AC	68 ms 27,264 KB
testcase_28	AC	65 ms 23,424 KB
testcase_29	AC	56 ms 23,168 KB
testcase_30	AC	60 ms 23,936 KB
testcase_31	AC	69 ms 27,392 KB
testcase_32	AC	350 ms 120,972 KB
testcase_33	AC	285 ms 74,768 KB
testcase_34	AC	494 ms 129,948 KB
testcase_35	AC	413 ms 126,164 KB
testcase_36	AC	191 ms 31,232 KB
testcase_37	AC	349 ms 121,492 KB
testcase_38	AC	241 ms 85,820 KB
testcase_39	AC	651 ms 175,780 KB
testcase_40	AC	302 ms 73,180 KB
testcase_41	AC	410 ms 121,668 KB
testcase_42	AC	464 ms 128,732 KB
testcase_43	AC	504 ms 124,644 KB
testcase_44	AC	519 ms 128,856 KB
testcase_45	AC	387 ms 52,096 KB
testcase_46	AC	606 ms 129,956 KB
testcase_47	AC	278 ms 75,528 KB
testcase_48	AC	441 ms 127,944 KB
testcase_49	AC	429 ms 130,048 KB
testcase_50	AC	333 ms 120,464 KB
testcase_51	AC	599 ms 128,080 KB
testcase_52	TLE	-
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	--	-
testcase_73	--	-
testcase_74	--	-
testcase_75	--	-
testcase_76	--	-
testcase_77	--	-
testcase_78	--	-
testcase_79	--	-
testcase_80	--	-
testcase_81	--	-
testcase_82	--	-
testcase_83	--	-
testcase_84	--	-
testcase_85	--	-
testcase_86	--	-
testcase_87	--	-
testcase_88	--	-
testcase_89	--	-
testcase_90	--	-
testcase_91	--	-
testcase_92	--	-
testcase_93	--	-

権限があれば一括ダウンロードができます

コンパイルメッセージ

Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

raw source code

using System;
using System.IO;
using System.Collections.Generic;
using System.Linq;
using System.Numerics;
using System.Diagnostics;
using System.Runtime.CompilerServices;

namespace FertiLib.Contest.E
{
	static class Program
	{
		public static void Solve(Scanner cin)
		{
			int n = cin.ReadInt();
			var a = cin.ReadIntArray(n);
			var q = cin.ReadInt();
			var prime = new List<int>();
			PrimeEnumerator.EnumPrime(2000, prime);
			var primedic = new SortedDictionary<int, int>();
			for (int i = 0; i < prime.Count; i++)
			{
				primedic.Add(prime[i], i);
			}
			var seg = new SegmentTree<int>[prime.Count];
			for (int i = 0; i < prime.Count; i++)
			{
				var div = new int[n];
				for (int j = 0; j < n; j++)
				{
					int t = 0;
					int tmp = a[j];
					while (tmp % prime[i] == 0)
					{
						t++;
						tmp /= prime[i];
					}
					div[j] = t;
				}
				seg[i] = new SegmentTree<int>(div, (x, y) => x + y, 0);
			}
			for (int i = 0; i < q; i++)
			{
				var (p, l, r) = cin.ReadValue<int, int, int>();
				bool ok = true;
				var factor = Factorizer.Factorize(p);
				foreach (var e in factor)
				{
					if (e.Key > 2000)
					{
						ok = false;
						break;
					}
					if (seg[primedic[(int)e.Key]].GetResult(l - 1, r) < e.Value)
					{
						ok = false;
						break;
					}
				}
				Console.WriteLine(ok ? "Yes" : "NO");
			}
		}

		public static void Main(string[] args)
		{
			var sw = new StreamWriter(Console.OpenStandardOutput()) { AutoFlush = false };
			Console.SetOut(sw);
			var cin = new Scanner();
			Solve(cin);
			Console.Out.Flush();
		}

		public static void YESNO(bool condition) => Console.WriteLine(condition ? "YES" : "NO");
		public static void YesNo(bool condition) => Console.WriteLine(condition ? "Yes" : "No");
		public static void yesno(bool condition) => Console.WriteLine(condition ? "yes" : "no");

		public static bool Chmax<T>(ref T a, T b) where T : IComparable<T>
		{
			if (a.CompareTo(b) >= 0) return false;
			a = b;
			return true;
		}
		public static bool Chmin<T>(ref T a, T b) where T : IComparable<T>
		{
			if (a.CompareTo(b) <= 0) return false;
			a = b;
			return true;
		}
	}

	public static class Factorizer
	{
		// http://miller-rabin.appspot.com/
		static readonly long[] baseSingle = { 126401071349994536 };                 // <= 291,831
		static readonly long[] baseDouble = { 336781006125, 9639812373923155 };     // <= 1,050,535,501
		static readonly long[] baseQuad = { 2, 2570940, 211991001, 3749873356 };    // <= 47,636,622,961,201
		static readonly long[] baseBig = { 2, 325, 9375, 28178, 450775, 9780504, 1795265022 };

		/// <summary>
		/// Factorizes n.
		/// </summary>
		/// <param name="n">The number to factorize.</param>
		/// <returns>Dictionary that has prime factors in Key, the number of each factor in Value.</returns>
		public static SortedDictionary<long, int> Factorize(long n)
		{
			var ret = new SortedDictionary<long, int>();
			var que = new Queue<long>();
			que.Enqueue(n);
			while (que.Count > 0)
			{
				var now = que.Dequeue();
				if (now == 1) continue;
				if (IsPrime(now))
				{
					if (ret.ContainsKey(now)) ret[now]++;
					else ret.Add(now, 1);
					continue;
				}

				long f = FindFactor(now);
				que.Enqueue(f);
				que.Enqueue(now / f);
			}

			return ret;
		}

		/// <summary>
		/// Tests if n is prime or not using Miller-Rabin Algorithm. Complexity: O(log^2 n)
		/// </summary>
		/// <param name="n">The number to test.</param>
		/// <returns>True if prime, False otherwise.</returns>
		public static bool IsPrime(long n)
		{
			if (n == 1) return false;
			if (n == 2) return true;

			long d = n - 1;
			int s = 0;
			while (d % 2 == 0)
			{
				d /= 2;
				s++;
			}

			long[] bases;
			if (n < 291831) bases = baseSingle;
			else if (n < 1050535501) bases = baseDouble;
			else if (n < 47636622961201) bases = baseQuad;
			else bases = baseBig;

			foreach (var e in bases)
			{
				if (!MillerRabinTest(e, d, n, s)) return false;
			}
			return true;
		}

		private static bool MillerRabinTest(long e, long d, long n, long s)
		{
			long pow = ModPow(e, d, n);
			if (pow == 1) return true;
			for (int i = 0; i < s; i++)
			{
				if (pow == n - 1) return true;
				pow = ModPow(pow, 2, n);
			}
			return false;
		}

		private static long ModPow(long a, long p, long mod)
		{
			if (mod <= int.MaxValue)
			{
				a %= mod;
				long ret = 1;
				for (long i = 1; i <= p; i *= 2)
				{
					if (p / i % 2 == 1)
					{
						ret *= a;
						ret %= mod;
					}

					a *= a;
					a %= mod;
				}
				return ret;
			}
			else
			{
				BigInteger retbig = 1;
				BigInteger abig = a % mod;
				for (BigInteger i = 1; i <= p; i *= 2)
				{
					if (p / i % 2 == 1)
					{
						retbig *= abig;
						retbig %= mod;
					}

					abig *= abig;
					abig %= mod;
				}
				return (long)retbig;
			}
		}

		/// <summary>
		/// Finds a factor of n by Pollard's ρ algorithm. Complexity: O(n^(1/4))
		/// </summary>
		/// <param name="n">The number to get a factor.</param>
		/// <returns>A prime factor of n.</returns>
		public static long FindFactor(long n)
		{
			if (n % 2 == 0) return 2;
			long i = 0;
			while (true)
			{
				i++;
				long x = i;
				long y = Next(i, n);
				while (true)
				{
					long g = GCD(x - y, n);
					if (g >= n) break;
					if (1 < g) return g;
					x = Next(x, n);
					y = Next(y, n);
					y = Next(y, n);
				}
			}
		}

		private static long Next(long x, long mod)
		{
			BigInteger tmp = x;
			tmp *= tmp;
			tmp++;
			tmp %= mod;
			return (long)tmp;
		}

		private static long GCD(long a, long b)
		{
			a = Math.Abs(a);
			b = Math.Abs(b);
			while (a > 0)
			{
				b %= a;
				var x = a;
				a = b;
				b = x;
			}
			return b;
		}
	}

	public static class PrimeEnumerator
	{
		public static bool[] EnumPrime(int max, List<int> prime = null)
		{
			var ret = new bool[max + 1];
			for (int i = 2; i < max + 1; i++)
			{
				ret[i] = true;
			}

			for (int i = 2; i * i <= max; i++)
			{
				if (!ret[i]) continue;
				for (int j = i * i; j <= max; j += i)
				{
					ret[j] = false;
				}
			}

			if (prime != null)
			{
				for (int i = 2; i < max + 1; i++)
				{
					if (ret[i]) prime.Add(i);
				}
			}

			return ret;
		}
	}

	public class SegmentTree<T>
	{
		public int Count { get; private set; }
		T[] tree;
		readonly Func<T, T, T> func;
		readonly T init;

		public SegmentTree(int size, Func<T, T, T> func, T init)
		{
			int n = 1;
			while (n <= size)
			{
				n *= 2;
			}
			Count = n;

			tree = Enumerable.Repeat(init, 2 * n - 1).ToArray();

			this.func = func;
			this.init = init;
		}

		public SegmentTree(IEnumerable<T> items, Func<T, T, T> func, T init) : this(items.Count(), func, init)
		{
			var array = items.ToArray();

			for (int i = 0; i < array.Length; i++)
			{
				tree[i + Count - 1] = array[i];
			}
			UpdateAll();
		}

		public T this[int index]
		{
			get => this.tree[index + Count - 1];
			private set => this.tree[index + Count - 1] = value;
		}

		public void Update(int index, T value)
		{
			index += Count - 1;
			tree[index] = value;
			while (index > 0)
			{
				index = (index - 1) / 2;
				tree[index] = func(tree[index * 2 + 1], tree[index * 2 + 2]);
			}
		}

		void UpdateAll()
		{
			int n = Count / 2;
			while (n > 0)
			{
				for (int i = 0; i < n; i++)
				{
					tree[n - 1 + i] = func(tree[2 * (n - 1 + i) + 1], tree[2 * (n - 1 + i) + 2]);
				}
				n /= 2;
			}
		}

		// Attention: [l, r). This result does not include r-th item.
		public T GetResult(int l, int r) => GetResult(l, r, 0, 0, Count);

		T GetResult(int l, int r, int k, int a, int b)
		{
			if (b <= l || r <= a) return init;
			if (l <= a && b <= r) return tree[k];

			var vl = GetResult(l, r, 2 * k + 1, a, (a + b) / 2);
			var vr = GetResult(l, r, 2 * k + 2, (a + b) / 2, b);

			return func(vl, vr);
		}
	}

	static class Extention
	{
		public static string Join<T>(this IEnumerable<T> x, string separator = "") => string.Join(separator, x);

		public static int UpperBound<T>(this IList<T> list, T value) => list.BinarySearch(value, true, 0, list.Count, Comparer<T>.Default);
		public static int LowerBound<T>(this IList<T> list, T value) => list.BinarySearch(value, false, 0, list.Count, Comparer<T>.Default);
		public static int BinarySearch<T>(this IList<T> list, T value, bool isUpperBound, int index, int length, Comparer<T> comparer)
		{
			var ng = index - 1;
			var ok = index + length;
			while (ok - ng > 1)
			{
				var mid = ng + (ok - ng) / 2;
				var res = comparer.Compare(list[mid], value);
				if (res < 0 || (isUpperBound && res == 0)) ng = mid;
				else ok = mid;
			}
			return ok;
		}
	}

	class Scanner
	{
		string[] s;
		int i;

		char[] separator = new char[] { ' ' };

		public Scanner()
		{
			s = new string[0];
			i = 0;
		}

		public string Read() => ReadString();

		public string ReadString()
		{
			if (i < s.Length) return s[i++];
			string st = Console.ReadLine();
			while (st == "") st = Console.ReadLine();
			s = st.Split(separator, StringSplitOptions.RemoveEmptyEntries);
			if (s.Length == 0) return ReadString();
			i = 0;
			return s[i++];
		}

		public string[] ReadStringArray(int N)
		{
			string[] Array = new string[N];
			for (int i = 0; i < N; i++)
			{
				Array[i] = ReadString();
			}
			return Array;
		}

		public int ReadInt()
		{
			return int.Parse(ReadString());
		}

		public int[] ReadIntArray(int N, int add = 0)
		{
			int[] Array = new int[N];
			for (int i = 0; i < N; i++)
			{
				Array[i] = ReadInt() + add;
			}
			return Array;
		}

		public long ReadLong()
		{
			return long.Parse(ReadString());
		}

		public long[] ReadLongArray(int N, long add = 0)
		{
			long[] Array = new long[N];
			for (int i = 0; i < N; i++)
			{
				Array[i] = ReadLong() + add;
			}
			return Array;
		}

		public double ReadDouble()
		{
			return double.Parse(ReadString());
		}

		public double[] ReadDoubleArray(int N, double add = 0)
		{
			double[] Array = new double[N];
			for (int i = 0; i < N; i++)
			{
				Array[i] = ReadDouble() + add;
			}
			return Array;
		}

		public T1 ReadValue<T1>() => (T1)Convert.ChangeType(ReadString(), typeof(T1));

		public (T1, T2) ReadValue<T1, T2>()
		{
			var inputs = ReadStringArray(2);
			var v1 = (T1)Convert.ChangeType(inputs[0], typeof(T1));
			var v2 = (T2)Convert.ChangeType(inputs[1], typeof(T2));
			return (v1, v2);
		}

		public (T1, T2, T3) ReadValue<T1, T2, T3>()
		{
			var inputs = ReadStringArray(3);
			var v1 = (T1)Convert.ChangeType(inputs[0], typeof(T1));
			var v2 = (T2)Convert.ChangeType(inputs[1], typeof(T2));
			var v3 = (T3)Convert.ChangeType(inputs[2], typeof(T3));
			return (v1, v2, v3);
		}

		public (T1, T2, T3, T4) ReadValue<T1, T2, T3, T4>()
		{
			var inputs = ReadStringArray(4);
			var v1 = (T1)Convert.ChangeType(inputs[0], typeof(T1));
			var v2 = (T2)Convert.ChangeType(inputs[1], typeof(T2));
			var v3 = (T3)Convert.ChangeType(inputs[2], typeof(T3));
			var v4 = (T4)Convert.ChangeType(inputs[3], typeof(T4));
			return (v1, v2, v3, v4);
		}

		public (T1, T2, T3, T4, T5) ReadValue<T1, T2, T3, T4, T5>()
		{
			var inputs = ReadStringArray(5);
			var v1 = (T1)Convert.ChangeType(inputs[0], typeof(T1));
			var v2 = (T2)Convert.ChangeType(inputs[1], typeof(T2));
			var v3 = (T3)Convert.ChangeType(inputs[2], typeof(T3));
			var v4 = (T4)Convert.ChangeType(inputs[3], typeof(T4));
			var v5 = (T5)Convert.ChangeType(inputs[4], typeof(T5));
			return (v1, v2, v3, v4, v5);
		}

		public (T1, T2, T3, T4, T5, T6) ReadValue<T1, T2, T3, T4, T5, T6>()
		{
			var inputs = ReadStringArray(6);
			var v1 = (T1)Convert.ChangeType(inputs[0], typeof(T1));
			var v2 = (T2)Convert.ChangeType(inputs[1], typeof(T2));
			var v3 = (T3)Convert.ChangeType(inputs[2], typeof(T3));
			var v4 = (T4)Convert.ChangeType(inputs[3], typeof(T4));
			var v5 = (T5)Convert.ChangeType(inputs[4], typeof(T5));
			var v6 = (T6)Convert.ChangeType(inputs[5], typeof(T6));
			return (v1, v2, v3, v4, v5, v6);
		}

		public (T1, T2, T3, T4, T5, T6, T7) ReadValue<T1, T2, T3, T4, T5, T6, T7>()
		{
			var inputs = ReadStringArray(7);
			var v1 = (T1)Convert.ChangeType(inputs[0], typeof(T1));
			var v2 = (T2)Convert.ChangeType(inputs[1], typeof(T2));
			var v3 = (T3)Convert.ChangeType(inputs[2], typeof(T3));
			var v4 = (T4)Convert.ChangeType(inputs[3], typeof(T4));
			var v5 = (T5)Convert.ChangeType(inputs[4], typeof(T5));
			var v6 = (T6)Convert.ChangeType(inputs[5], typeof(T6));
			var v7 = (T7)Convert.ChangeType(inputs[6], typeof(T7));
			return (v1, v2, v3, v4, v5, v6, v7);
		}

		public T1[] ReadValueArray<T1>(int N)
		{
			var v1 = new T1[N];
			for (int i = 0; i < N; i++)
			{
				v1[i] = ReadValue<T1>();
			}
			return v1;
		}

		public (T1[], T2[]) ReadValueArray<T1, T2>(int N)
		{
			var (v1, v2) = (new T1[N], new T2[N]);
			for (int i = 0; i < N; i++)
			{
				var (t1, t2) = ReadValue<T1, T2>();
				v1[i] = t1;
				v2[i] = t2;
			}
			return (v1, v2);
		}

		public (T1[], T2[], T3[]) ReadValueArray<T1, T2, T3>(int N)
		{
			var (v1, v2, v3) = (new T1[N], new T2[N], new T3[N]);
			for (int i = 0; i < N; i++)
			{
				var (t1, t2, t3) = ReadValue<T1, T2, T3>();
				v1[i] = t1;
				v2[i] = t2;
				v3[i] = t3;
			}
			return (v1, v2, v3);
		}

		public (T1[], T2[], T3[], T4[]) ReadValueArray<T1, T2, T3, T4>(int N)
		{
			var (v1, v2, v3, v4) = (new T1[N], new T2[N], new T3[N], new T4[N]);
			for (int i = 0; i < N; i++)
			{
				var (t1, t2, t3, t4) = ReadValue<T1, T2, T3, T4>();
				v1[i] = t1;
				v2[i] = t2;
				v3[i] = t3;
				v4[i] = t4;
			}
			return (v1, v2, v3, v4);
		}

		public (T1[], T2[], T3[], T4[], T5[]) ReadValueArray<T1, T2, T3, T4, T5>(int N)
		{
			var (v1, v2, v3, v4, v5) = (new T1[N], new T2[N], new T3[N], new T4[N], new T5[N]);
			for (int i = 0; i < N; i++)
			{
				var (t1, t2, t3, t4, t5) = ReadValue<T1, T2, T3, T4, T5>();
				v1[i] = t1;
				v2[i] = t2;
				v3[i] = t3;
				v4[i] = t4;
				v5[i] = t5;
			}
			return (v1, v2, v3, v4, v5);
		}

		public (T1[], T2[], T3[], T4[], T5[], T6[]) ReadValueArray<T1, T2, T3, T4, T5, T6>(int N)
		{
			var (v1, v2, v3, v4, v5, v6) = (new T1[N], new T2[N], new T3[N], new T4[N], new T5[N], new T6[N]);
			for (int i = 0; i < N; i++)
			{
				var (t1, t2, t3, t4, t5, t6) = ReadValue<T1, T2, T3, T4, T5, T6>();
				v1[i] = t1;
				v2[i] = t2;
				v3[i] = t3;
				v4[i] = t4;
				v5[i] = t5;
				v6[i] = t6;
			}
			return (v1, v2, v3, v4, v5, v6);
		}

		public (T1[], T2[], T3[], T4[], T5[], T6[], T7[]) ReadValueArray<T1, T2, T3, T4, T5, T6, T7>(int N)
		{
			var (v1, v2, v3, v4, v5, v6, v7) = (new T1[N], new T2[N], new T3[N], new T4[N], new T5[N], new T6[N], new T7[N]);
			for (int i = 0; i < N; i++)
			{
				var (t1, t2, t3, t4, t5, t6, t7) = ReadValue<T1, T2, T3, T4, T5, T6, T7>();
				v1[i] = t1;
				v2[i] = t2;
				v3[i] = t3;
				v4[i] = t4;
				v5[i] = t5;
				v6[i] = t6;
				v7[i] = t7;
			}
			return (v1, v2, v3, v4, v5, v6, v7);
		}
	}
}

yukicoder

結果

テストケース

コンパイルメッセージ

ソースコード