結果

問題 No.732 3PrimeCounting
ユーザー 37zigen37zigen
提出日時 2020-04-04 16:30:39
言語 Java
(openjdk 23)
結果
AC  
実行時間 2,404 ms / 3,000 ms
コード長 1,495 bytes
コンパイル時間 2,347 ms
コンパイル使用メモリ 79,352 KB
実行使用メモリ 45,484 KB
最終ジャッジ日時 2024-07-03 07:31:29
合計ジャッジ時間 51,742 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 89
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.util.ArrayList;
import java.util.Arrays;
import java.util.Scanner;

public class Main {

	public static void main(String[] args) {
		long curtime=System.currentTimeMillis();
		new Main().run();
		System.err.println(System.currentTimeMillis()-curtime);
	}
	
	int MAX=300001;
	boolean[] isPrime=new boolean[MAX];
	ArrayList<Integer> primes=new ArrayList<>();
	{
		Arrays.fill(isPrime, true);
		isPrime[0]=isPrime[1]=false;
		for(int i=2;i<isPrime.length;++i) {
			if(!isPrime[i])continue;
			primes.add(i);
			for(long j=(long)i*i;j<isPrime.length;j+=i)
				isPrime[(int)j]=false;
		}
	}
	
	void run() {
		Scanner sc=new Scanner(System.in);
		int N=sc.nextInt();//5<=N<=1e5
		long[] dp=new long[MAX];
		for(int i=0;i<primes.size();++i) {
			int p=primes.get(i);
			if(p>N)break;
			++dp[p+p];
			for(int j=i+1;j<primes.size();++j) {
				int q=primes.get(j);
				if(q>N)break;
				dp[p+q]+=2;
			}
		}
		
		long ans=0;
		int src=0;
		for(int p:primes) {
			if(p>N)break;
			while(src<primes.size()&&primes.get(src)-p<=0)++src;
			for(int i=src;i<primes.size();++i) {
				int target=primes.get(i);
				ans+=dp[target-p];
			}
		}
		//a+b+c が　素数　かつ　a==b==c　は存在しない。
		for(int p:primes) {
			if(p>N)break;
			for(int q:primes) {
				if(q>N)break;
				if(isPrime[p+p+q]) {
					ans-=3;
				}
			}
		}
		if(ans%6!=0)throw new AssertionError();
		System.out.println(ans/6);
	}
	
	void tr(Object...objects) {System.out.println(Arrays.deepToString(objects));}
}
0