結果

問題 No.2230 Good Omen of White Lotus
ユーザー tailstails
提出日時 2023-02-25 17:06:17
言語 C90
(gcc 11.4.0)
結果
AC  
実行時間 22 ms / 2,000 ms
コード長 1,763 bytes
コンパイル時間 565 ms
コンパイル使用メモリ 29,644 KB
実行使用メモリ 8,136 KB
最終ジャッジ日時 2023-10-11 14:35:25
合計ジャッジ時間 3,710 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,348 KB
testcase_01 AC 0 ms
4,348 KB
testcase_02 AC 1 ms
4,352 KB
testcase_03 AC 1 ms
4,352 KB
testcase_04 AC 1 ms
4,348 KB
testcase_05 AC 1 ms
4,352 KB
testcase_06 AC 0 ms
4,352 KB
testcase_07 AC 1 ms
4,348 KB
testcase_08 AC 1 ms
4,348 KB
testcase_09 AC 1 ms
4,356 KB
testcase_10 AC 0 ms
4,352 KB
testcase_11 AC 1 ms
4,352 KB
testcase_12 AC 1 ms
4,348 KB
testcase_13 AC 1 ms
4,348 KB
testcase_14 AC 8 ms
4,380 KB
testcase_15 AC 1 ms
4,352 KB
testcase_16 AC 14 ms
6,128 KB
testcase_17 AC 14 ms
6,192 KB
testcase_18 AC 14 ms
6,084 KB
testcase_19 AC 14 ms
6,136 KB
testcase_20 AC 2 ms
4,352 KB
testcase_21 AC 1 ms
4,352 KB
testcase_22 AC 1 ms
4,352 KB
testcase_23 AC 2 ms
4,352 KB
testcase_24 AC 1 ms
4,352 KB
testcase_25 AC 1 ms
4,352 KB
testcase_26 AC 1 ms
4,348 KB
testcase_27 AC 15 ms
6,952 KB
testcase_28 AC 15 ms
7,920 KB
testcase_29 AC 14 ms
7,012 KB
testcase_30 AC 15 ms
8,136 KB
testcase_31 AC 18 ms
7,892 KB
testcase_32 AC 15 ms
7,104 KB
testcase_33 AC 22 ms
7,100 KB
testcase_34 AC 22 ms
7,136 KB
testcase_35 AC 22 ms
7,044 KB
testcase_36 AC 22 ms
7,104 KB
testcase_37 AC 21 ms
7,040 KB
testcase_38 AC 22 ms
7,140 KB
testcase_39 AC 22 ms
7,044 KB
testcase_40 AC 6 ms
4,348 KB
testcase_41 AC 6 ms
4,352 KB
testcase_42 AC 14 ms
5,164 KB
testcase_43 AC 1 ms
4,348 KB
testcase_44 AC 19 ms
6,316 KB
testcase_45 AC 7 ms
4,348 KB
testcase_46 AC 11 ms
4,444 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function ‘main’:
main.c:112:7: warning: passing argument 1 of ‘sort’ from incompatible pointer type [-Wincompatible-pointer-types]
  sort(xy,n);
       ^~
main.c:66:17: note: expected ‘ulong *’ {aka ‘long unsigned int *’} but argument is of type ‘XY *’ {aka ‘struct <anonymous> *’}
 void sort(ulong*a,int n){
           ~~~~~~^
main.c:7:59: warning: implicit declaration of function ‘write’ [-Wimplicit-function-declaration]
 #define wt1(v) ({char wbuf[64],*wp=wbuf+sizeof wbuf;wt(v);write(1,wp,wbuf+sizeof wbuf-wp);})
                                                           ^~~~~
main.c:116:2: note: in expansion of macro ‘wt1’
  wt1(y);
  ^~~
main.c:117:2: warning: implicit declaration of function ‘_exit’ [-Wimplicit-function-declaration]
  _exit(0);
  ^~~~~
main.c:117:2: warning: incompatible implicit declaration of built-in function ‘_exit’

ソースコード

diff # 

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

#define getrp() ({char*mmap();mmap(0l,1l<<25,1,2,0,0ll);})
#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 wt1(v) ({char wbuf[64],*wp=wbuf+sizeof wbuf;wt(v);write(1,wp,wbuf+sizeof wbuf-wp);})
#define rep(v,e) for(long v=0;v<e;++v)
#define MD 998244353

int modpow(int a, unsigned b){
	int r=1;
	while(b){
		if(b&1){
			r=1ll*r*a%MD;
		}
		b>>=1;
		a=1ll*a*a%MD;
	}
	return r;
}

int inverse(int a){
	int b=MD;
	int u=1;
	int v=0;
	while(b){
		int q=a/b,t;
		t=b, b=a-q*b, a=t;
		t=v, v=u-q*v, u=t;
	}
	if(u<0){
		u+=MD;
	}
	return u;
}

typedef unsigned long ulong;

void sort_aux(ulong*a,ulong*b,int n,int ph){
	int c[512];
	for(int i=0;i<512;++i){
		c[i]=0;
	}
	for(int i=0;i<n;++i){
		++c[a[i]&511];
	}
	int t=0;
	for(int i=0;i<512;++i){
		int u=c[i];
		c[i]=t;
		t+=u;
	}
	if(ph==0){
		for(int i=0;i<n;++i){
			b[c[a[i]&511]++]=a[i]>>9|a[i]<<55;
		}
	}
	if(ph==1){
		for(int i=0;i<n;++i){
			b[c[a[i]&511]++]=a[i]>>23|a[i]<<41;
		}
	}
}

void sort(ulong*a,int n){
	ulong*b=a+n;
	sort_aux(a,b,n,0);
	sort_aux(b,a,n,1);
	sort_aux(a,b,n,0);
	sort_aux(b,a,n,1);
}

typedef struct{
	unsigned x,y;
}XY;

XY xy[200000*2];

int lis[200000];

long f1(long n){
	long z=0;
	rep(i,n){
		int x=xy[i].x;
		long l=-1,r=z;
		while(l+1<r){
			long m=l+r>>1;
			if(lis[m]<=x){
				l=m;
			}else{
				r=m;
			}
		}
		lis[r]=x;
		if(r==z){
			++z;
		}
	}
	return z;
}

int main(){
	char*rp=getrp();
	long h=rd()+rd()-3;
	long n=rd();
	long p=rd();
	rep(i,n){
		xy[i].x=rd();
		xy[i].y=rd();
	}
	sort(xy,n);
	long z=f1(n);
	long x=1+MD-inverse(p);
	long y=(1+MD-1ll*modpow(x,h-z)*modpow((x*2-1)%MD,z)%MD)%MD;
	wt1(y);
	_exit(0);
}
0