結果

問題 No.1514 Squared Matching
ユーザー tails
提出日時 2021-05-22 11:15:24
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 65 ms / 4,000 ms
コード長 774 bytes
コンパイル時間 147 ms
コンパイル使用メモリ 31,104 KB
実行使用メモリ 7,536 KB
最終ジャッジ日時 2024-10-10 10:52:06
合計ジャッジ時間 2,171 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:11:1: warning: return type defaults to 'int' [-Wimplicit-int]
   11 | main(){
      | ^~~~
main.c: In function 'main':
main.c:41:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   41 |         write(1,wp,wbuf+sizeof wbuf-wp);
      |         ^~~~~
main.c:42:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   42 |         _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);}

unsigned long vis[50000000/64+4];

main(){
	char*rp=mmap(0l,1l<<25,1,2,0,0ll);
	rd(n);
	long res=0;

	long n4=n/4+1;
	for(long i=1;i<n4;++i){
		if(!(vis[i>>6]&1ul<<i)){
			long count=1;
			for(int j=2,ijj;ijj=i*j*j,ijj<=n;++j){
				vis[ijj>>6]|=1ul<<ijj;
				count+=1;
			}
			res+=count*count;
		}
	}

	if(n4&63){
		vis[n4>>6]|=~0ul>>-n4;
	}
	vis[n+1>>6]|=~0ul<<n+1;

	long e=(n>>6)+1;
	res+=e-(n4>>6)<<6;
	for(long i=n4>>6;i<e;++i){
		res-=__builtin_popcountl(vis[i]);
	}

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