結果
問題 | No.1086 桁和の桁和2 |
ユーザー | 沙耶花 |
提出日時 | 2020-06-19 22:39:55 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 385 ms / 3,000 ms |
コード長 | 1,752 bytes |
コンパイル時間 | 2,777 ms |
コンパイル使用メモリ | 204,860 KB |
最終ジャッジ日時 | 2025-01-11 07:13:35 |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 5 |
other | AC * 31 |
ソースコード
#include <bits/stdc++.h> using namespace std; #define modulo 1000000007 #define mod(mod_x) ((((long long)mod_x+modulo))%modulo) #define Inf 1000000000LL int beki(int a,long long b,int M = modulo){ int x = 1; while(b!=0){ if(b&1){ x=((long long)x*a)%M; } a=((long long)a*a)%M; b>>=1; } return x; } int gyakugen(int a){ return beki(a,modulo-2); } vector<int> Merge2(vector<int> A,vector<int> B){ vector<int> ret(9,0); for(int i=0;i<A.size();i++){ for(int j=0;j<B.size();j++){ int nxt = (i+j)%9; int D = mod(A[i]*B[j]); ret[nxt] = mod(ret[nxt] + D); } } return ret; } vector<int> calc(long long X){ vector<int> ret(9,0); if(X==0)return ret; int t = beki(10,X); t = mod(t - 1); t = mod(t * gyakugen(9)); for(int i=0;i<9;i++)ret[i] = t; return ret; } vector<int> calc(long long L,long long R){ vector<int> ret = calc(R); vector<int> t = calc(L); for(int i=0;i<ret.size();i++)ret[i] = mod(ret[i] - t[i]); ret[0] = mod(ret[0] + 1); return ret; } int main(){ int N; cin>>N; deque<long long> L(N),R(N); deque<int> D(N); for(int i=0;i<N;i++)cin>>L[i]; for(int i=0;i<N;i++)cin>>R[i]; for(int i=0;i<N;i++)cin>>D[i]; while(D.size()>0&&D.front()==0){ D.pop_front(); L.pop_front(); R.pop_front(); } if(D.size()==0){ cout<<1<<endl; return 0; } N = D.size(); bool zero = false; for(int i=0;i<N;i++){ if(D[i]==0)zero=true; } if(zero){ cout<<0<<endl; return 0; } vector<int> now = {1,0,0,0,0,0,0,0,0}; for(int i=0;i<N;i++){ now = Merge2(now,calc(L[i],R[i])); for(int j=0;j<9;j++){ if(j!=D[i]%9)now[j] = 0; } if(D[i]==9&&i==0){ now[0] = mod(now[0] - 1); } } int ans= 0; for(int i=0;i<now.size();i++)ans =mod(ans + now[i]); cout<<ans<<endl; return 0; }