結果

問題 No.260 世界のなんとか3
コンテスト
ユーザー bal4u
提出日時 2019-08-30 20:41:45
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 13 ms / 2,000 ms
コード長 2,288 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 312 ms
コンパイル使用メモリ 34,116 KB
実行使用メモリ 5,632 KB
最終ジャッジ日時 2024-11-21 20:09:24
合計ジャッジ時間 1,621 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.c: In function 'ins':
main.c:10:14: warning: implicit declaration of function 'getchar_unlocked' [-Wimplicit-function-declaration]
   10 | #define gc() getchar_unlocked()
      |              ^~~~~~~~~~~~~~~~
main.c:17:16: note: in expansion of macro 'gc'
   17 |         do c = gc(), *s++ = c & 0xf;
      |                ^~

ソースコード

diff #
raw source code 

// yukicoder: 260 世界のなんとか3
// 2019.8.30 bal4u

#include <stdio.h>
#include <string.h>

typedef long long ll;

#if 1
#define gc() getchar_unlocked()
#else
#define gc() getchar()
#endif

int ins(char *s) { // 文字列の入力　スペース以下の文字で入力終了
	char *p = s; int c;
	do c = gc(), *s++ = c & 0xf;
	while (c > ' ');
	*--s = 0;
	return s - p;
}

#define MOD 1000000007

char AA[10005], *A; int WA;
char B[10005]; int WB;
int dp[10005][2][3][8][2];   // 3である桁, mod 3, mod 8, 元の数
int C[4] = {0,1,2,4};

inline static void chmod(int *a, int b) { *a = ((ll)(*a) + b) % MOD; }

int calc(int w, char *s) {
	int i, j, i3, m3, m8, k;
	
	memset(dp, 0, sizeof(dp));
	dp[0][0][0][0][1] = 1;
    for (i = 0; i < w; i++) for (i3 = 0; i3 < 2; i3++) for (m3 = 0; m3 < 3; m3++) for (m8 = 0; m8 < 8; m8++)
		for (k = 0; k < 2; k++)	if (dp[i][i3][m3][m8][k]) {
			if (k) {
               	for (j = s[i]-1; j >= 0; j--) {
               		if (w-i <= 3) chmod(&dp[i+1][i3|(j==3)][(m3+j)%3][(m8+C[w-i]*j)&7][0], dp[i][i3][m3][m8][k]);
               		else          chmod(&dp[i+1][i3|(j==3)][(m3+j)%3][             m8][0], dp[i][i3][m3][m8][k]);
               	}
           		if (w-i <= 3) chmod(&dp[i+1][i3|(s[i]==3)][(m3+s[i])%3][(m8+C[w-i]*s[i])&7][k], dp[i][i3][m3][m8][k]);
           		else 		  chmod(&dp[i+1][i3|(s[i]==3)][(m3+s[i])%3][                m8][k], dp[i][i3][m3][m8][k]);
           	} else {
               	for (j = 9; j >= 0; j--) {
              		if (w-i <= 3) chmod(&dp[i+1][i3|(j==3)][(m3+j)%3][(m8+C[w-i]*j)&7][k], dp[i][i3][m3][m8][k]);
                    else          chmod(&dp[i+1][i3|(j==3)][(m3+j)%3][             m8][k], dp[i][i3][m3][m8][k]);
                }
            }
        }
    k = 0;
	for (m8 = 1; m8 < 8; m8++) chmod(&k, dp[w][0][0][m8][0] + dp[w][0][0][m8][1]);
	for (m3 = 0; m3 < 3; m3++) for (m8 = 1; m8 < 8; m8++) chmod(&k, dp[w][1][m3][m8][0] + dp[w][1][m3][m8][1]);
	return k;
}

int main()
{
	int i, ans;
	
	WA = ins(A = AA), WB = ins(B);
	if (WA == 1 && *A == 1) ans = calc(WB, B);
	else {
		for (i = WA-1; i >= 0; i--) {
			if (A[i] > 0) { A[i]--; break; }
			A[i] = 9;
		}
		if (A[0] == 0) WA--, A++;
		ans = calc(WB, B) - calc(WA, A);
	}
	if (ans < 0) ans += MOD;
	printf("%d\n", ans);
	return 0;
}
0