結果

問題 No.1581 Multiple Sequence
ユーザー tails
提出日時 2021-07-06 10:49:06
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 4 ms / 2,000 ms
コード長 1,130 bytes
コンパイル時間 324 ms
コンパイル使用メモリ 37,504 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-07-01 12:06:16
合計ジャッジ時間 1,166 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 21
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c:23:1: warning: return type defaults to 'int' [-Wimplicit-int]
   23 | main(){
      | ^~~~
main.c: In function 'main':
main.c:26:9: warning: implicit declaration of function 'read'; did you mean 'rd'? [-Wimplicit-function-declaration]
   26 |         read(0,rp,sizeof buf);
      |         ^~~~
      |         rd
main.c:29:16: warning: initialization of 'long int *' from incompatible pointer type 'int *' [-Wincompatible-pointer-types]
   29 |         long*p=dp+1;
      |                ^~
main.c:52:9: warning: implicit declaration of function 'write' [-Wimplicit-function-declaration]
   52 |         write(1,wp,buf+sizeof buf-wp);
      |         ^~~~~
main.c:53:9: warning: implicit declaration of function '_exit'; did you mean '_Exit'? [-Wimplicit-function-declaration]
   53 |         _exit(0);
      |         ^~~~~
      |         _Exit

ソースコード

diff # 

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

char*mmap();
#define rd(v) long v=0;{int _c;while(_c=*rp++-48,_c>=0)v=v*10+_c;}
#define wt(v) {long _z=v;do*--wp=_z%10+48;while(_z/=10);}

#define MD 1000000007

#define LIM 100000

#define A1  1, 2, 3, 5, 6,10, 1, 5, 3, 7, 1,13,
#define A2  1, 2, 8, 5, 1,10, 1,10, 3, 2, 1,13,
#define A3  6, 2, 3, 5, 1,15, 1, 5, 3, 2, 6,13,
#define A4  1, 2, 3,10, 1,10, 1, 5, 8, 2, 1,13,
#define A5  1, 7, 3, 5, 1,10, 6, 5, 3, 2, 1,18,
#define A A1 A2 A3 A4 A5
#define B A A A A A A A A A A
#define C B B B B B B B B B B
#define D C C C C C C C C C C
int dp[]={0,D};

main(){
	char buf[64];
	char*rp=buf;
	read(0,rp,sizeof buf);
	rd(m);

	long*p=dp+1;
	long r=0;
	r+=1*(m/1-LIM/2/1);
	r+=1*(m/2-LIM/2/2);
	r+=2*(m/3-LIM/2/3);
	r+=3*(m/4-LIM/2/4);
	r+=5*(m/5-LIM/2/5);
	r+=6*(m/6-LIM/2/6);
	for(long i=7;i<=LIM/2;++i){
		long a=dp[i-1];
		for(long j=i;j<=LIM/2;j+=i){
			long b=a+dp[j];
			dp[j]=b+(b>=MD?-MD:0);
		}
		r+=a*((m+0)/i-(LIM/2+0)/i);
	}
	if(m<=LIM/2){
		r=dp[m];
	}else{
		r=(r+dp[LIM/2])%MD;
	}
	char*wp=buf+sizeof buf;
	wt(r);
	write(1,wp,buf+sizeof buf-wp);
	_exit(0);
}
0