結果

問題 No.3651 K-th Sum of Divisors
コンテスト
ユーザー tobisatis
提出日時 2026-08-28 21:42:21
言語 C#
(.NET 10.0.201)
コンパイル:
dotnet_c
実行:
/usr/bin/dotnet_wrap
結果
AC  
実行時間 91 ms / 2,000 ms
+ 599µs
コード長 1,905 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 12,543 ms
コンパイル使用メモリ 171,216 KB
実行使用メモリ 219,288 KB
最終ジャッジ日時 2026-08-28 21:42:43
合計ジャッジ時間 19,656 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 55
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (111 ミリ秒)。
  main -> /home/judge/data/code/bin/Release/net10.0/main.dll
  main -> /home/judge/data/code/bin/Release/net10.0/publish/

ソースコード

diff #
raw source code 

#nullable enable

#region
var (_input, _iter) = (Array.Empty<string>(), 0);
T I<T>() where T : IParsable<T>
{
    while (_iter >= _input.Length) (_input, _iter) = (Console.ReadLine()!.Trim().Split(' '), 0);
    return T.Parse(_input[_iter++], null);
}
#endregion

const int P = 100003;

var n = I<int>();
var k = I<long>();

if (k == 1)
{
    Console.WriteLine(n);
    return;
}

var prime = new Prime(P);
var divall = prime.DivisorsAll();
var fz = new int[P + 1];
for (var i = 1; i <= P; i++)
{
    var divz = divall[i];
    foreach (var div in divz) fz[i] = (fz[i] + div) % P;
}
var ans = 0;
for (var i = 1; i * i <= n; i++)
{
    if (n % i != 0) continue;
    ans += i;
    if (i * i != n) ans += n / i;
    ans %= P;
}
k -= 2;
while (k > 0)
{
    if ((k & 1) != 0) ans = fz[ans];
    var nfz = new int[P + 1];
    for (var i = 1; i <= P; i++) nfz[i] = fz[fz[i]];
    fz = nfz;
    k >>= 1;
}
Console.WriteLine(ans);

class Prime
{
    public readonly int[] Sieve, Primes;

    public Prime(int n)
    {
        if (n <= 3) n = 3;
        var sieve = new int[n + 1];
        var (d, i) = (2, 5);
        var primes = new List<int>(){ 2, 3 };
        sieve[1] = 1;
        for (var j = 2; j <= n; j += 2) sieve[j] = 2;
        for (var j = 3; j <= n; j += 3) sieve[j] = 3;
        while (i <= n)
        {
            if (sieve[i] == 0)
            {
                primes.Add(i);
                for (var j = i; j <= n; j += i) sieve[j] = i;
            }
            i += d;
            d ^= 6;
        }
        (Sieve, Primes) = (sieve, primes.ToArray());
    }
}

static class PrimeExtensions
{
    public static List<int>[] DivisorsAll(this Prime prime)
    {
        var l = prime.Sieve.Length;
        var res = new List<int>[l];
        for (var i = 0; i < l; i++) res[i] = new();
        for (var i = 1; i < l; i++) for (var j = i; j < l; j += i) res[j].Add(i);
        return res;
    }

}
0