結果
問題 | No.315 世界のなんとか3.5 |
ユーザー |
![]() |
提出日時 | 2023-04-22 17:35:57 |
言語 | Go (1.23.4) |
結果 |
AC
|
実行時間 | 1,762 ms / 2,000 ms |
コード長 | 5,751 bytes |
コンパイル時間 | 13,454 ms |
コンパイル使用メモリ | 224,172 KB |
実行使用メモリ | 8,960 KB |
最終ジャッジ日時 | 2024-11-07 05:46:32 |
合計ジャッジ時間 | 38,380 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 36 |
ソースコード
package mainimport ("bufio""fmt""io""os""strconv")var iost *Iosttype Iost struct {Scanner *bufio.ScannerWriter *bufio.Writer}func NewIost(fp io.Reader, wfp io.Writer) *Iost {const BufSize = 2000005scanner := bufio.NewScanner(fp)scanner.Split(bufio.ScanWords)scanner.Buffer(make([]byte, BufSize), BufSize)return &Iost{Scanner: scanner, Writer: bufio.NewWriter(wfp)}}func (i *Iost) Text() string {if !i.Scanner.Scan() {panic("scan failed")}return i.Scanner.Text()}func (i *Iost) Atoi(s string) int { x, _ := strconv.Atoi(s); return x }func (i *Iost) GetNextInt() int { return i.Atoi(i.Text()) }func (i *Iost) Atoi64(s string) int64 { x, _ := strconv.ParseInt(s, 10, 64); return x }func (i *Iost) GetNextInt64() int64 { return i.Atoi64(i.Text()) }func (i *Iost) Atof64(s string) float64 { x, _ := strconv.ParseFloat(s, 64); return x }func (i *Iost) GetNextFloat64() float64 { return i.Atof64(i.Text()) }func (i *Iost) Print(x ...interface{}) { fmt.Fprint(i.Writer, x...) }func (i *Iost) Printf(s string, x ...interface{}) { fmt.Fprintf(i.Writer, s, x...) }func (i *Iost) Println(x ...interface{}) { fmt.Fprintln(i.Writer, x...) }func isLocal() bool { return os.Getenv("NICKEL") == "BACK" }func main() {fp := os.Stdinwfp := os.Stdoutif isLocal() {fp, _ = os.Open(os.Getenv("WELL_EVERYBODY_LIES_TOO_MUCH"))}iost = NewIost(fp, wfp)defer func() {iost.Writer.Flush()}()solve()}func solve() {SetMod(Mod1000000007)a := iost.Text()b := iost.Text()p := iost.GetNextInt()iost.Println(g(b, 0).Sub(g(a, 1)).Sub(f(b, 0, p)).Add(f(a, 1, p)))}func f(s string, sub, p int) Mint {n := len(s)ss := make([]int, n)for i := 0; i < n; i++ {ss[i] = int(s[i] - '0')}ss[n-1] -= subfor i := n - 2; i >= 0; i-- {if ss[i+1] < 0 {ss[i+1] += 10ss[i]--}}if ss[0] == 0 {ss = ss[1:]}for p%10 == 0 && len(ss) > 0 {p /= 10ss = ss[:len(ss)-1]}n = len(ss)if n == 0 {return 0}dp := [2][2][24]Mint{}dp[0][0][0] = 1var ans Mintm := 1for i := n - 1; i >= 0; i-- {dp2 := [2][2][24]Mint{}for j := 0; j < 2; j++ {for k := 0; k < 2; k++ {for l := 0; l < 24; l++ {a := dp[j][k][l] % Mod1000000007tl := l - mfor d := 0; d < 10; d++ {tl += mif tl >= 24 {tl -= 24}tj := 0if ss[i]-j < d {tj = 1}tk := kif d == 3 {tk = 1}dp2[tj][tk][tl] += aif d == 0 {continue}if tk == 0 && tl%3 != 0 {continue}if tl%8 != 0 {continue}if i > 0 || tj == 0 {ans += a}}}}}dp, dp2 = dp2, dpm = m * 10 % 24}return ans % Mod1000000007}func g(s string, sub int) Mint {n := len(s)ss := make([]int, n)for i := 0; i < n; i++ {ss[i] = int(s[i] - '0')}ss[n-1] -= subfor i := n - 2; i >= 0; i-- {if ss[i+1] < 0 {ss[i+1] += 10ss[i]--}}if ss[0] == 0 {ss = ss[1:]}n = len(ss)if n == 0 {return 0}dp := [2][2][3]Mint{}dp[0][0][0] = 1var ans Mintm := 1for i := n - 1; i >= 0; i-- {dp2 := [2][2][3]Mint{}for j := 0; j < 2; j++ {for k := 0; k < 2; k++ {for l := 0; l < 3; l++ {a := dp[j][k][l] % Mod1000000007tl := l - mfor d := 0; d < 10; d++ {tl += mif tl >= 3 {tl -= 3}tj := 0if ss[i]-j < d {tj = 1}tk := kif d == 3 {tk = 1}tl := (l + d*m) % 3dp2[tj][tk][tl] += aif d == 0 {continue}if tk == 0 && tl%3 != 0 {continue}if i > 0 || tj == 0 {ans += a}}}}}dp, dp2 = dp2, dpm = m * 10 % 3}return ans % Mod1000000007}// Mod constants.const (Mod1000000007 = 1000000007Mod998244353 = 998244353)var (mod Mintfmod func(Mint) Mint)// Mint treats the modular arithmetictype Mint int64// SetMod sets the mod. It must be called first.func SetMod(newmod Mint) {switch newmod {case Mod1000000007:fmod = staticMod1000000007case Mod998244353:fmod = staticMod998244353default:mod = newmodfmod = dynamicMod}}func dynamicMod(m Mint) Mint {m %= modif m < 0 {return m + mod}return m}func staticMod1000000007(m Mint) Mint {m %= Mod1000000007if m < 0 {return m + Mod1000000007}return m}func staticMod998244353(m Mint) Mint {m %= Mod998244353if m < 0 {return m + Mod998244353}return m}// Mod returns m % mod.func (m Mint) Mod() Mint {return fmod(m)}// Inv returns modular multiplicative inversefunc (m Mint) Inv() Mint {return m.Pow(Mint(0).Sub(2))}// Pow returns m^nfunc (m Mint) Pow(n Mint) Mint {p := Mint(1)for n > 0 {if n&1 == 1 {p.MulAs(m)}m.MulAs(m)n >>= 1}return p}// Add returns m+xfunc (m Mint) Add(x Mint) Mint {return (m + x).Mod()}// Sub returns m-xfunc (m Mint) Sub(x Mint) Mint {return (m - x).Mod()}// Mul returns m*xfunc (m Mint) Mul(x Mint) Mint {return (m * x).Mod()}// Div returns m/xfunc (m Mint) Div(x Mint) Mint {return m.Mul(x.Inv())}// AddAs assigns *m + x to *m and returns mfunc (m *Mint) AddAs(x Mint) *Mint {*m = m.Add(x)return m}// SubAs assigns *m - x to *m and returns mfunc (m *Mint) SubAs(x Mint) *Mint {*m = m.Sub(x)return m}// MulAs assigns *m * x to *m and returns mfunc (m *Mint) MulAs(x Mint) *Mint {*m = m.Mul(x)return m}// DivAs assigns *m / x to *m and returns mfunc (m *Mint) DivAs(x Mint) *Mint {*m = m.Div(x)return m}