結果

問題 No.8030 ミラー・ラビン素数判定法のテスト
ユーザー iwkjosec
提出日時 2019-10-11 00:19:51
言語 C#(csc)
(csc 3.9.0)
結果
WA  
(最新)
AC  
(最初)
実行時間 -
コード長 1,265 bytes
コンパイル時間 2,551 ms
コンパイル使用メモリ 113,340 KB
実行使用メモリ 30,224 KB
最終ジャッジ日時 2024-11-18 17:08:09
合計ジャッジ時間 7,627 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 9 WA * 1
権限があれば一括ダウンロードができます
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc)
Copyright (C) Microsoft Corporation. All rights reserved.

ソースコード

diff # 

using System;
using System.Collections.Generic;
using System.Linq;
using System.Numerics;
using static System.Console;
using static System.Math;

class Program
{
    static void Main()
    {
        var n = long.Parse(ReadLine());
        for (int i = 0; i < n; i++)
        {
            var x = long.Parse(ReadLine());
            WriteLine(x.ToString() + " " + (MillerRabin(x) ? "1" : "0"));
        }
    }

    public static int[] b64 = new int[] { 2, 325, 9375, 28178, 450775, 9780504, 1795265022 };// n < at least 2^64
    public static bool MillerRabin(long n)
    {
        if (n < 2) return false;
        if ((n & 1) == 0) return n == 2;
        var d = n - 1;
        var s = 0;
        while ((d & 1) == 0)
        {
            s++;
            d >>= 1;
        }
        for (int i = 0; i < b64.Length && b64[i] < n; i++)
        {
            var c = true;//このへんメソッド抽出すべきか したらアーリーリターンできるか
            var x = BigInteger.ModPow(b64[i], d, n);
            if (x == 1) c = false;
            for (int r = 0; r < s; r++)
            {
                if (x == n - 1) c = false;
                x = x * x % n;
            }
            if (c) return false;
        }
        return true;
    }
}
0