結果

問題 No.3688 LCM Sum
コンテスト
ユーザー tobisatis
提出日時 2026-09-06 19:59:44
言語 C#
(.NET 10.0.400 + ACL)
コンパイル:
dotnet_c
実行:
/usr/bin/dotnet_wrap
結果
AC  
実行時間 2,822 ms / 3,000 ms
+ 757µs
コード長 2,289 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 12,459 ms
コンパイル使用メモリ 172,536 KB
実行使用メモリ 690,872 KB
最終ジャッジ日時 2026-09-06 20:00:34
合計ジャッジ時間 38,556 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 11
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (102 ミリ秒)。
  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 Mod = 998244353;

var n = I<long>();
var m = I<long>();
var s = (int)long.Max(n, m);

var prime = new Prime(s);
var f = prime.Moebius();
var fi = new long[s + 1];
for (var i = 1L; i <= s; i++) fi[i] = i * f[i] % Mod;
prime.DivisorTransform(fi);
var ans = 0L;
for (var t = 1; t <= s; t++)
{
    var nc = n / t;
    var mc = m / t;
    var ns = nc * (nc + 1) / 2 % Mod;
    var ms = mc * (mc + 1) / 2 % Mod;
    var z = ns * ms % Mod * t % Mod * fi[t] % Mod;
    ans = (ans + Mod + z) % Mod;
}
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 int[] Moebius(this Prime prime)
    {
        var sieve = prime.Sieve;
        var l = sieve.Length;
        var res = new int[l];
        for (var i = 0; i < l; i++) res[i] = 1;
        var primes = prime.Primes;
        foreach (var p in primes)
        {
            for (var i = p; i < l; i += p) res[i] = -res[i];
            if ((long)p * p >= l) continue;
            var q = p * p;
            for (var i = q; i < l; i += q) res[i] = 0;
        }
        return res;
    }

    public static void DivisorTransform(this Prime prime, IList<long> values)
    {
        var sieve = prime.Sieve;
        for (var p = 2; p < sieve.Length; p++)
        {
            if (sieve[p] != p) continue;
            for (var i = 1; i * p < values.Count; i++) values[i * p] += values[i];
        }
    }
}
0