結果

問題 No.1884 Sequence
ユーザー tails
提出日時 2022-03-26 22:14:18
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 22 ms / 2,000 ms
コード長 1,453 bytes
コンパイル時間 415 ms
コンパイル使用メモリ 31,104 KB
実行使用メモリ 7,680 KB
最終ジャッジ日時 2024-10-15 14:13:12
合計ジャッジ時間 2,647 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 40
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function 'yes':
main.c:59:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   59 |         write(1,"Yes",3);
      |         ^~~~~
main.c:60:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   60 |         _exit(0);
      |         ^~~~~
      |         _Exit
main.c:61:1: warning: 'noreturn' function does return
   61 | }
      | ^
main.c: In function 'no':
main.c:67:1: warning: 'noreturn' function does return
   67 | }
      | ^

ソースコード

diff # 

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

#define rd(v) long v=0;{long _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}
#define rep(v,e) for(long v=0;v<e;++v)

typedef unsigned long ulong;

void radix_sort_aux(ulong*a,ulong*b,int n){
	int c[256];
	for(int i=0;i<256;++i){
		c[i]=0;
	}
	for(int i=0;i<n;++i){
		++c[a[i]&255];
	}
	int t=0;
	for(int i=0;i<256;++i){
		int u=c[i];
		c[i]=t;
		t+=u;
	}
	for(int i=0;i<n;++i){
		b[c[a[i]&255]++]=a[i]>>8|a[i]<<56;
	}
}

void radix_sort(ulong*a,int n){
	static ulong b[200000];
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
}

long gcd(long a,long b){
	long x=__builtin_ctzl(a);
	long y=__builtin_ctzl(b);
	a>>=x;
	b>>=y;
	long c;
	while(c=a-b){
		c>>=__builtin_ctzl(c);
		if(c<0){
			b=-c;
		}else{
			a=c;
		}
	}
	return a<<(x<=y?x:y);
}

void yes() __attribute__((noreturn));
void yes(){
	write(1,"Yes",3);
	_exit(0);
}

void no() __attribute__((noreturn));
void no(){
	write(1,"No",2);
	_exit(0);
}

ulong a[200000];

int main(){
	char*mmap();
	char*rp=mmap(0l,1l<<25,1,2,0,0ll);
	rd(n);
	long l=0;
	while(*rp){
		rd(x);
		if(x){
			a[l++]=x;
		}
	}
	if(l<2){
		yes();
	}
	radix_sort(a,l);
	long g=a[l-1]-a[0];
	if(!g){
		yes();
	}
	rep(i,l-1){
		long d=a[i+1]-a[i];
		if(!d){
			no();
		}
		g=gcd(g,d);
	}
	if((n-1)*g<a[l-1]-a[0]){
		no();
	}
	yes();
}
0