結果

問題 No.260 世界のなんとか3
ユーザー E31415926
提出日時 2017-02-21 21:57:30
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,653 bytes
コンパイル時間 1,972 ms
コンパイル使用メモリ 171,640 KB
実行使用メモリ 11,136 KB
最終ジャッジ日時 2024-12-30 15:13:23
合計ジャッジ時間 4,062 ms
ジャッジサーバーID
(参考情報)
judge2 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 2 WA * 24 RE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #

#include<bits/stdc++.h>
#define int long long
using namespace std;
#define rep(i,n) for(int i=0;i<(n);++i)
#define INF (1ll<<60)
typedef pair<int, int> pii;
#define FI first
#define SE second
#define all(s) s.begin(),s.end()
#define RREP(i,n) for(int i=(n)-1;i>=0;--i)

int A;
string B;
const int mod = 1000000007;
int dp[10001][2][2][3][8]; //桁,未満確定(0未定1確定),3つく(0つかない1つく),mod3,mod8

signed main() {
	cin >> A >> B;
	dp[0][0][0][0][0] = 1;
	rep(i,B.length())
	{
		rep(j,2)
		{
			rep(k,2)
			{
				rep(l,3)
				{
					rep(m,8)
					{
						int lim = j ? 9ll : B[i] - '0';
						rep(d,lim+1)
						{
							(dp[i + 1][j || d < lim][k || d == 3][(l + d) % 3][(m
									* 10 + d) % 8] += dp[i][j][k][l][m]) %= mod;
						}
					}
				}
			}
		}
	}
	int ansB = 0;
	rep(j,2)
	{
		rep(l,3)
		{
			rep(m,7)
			{
				ansB += dp[B.length()][j][1][l][m + 1] % mod;
			}
		}
	}
	rep(j,2)
	{
		rep(m,7)
		{
			ansB += dp[B.length()][j][0][0][m + 1] % mod;
		}
	}
	A--;
	string AA = to_string(A);
	memset(dp, 0, sizeof(dp));
	dp[0][0][0][0][0] = 1;
	rep(i,AA.length())
	{
		rep(j,2)
		{
			rep(k,2)
			{
				rep(l,3)
				{
					rep(m,8)
					{
						int lim = j ? 9ll : AA[i] - '0';
						rep(d,lim+1)
						{
							(dp[i + 1][j || d < lim][k || d == 3][(l + d) % 3][(m
									* 10 + d) % 8] += dp[i][j][k][l][m]) %= mod;
						}
					}
				}
			}
		}
	}
	int ansA = 0;
	rep(j,2)
	{
		rep(l,3)
		{
			rep(m,7)
			{
				ansA += dp[AA.length()][j][1][l][m + 1] % mod;
			}
		}
	}
	rep(j,2)
	{
		rep(m,7)
		{
			ansA += dp[AA.length()][j][0][0][m + 1] % mod;
		}
	}
	int ans = ((ansB + mod) - ansA) % mod;
	cout << ans << endl;
}
0