結果

問題 No.2012 Largest Triangle
ユーザー tails
提出日時 2022-08-03 21:08:20
言語 C90
(gcc 12.3.0)
結果
AC  
実行時間 12 ms / 2,500 ms
コード長 1,626 bytes
コンパイル時間 448 ms
コンパイル使用メモリ 26,880 KB
実行使用メモリ 7,180 KB
最終ジャッジ日時 2024-09-13 14:50:39
合計ジャッジ時間 2,664 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 41
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function ‘f1’:
main.c:59:20: warning: passing argument 1 of ‘radix_sort’ from incompatible pointer type [-Wincompatible-pointer-types]
   59 |         radix_sort(a,n);
      |                    ^
      |                    |
      |                    P *
main.c:29:23: note: expected ‘ulong *’ {aka ‘long unsigned int *’} but argument is of type ‘P *’
   29 | void radix_sort(ulong*a,int n){
      |                 ~~~~~~^
main.c: In function ‘f4’:
main.c:97:9: warning: implicit declaration of function ‘write’ [-Wimplicit-function-declaration]
   97 |         write(1,wp,wbuf+sizeof wbuf-wp);
      |         ^~~~~
main.c: In function ‘main’:
main.c:102:9: warning: implicit declaration of function ‘_exit’ [-Wimplicit-function-declaration]
  102 |         _exit(0);
      |         ^~~~~
main.c:102:9: warning: incompatible implicit declaration of built-in function ‘_exit’ [-Wbuiltin-declaration-mismatch]

ソースコード

diff # 

#pragma GCC optimize("Ofast")

#define rd() ({long _v=0,_c;while(_c=*rp++-48,_c>=0)_v=_v*10+_c;_v;})
#define wt(v) {ulong _z=v;do*--wp=_z%10+48;while(_z/=10);}
#define rep(v,e) for(long v=0;v<e;++v)

typedef unsigned long ulong;
typedef struct { int x,y; } P;

void radix_sort_aux(ulong*a,ulong*b,int n){
	int c[2048];
	for(int i=0;i<2048;++i){
		c[i]=0;
	}
	for(int i=0;i<n;++i){
		++c[a[i]&2047];
	}
	int t=0;
	for(int i=2048;i--;){
		int u=c[i];
		c[i]=t;
		t+=u;
	}
	for(int i=0;i<n;++i){
		b[c[a[i]&2047]++]=a[i]>>11|a[i]<<53;
	}
}

void radix_sort(ulong*a,int n){
	static ulong b[200000];
	radix_sort_aux(a,b,n);
	radix_sort_aux(b,a,n);
	for(int i=0;i<n;++i){
		a[i]=a[i]>>42|a[i]<<22;
	}
}

P a[200000];

long f1(){
	char*mmap();
	char*rp=mmap(0l,1l<<25,1,2,0,0ll);
	long n=rd();
	rep(i,n){
		long neg=0;
		if(*rp=='-'){
			++rp;
			neg=!neg;
		}
		long x=rd();
		if(*rp=='-'){
			++rp;
			neg=!neg;
		}
		long y=rd();
		a[i].y=y;
		a[i].x=(neg?-x:x)+(1<<20);
	}
	radix_sort(a,n);
	return n;
}

long f2(long n){
	long m=0;
	rep(i,n){
		long y=a[i].y;
		long x=a[i].x-(1<<20);
		while(m>=2&&(a[m-1].x-x)*(a[m-2].y-y)-(a[m-1].y-y)*(a[m-2].x-x)>=0){
			--m;
		}
		a[m].y=y;
		a[m].x=x;
		++m;
	}
	return m;
}

long f3(long m){
	long z=0;
	long r=0;
	rep(i,m){
		r=r<i?i:r;
		long u=(long)a[i].x*(long)a[r].y-(long)a[i].y*(long)a[r].x;
		for(long v;v=(long)a[i].x*(long)a[r+1].y-(long)a[i].y*(long)a[r+1].x,u<v;){
			u=v;
			++r;
			r=r==m?0:r;
		}
		z=z<u?u:z;
	}
	return z;
}

void f4(ulong z){
	char wbuf[64],*wp=wbuf+sizeof wbuf;
	wt(z);
	write(1,wp,wbuf+sizeof wbuf-wp);
}

int main(){
	f4(f3(f2(f1())));
	_exit(0);
}
0