結果

問題 No.1409 Simple Math in yukicoder
ユーザー kakel-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;
    }
}
0