結果
| 問題 |
No.36 素数が嫌い!
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2015-10-05 22:52:40 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,492 bytes |
| コンパイル時間 | 402 ms |
| コンパイル使用メモリ | 55,800 KB |
| 実行使用メモリ | 5,376 KB |
| 最終ジャッジ日時 | 2024-07-20 01:12:58 |
| 合計ジャッジ時間 | 1,288 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | WA * 4 |
| other | AC * 10 WA * 16 |
ソースコード
#include <iostream>
using namespace std;
class Math
{
private:
public:
static int calcPowMod(int x, int n, int mod);
static bool isPrime(unsigned long n);
static int countPrime(int max);
};
/**
* calcPowMod
* x^n % mod の 計算結果を返します
* (二分累乗法)
**/
int Math::calcPowMod(int x, int n, int mod)
{
unsigned long sum;
if ( n == 0 ) { return 1; }
if ( n == 1 ) { return x % mod; }
sum = Math::calcPowMod(x, n/2, mod);
return n%2 == 0 ? (sum * sum) % mod : (x * sum * sum) % mod;
}
/**
* isPrime
* 素数かどうか判定します
**/
bool Math::isPrime(unsigned long n)
{
int k;
if ( n == 1 ) { return false; }
// if ( n < 2 ) { return true; }
else if ( n == 2 ) { return true; }
if ( n%2 == 0 ) { return true; }
for ( k = 3; k <= n/k; k+=2 ) { if ( n%k == 0 ) { return true; } }
return false;
}
/**
* countPrime
* maxまでの素数の数を数えます
**/
int Math::countPrime(int max)
{
int i, j, k;
int count = 1;
int *prime;
max = (max - 3) / 2;
prime = new int(max);
if ( prime == NULL ) { return -1; }
for ( k = 0; k < max; k++ ) { prime[k] = 1; }
for ( i = 0; i < max; i++ ) {
int prime_num = i+i+3;
if ( prime[i] == 0 ) { continue; }
count++;
for ( j = i+prime_num; j < max; j += prime_num ) { prime[j] = 0; }
}
delete[] prime;
return count;
}
int main(void)
{
unsigned long max;
cin >> max;
string ans = ! Math::isPrime(max) ? "YES" : "NO";
cout << ans << endl;
}