結果

問題 No.1409 Simple Math in yukicoder
ユーザー kakel-sankakel-san
提出日時 2024-07-15 23:33:12
言語 C#
(.NET 8.0.404)
結果
AC  
実行時間 148 ms / 2,000 ms
コード長 3,137 bytes
コンパイル時間 10,188 ms
コンパイル使用メモリ 171,308 KB
実行使用メモリ 194,140 KB
最終ジャッジ日時 2024-07-15 23:33:57
合計ジャッジ時間 38,814 ms
ジャッジサーバーID
(参考情報)
judge1 / judge4
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 2
other AC * 58
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (106 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/

ソースコード

diff #
プレゼンテーションモードにする

using System;
using static System.Console;
using System.Linq;
using System.Collections.Generic;
class Program
{
static int NN => int.Parse(ReadLine());
static int[] NList => ReadLine().Split().Select(int.Parse).ToArray();
static int[][] NArr(long n) => Enumerable.Repeat(0, (int)n).Select(_ => NList).ToArray();
static int[] LList(long n) => Enumerable.Repeat(0, (int)n).Select(_ => int.Parse(ReadLine())).ToArray();
public static void Main()
{
Solve();
}
static void Solve()
{
var r = new Random();
var t = NN;
var ans = new int[t][];
for (var u = 0; u < t; ++u)
{
var c = NList;
var (v, x) = (c[0], c[1]);
ans[u] = Math(r, v, x);
}
WriteLine(string.Join("\n", ans.Select(ai => string.Join(" ", ai))));
}
static int[] Math(Random r, int v, int x)
{
if (x == 1) return new int[] { 1 };
if (v == 1)
{
var _ans = new int[x];
for (var i = 0; i < x; ++i) _ans[i] = i + 1;
return _ans;
}
var p = PDiv(v);
foreach (var kv in PDiv(x)) if (!p.ContainsKey(kv.Key)) p[kv.Key] = kv.Value;
var b = 0;
while (true)
{
b = r.Next(v * x) + 1;
var flg = true;
foreach (var kv in p)
{
if (Exp(b, v * x / kv.Key, v * x + 1) == 1)
{
flg = false;
break;
}
}
if (flg) break;
}
var ans = new int[x];
for (var i = 0; i < x; ++i)
{
ans[i] = (int)Exp(b, i * v, v * x + 1);
}
Array.Sort(ans);
return ans;
}
static Dictionary<long, long> PDiv(long n)
{
var dic = new Dictionary<long, long>();
var tmp = n;
while (tmp % 2 == 0)
{
tmp /= 2;
if (dic.ContainsKey(2)) ++dic[2];
else dic.Add(2, 1);
}
for (var p = 3L; p * p <= n; p += 2)
{
while (tmp % p == 0)
{
tmp /= p;
if (dic.ContainsKey(p)) ++dic[p];
else dic.Add(p, 1);
}
if (tmp == 1) break;
}
if (tmp > 1) dic.Add(tmp, 1);
return dic;
}
static long Exp(long n, long p, int mod)
{
long _n = n % mod;
var _p = p;
var result = 1L;
if ((_p & 1) == 1) result *= _n;
while (_p > 0)
{
_n = _n * _n % mod;
_p >>= 1;
if ((_p & 1) == 1) result = result * _n % mod;
}
return result;
}
static bool IsPrime(int n)
{
if (n == 1) return false;
if (n == 2) return true;
if (n % 2 == 0) return false;
int i = 3;
while (i * i <= n)
{
if (n % i == 0) return false;
i = i + 2;
}
return true;
}
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0