結果
問題 | 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 1000000000LLint 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;}