結果

問題 No.1322 Totient Bound
ユーザー chocoruskchocorusk
提出日時 2020-12-13 20:03:46
言語 Java21
(openjdk 21)
結果
AC  
実行時間 1,436 ms / 5,000 ms
コード長 2,456 bytes
コンパイル時間 2,541 ms
コンパイル使用メモリ 77,428 KB
実行使用メモリ 95,456 KB
最終ジャッジ日時 2023-10-20 04:06:57
合計ジャッジ時間 29,334 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 176 ms
82,412 KB
testcase_01 AC 175 ms
82,456 KB
testcase_02 AC 173 ms
82,480 KB
testcase_03 AC 174 ms
82,440 KB
testcase_04 AC 170 ms
82,448 KB
testcase_05 AC 170 ms
82,400 KB
testcase_06 AC 174 ms
82,432 KB
testcase_07 AC 171 ms
82,408 KB
testcase_08 AC 171 ms
82,436 KB
testcase_09 AC 171 ms
82,436 KB
testcase_10 AC 171 ms
82,456 KB
testcase_11 AC 173 ms
82,420 KB
testcase_12 AC 186 ms
82,448 KB
testcase_13 AC 188 ms
82,548 KB
testcase_14 AC 187 ms
82,396 KB
testcase_15 AC 181 ms
82,400 KB
testcase_16 AC 185 ms
82,472 KB
testcase_17 AC 181 ms
82,432 KB
testcase_18 AC 429 ms
83,104 KB
testcase_19 AC 569 ms
84,988 KB
testcase_20 AC 1,133 ms
84,964 KB
testcase_21 AC 1,131 ms
85,096 KB
testcase_22 AC 1,151 ms
85,132 KB
testcase_23 AC 1,324 ms
95,208 KB
testcase_24 AC 1,327 ms
95,444 KB
testcase_25 AC 1,285 ms
95,456 KB
testcase_26 AC 1,307 ms
95,448 KB
testcase_27 AC 1,294 ms
95,216 KB
testcase_28 AC 1,312 ms
95,368 KB
testcase_29 AC 1,310 ms
95,344 KB
testcase_30 AC 173 ms
82,432 KB
testcase_31 AC 172 ms
82,452 KB
testcase_32 AC 188 ms
82,456 KB
testcase_33 AC 1,425 ms
95,280 KB
testcase_34 AC 1,436 ms
95,236 KB
testcase_35 AC 1,348 ms
95,316 KB
testcase_36 AC 1,431 ms
95,332 KB
testcase_37 AC 1,435 ms
95,344 KB
testcase_38 AC 1,323 ms
95,368 KB
権限があれば一括ダウンロードができます

ソースコード

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 m-(int)(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<=isq; 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)isq; 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