結果

問題 No.3687 Coprime Count
コンテスト
ユーザー tobisatis
提出日時 2026-09-06 18:12:52
言語 C#
(.NET 10.0.400 + ACL)
コンパイル:
dotnet_c
実行:
/usr/bin/dotnet_wrap
結果
AC  
実行時間 467 ms / 2,000 ms
+ 824µs
コード長 1,668 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 5,274 ms
コンパイル使用メモリ 172,388 KB
実行使用メモリ 282,044 KB
最終ジャッジ日時 2026-09-06 18:13:15
合計ジャッジ時間 10,772 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge3_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (89 ミリ秒)。
  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

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 ans = 0L;
for (var i = 1; i <= s; i++) ans += f[i] * (n / i) * (m / i);
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;
    }

}
0