結果
| 問題 | No.732 3PrimeCounting |
| コンテスト | |
| ユーザー |
dama_math
|
| 提出日時 | 2018-09-07 23:08:31 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
RE
|
| 実行時間 | - |
| コード長 | 997 bytes |
| コンパイル時間 | 1,644 ms |
| コンパイル使用メモリ | 168,864 KB |
| 実行使用メモリ | 5,248 KB |
| 最終ジャッジ日時 | 2024-11-29 19:28:37 |
| 合計ジャッジ時間 | 15,405 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 68 WA * 4 RE * 17 |
ソースコード
#include <bits/stdc++.h>
#define ll long long
#define PB push_back
#define FOR(i,a,b) for(int i=a;i<b;i++)
using namespace std;
//素数判定
bool pj(int n){
if(n==1) return false;
int d=2;
while(d*d<=n){
if(n%d==0)return false;
d++;
}
return true;
}
int m[300000];
vector<int> v;
int main(){
int n; cin>>n;
memset(m,0,sizeof(m));
FOR(i,1,3*n+1) if(pj(i)) v.PB(i);
for(int i=0;v[i]<=n;i++){
for(int j=0;v[j]<=n;j++){
m[v[i]+v[j]]++;
}
}
ll ans1=0,ans2=0;
for(int i=0;v[i]<=3*n;i++){
for(int j=0;v[j]<=n;j++){
if(v[i]>v[j]){
ans1+=m[v[i]-v[j]];
}
}
}
memset(m,0,sizeof(m));
for(int i=0;v[i]<=n;i++){
m[2*v[i]]++;
}
for(int i=0;v[i]<=3*n;i++){
for(int j=0;v[j]<=n;j++){
if(v[i]>v[j]){
ans2+=m[v[i]-v[j]];
}
}
}
cout<<(ans1-ans2*3)/6<<endl;
return 0;
}
dama_math