結果

問題 No.1514 Squared Matching
コンテスト
ユーザー tails
提出日時 2021-05-21 23:08:39
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 3,952 ms / 4,000 ms
コード長 844 bytes
コンパイル時間 196 ms
コンパイル使用メモリ 31,360 KB
実行使用メモリ 489,760 KB
最終ジャッジ日時 2024-10-10 10:03:13
合計ジャッジ時間 82,888 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:26:1: warning: return type defaults to 'int' [-Wimplicit-int]
   26 | main(){
      | ^~~~
main.c: In function 'main':
main.c:30:9: warning: implicit declaration of function 'read'; did you mean 'rd'? [-Wimplicit-function-declaration]
   30 |         read(0,rbuf,sizeof rbuf);
      |         ^~~~
      |         rd
main.c:54:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   54 |         write(1,wp,wbuf+sizeof wbuf-wp);
      |         ^~~~~
main.c:55:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   55 |         _exit(0);
      |         ^~~~~
      |         _Exit

ソースコード

diff # 

#pragma GCC optimize("Ofast")
#pragma GCC target("avx2")

char*mmap();
#define rd(v) int v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}

#define wt_rev(v) {long _z=v;do*--wp=_z%10+48;while(_z/=10);}

#define SIEVE_N 50000001
int sieve[SIEVE_N];
void mksieve(){
	for(int i=2;i<SIEVE_N;++i){
		if(sieve[i]==0){
			for(int j=i;j<SIEVE_N;j+=i){
				if(sieve[j]==0){
					sieve[j]=i;
				}
			}
		}
	}
}

int cac[50000001];
short count[50000001];

main(){
	mksieve();

	char rbuf[64],*rp=rbuf;
	read(0,rbuf,sizeof rbuf);
	rd(n);
	
	int res=0;
	for(int i=1;i<=n;++i){
		int j=i;
		int l=1;
		while(j>1){
			int p=sieve[j];
			j/=p;
			if(sieve[j]==p){
				l*=cac[j/p];
				break;
			}
			l*=p;
		}

		cac[i]=l;
		count[l]+=1;
		res+=count[l]*2-1;
	}

	char wbuf[64],*wp=wbuf+sizeof wbuf;
	wt_rev(res);
	write(1,wp,wbuf+sizeof wbuf-wp);
	_exit(0);
}
0