結果

問題 No.1322 Totient Bound
ユーザー chocoruskchocorusk
提出日時 2020-12-13 19:58:35
言語 Java21
(openjdk 21)
結果
RE  
実行時間 -
コード長 2,454 bytes
コンパイル時間 2,320 ms
コンパイル使用メモリ 77,132 KB
実行使用メモリ 94,136 KB
最終ジャッジ日時 2024-09-19 23:51:17
合計ジャッジ時間 16,970 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 167 ms
66,016 KB
testcase_01 AC 166 ms
66,020 KB
testcase_02 AC 167 ms
65,956 KB
testcase_03 AC 163 ms
66,232 KB
testcase_04 AC 168 ms
65,920 KB
testcase_05 AC 166 ms
66,016 KB
testcase_06 AC 164 ms
66,208 KB
testcase_07 AC 168 ms
66,084 KB
testcase_08 AC 172 ms
66,112 KB
testcase_09 AC 169 ms
66,108 KB
testcase_10 AC 169 ms
66,080 KB
testcase_11 AC 169 ms
65,924 KB
testcase_12 AC 167 ms
66,004 KB
testcase_13 AC 194 ms
66,364 KB
testcase_14 AC 183 ms
66,132 KB
testcase_15 AC 173 ms
65,928 KB
testcase_16 AC 177 ms
66,116 KB
testcase_17 AC 172 ms
66,028 KB
testcase_18 AC 413 ms
67,524 KB
testcase_19 AC 543 ms
67,664 KB
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 RE -
testcase_25 RE -
testcase_26 RE -
testcase_27 RE -
testcase_28 RE -
testcase_29 RE -
testcase_30 AC 164 ms
66,124 KB
testcase_31 AC 167 ms
66,004 KB
testcase_32 AC 165 ms
65,956 KB
testcase_33 RE -
testcase_34 RE -
testcase_35 RE -
testcase_36 RE -
testcase_37 RE -
testcase_38 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import java.util.Scanner;

public class Main {
	public static final int MAX=1000000;
	public static boolean[] isprime=new boolean[MAX];
	public static void sieve() {
		for(int i=2; i<MAX; i++) isprime[i]=true;
		for(int i=2; i<MAX; i++) {
			if(isprime[i]) {
				for(int j=(i<<1); j<MAX; j+=i) {
					isprime[j]=false;
				}
			}
		}
	}
	public static long n, sq;
	public static int isq;
	public static long[] p=new long[MAX];
	public static long[] divs=new long[MAX];
	public static long[] pc=new long[MAX];
	public static int m;
	public static void calc() {
		while((sq+1)*(sq+1)<=n) sq++;
		for(long i=1; i<=sq; i++) divs[m++]=i;
		for(long i=sq; i>=1; i--) if(n/i>sq) divs[m++]=n/i;
		int k=0;
		isq=-1;
		for(int i=2; i<MAX; i++) {
			if(isprime[i]) {
				p[k]=i;
				if(p[k]>sq && isq==-1) isq=k;
				k++;
			}
		}
	}
	public static int idx(long x) {
		if(x<=sq) return (int)(x-1);
		else return (int)(m-n/x);
	}
	public static void primecount(){
		long[] dp=new long[m];
		for(int i=0; i<m; i++) dp[i]=divs[i];
		int l=0, r=0;
		while(r<m && divs[r]<p[0]*p[0]) r++;
		for(int i=l; i<r; i++) {
			if(divs[i]==1) pc[i]=0;
			else if(divs[i]==2) pc[i]=1;
			else pc[i]=2;
		}
		for(int i=1; i<=sq; i++) {
			l=r;
			while(r<m && divs[r]<p[i]*p[i]) r++;
			for(int j=m-1; j>=l; j--) {
				long a=dp[j];
				if(divs[j]<p[i-1]*p[i-1]) a=pc[j]-(i-1)+1;
				int k=idx(divs[j]/p[i-1]);
				long b=dp[k];
				if(divs[k]<p[i-1]*p[i-1]) b=pc[k]-(i-1)+1;
				dp[j]=a-b;
			}
			for(int j=l; j<r; j++) {
				pc[j]=dp[j]-1+i;
			}
		}
	}
	public static boolean is_prime(long x) {
		if(x<MAX) return isprime[(int)x];
		for(long i=2; i*i<=x; i++) {
			if(x%i==0) return false;
		}
		return true;
	}
	public static void main(String[] args) {
		Scanner scanner=new Scanner(System.in);
		n=scanner.nextLong();
		sieve();
		calc();
		primecount();
		for(int i=0; i<m; i++) {
			if(is_prime(divs[i]+1)) pc[i]++;
		}
		long[] dp=new long[m];
		int l=m-1;
		for(int i=(int)sq; i>=0; i--) {
			while(l>=0 && divs[l]>=p[i]*(p[i]-1)) l--;
			for(int j=m-1; j>=l+1; j--) {
				long a=dp[j];
				if(divs[j]<p[i+1]-1) a=1;
				else if(divs[j]<p[i+1]*(p[i+1]-1)) a=pc[j]-(i+1)+1;
				long q=p[i]-1;
				while(q<=divs[j]) {
					long x=divs[j]/q;
					int k=idx(x);
					long b=dp[k];
					if(x<p[i+1]-1) b=1;
					else if(x<p[i+1]*(p[i+1]-1)) b=pc[k]-(i+1)+1;
					a+=b;
					q*=p[i];
				}
				dp[j]=a;
			}
		}
		dp[0]=2;
		System.out.println(dp[m-1]);
		scanner.close();
	}
}
0