結果

問題 No.1161 Many Powers
ユーザー tails
提出日時 2020-11-18 20:01:00
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 12 ms / 2,000 ms
コード長 882 bytes
コンパイル時間 237 ms
コンパイル使用メモリ 31,488 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-23 09:09:48
合計ジャッジ時間 942 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 23
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:8:1: warning: return type defaults to 'int' [-Wimplicit-int]
    8 | main(){
      | ^~~~
main.c: In function 'main':
main.c:10:9: warning: implicit declaration of function 'read' [-Wimplicit-function-declaration]
   10 |         read(0,rbuf,sizeof rbuf);
      |         ^~~~
main.c:67:17: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   67 |                 write(1,&d,n);
      |                 ^~~~~
main.c:69:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   69 |         _exit(0);
      |         ^~~~~
      |         _Exit

ソースコード

diff # 

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

#define RDL(v) long v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}

int w[100001];

main(){
	char rbuf[64];
	read(0,rbuf,sizeof rbuf);
	char*rp=rbuf;
	RDL(a);
	RDL(b);
	RDL(c);
	{
		long p=c/(c&-c);
		long e=p;
		long i=3;
		while(e>1){
			if(e%i==0){
				p-=p/i;
				while(e%i==0){
					e/=i;
				}
			}
			i+=2;
			if(i*i>e){
				i=e;
			}
		}
		if(b>=p){
			b=b%p+p;
		}
	}
	long z=0,y=0;
	long ac=a%c;
	for(long i=1;i<=c;++i){
		long x;
		if((i&1)==0&&i!=2){
			x=(long)w[2]*w[i/2]%c;
		}else if((int)i%3==0&&i!=3){
			x=(long)w[3]*w[i/3]%c;
		}else{
			x=1;
			long e=i;
			long d=b;
			while(d){
				if(d&1){
					x=x*e%c;
				}
				d>>=1;
				e=e*e%c;
			}
		}
		w[i]=x;
		y+=x;
		if(i==ac){
			z=y;
		}
	}
	z+=a/c*y;
	z=z%c;

	{
		long n=0,d=0;
		while(++n,d=d<<8|0x30|z%10,z/=10);
		write(1,&d,n);
	}
	_exit(0);
}
0