結果

問題 No.3478 XOR-Folding Primes
コンテスト
ユーザー kakel-san
提出日時 2026-03-20 23:40:26
言語 C#
(.NET 10.0.102)
コンパイル:
dotnet_c
実行:
/usr/bin/dotnet_wrap
結果
AC  
実行時間 1,135 ms / 4,000 ms
コード長 2,954 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 12,894 ms
コンパイル使用メモリ 172,404 KB
実行使用メモリ 324,504 KB
最終ジャッジ日時 2026-03-20 23:40:52
合計ジャッジ時間 18,071 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 8
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (85 ミリ秒)。
  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 

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();
    public static void Main()
    {
        Solve();
    }
    static void Solve()
    {
        var prime = Enumerable.Repeat(true, 10_000_001).ToArray();
        prime[0] = false;
        prime[1] = false;
        for (var i = 2; i < prime.Length; ++i)
        {
            if (!prime[i]) continue;
            for (var j = i * 2; j < prime.Length; j += i) prime[j] = false;
        }
        var pcount = new int[prime.Length];
        for (var i = 2; i < prime.Length; ++i) pcount[i] = pcount[i - 1] + (prime[i] ? 1 : 0);
        var a2 = new int[prime.Length];
        for (var i = 5; i < prime.Length; ++i)
        {
            a2[i] = a2[i - 1];
            if ((i ^ 2) < i && prime[i] && prime[i - 2]) a2[i] += 2;
        }
        var t = NN;
        var ans = new long[t];
        for (var u = 0; u < t; ++u)
        {
            var c = NList;
            var (n, m) = (c[0], c[1]);
            ans[u] = Fold(n, m, pcount, a2);
        }
        WriteLine(string.Join("\n", ans));
    }
    static long Fold(int n, int m, int[] pcount, int[] a2)
    {
        if (n == 1) return pcount[m];
        if (m < 7) return 0;
        var a = new long[][] { new long[] { 1, a2[m] }};
        var x = new long[][]
        {
            new long[] { 0, a2[m] },
            new long[] { 1, 1 }
        };
        var b = Mul(a, Pow(x, n - 1));
        return (b[0][0] + b[0][1]) % mod;
    }
    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 int mod = 998_244_353;
    // 行列の累乗
    static long[][] Pow(long[][] m, long k)
    {
        var multi = m;
        var r = new long[m.Length][];
        for (var i = 0; i < m.Length; ++i)
        {
            r[i] = new long[m.Length];
            r[i][i] = 1;
        }
        while (k > 0)
        {
            if ((k & 1) == 1) r = Mul(r, multi);
            multi = Mul(multi, multi);
            k >>= 1;
        }
        return r;
    }
    // 行列の積
    static long[][] Mul(long[][] x, long[][] y)
    {
        var r = new long[x.Length][];
        for (var i = 0; i < x.Length; ++i) r[i] = new long[y[0].Length];
        for (var i = 0; i < x.Length; ++i) for (var k = 0; k < x[0].Length; ++k)
        {
            for (var j = 0; j < y[0].Length; ++j) r[i][j] = (r[i][j] + x[i][k] * y[k][j]) % mod;
        }
        return r;
    }
}
0