結果

問題 No.1884 Sequence
コンテスト
ユーザー 👑 tails
提出日時 2022-03-26 22:14:18
言語 C11(gcc12 gnu拡張)
(gcc 12.4.0)
コンパイル:
gcc-12 -O2 -std=gnu11 -DONLINE_JUDGE -o a.out _filename_ -lm
実行:
./a.out
結果
AC  
実行時間 20 ms / 2,000 ms
コード長 1,453 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 102 ms
コンパイル使用メモリ 33,480 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2026-03-08 16:30:54
合計ジャッジ時間 2,375 ms
ジャッジサーバーID
（参考情報） 		judge3 / 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’ [-Wimplicit-function-declaration]
   60 |         _exit(0);
      |         ^~~~~
main.c:60:9: warning: incompatible implicit declaration of built-in function ‘_exit’ [-Wbuiltin-declaration-mismatch]
main.c: In function ‘no’:
main.c:66:9: warning: incompatible implicit declaration of built-in function ‘_exit’ [-Wbuiltin-declaration-mismatch]
   66 |         _exit(0);
      |         ^~~~~

ソースコード

diff #
raw source code 

#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