結果
問題 | No.528 10^9と10^9+7と回文 |
ユーザー |
![]() |
提出日時 | 2021-04-02 16:48:53 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 325 ms / 1,000 ms |
コード長 | 3,849 bytes |
コンパイル時間 | 3,092 ms |
コンパイル使用メモリ | 202,748 KB |
最終ジャッジ日時 | 2025-01-20 07:56:04 |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 28 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define ALL(x) begin(x),end(x)#define rep(i,n) for(int i=0;i<(n);i++)#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;using ll=long long;const int INF=1000000000;const ll LINF=1001002003004005006ll;int dx[]={1,0,-1,0},dy[]={0,1,0,-1};// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}struct IOSetup{IOSetup(){cin.tie(0);ios::sync_with_stdio(0);cout<<fixed<<setprecision(12);}} iosetup;template<typename T>ostream &operator<<(ostream &os,const vector<T>&v){for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");return os;}template<typename T>istream &operator>>(istream &is,vector<T>&v){for(T &x:v)is>>x;return is;}// a^n (mod m)ll pow_mod(ll a,ll n,ll m){ll ret=1;while(n){if(n&1) ret=ret*a%m;a=(a*a)%m;n/=2;}return ret;}ll inv_mod(ll a,ll m){ll b=m,u=1,v=0,t;while(b){t=a/b;swap(a-=t*b,b);swap(u-=t*v,v);}u%=m;if(u<0) u+=m;return u;}ll solve(string s,ll mod){int n=(int)s.size();int m=n/2;// 0:up, 1:eq, 2:lwauto left=[&](int flag,int put,int dig){if(flag==0){return flag;}if(flag==1){if(put==dig) return 1;else if(put<dig) return 2;return 0;}if(flag==2){return 2;}return -1;};auto right=[&](int flag,int put,int dig){if(dig==put){return flag;}if(put<dig) return 2;return 0;};vector<vector<vector<ll>>> dp(m+1,vector<vector<ll>>(3,vector<ll>(3,0)));dp[0][1][1]=1;for(int i=0;i<m;i++){int u=s[i]-'0',d=s[n-1-i]-'0';for(int p=0;p<10;p++){if(i==0 and p==0) continue;rep(j,3){rep(k,3){int lt=left(j,p,u);int rt=right(k,p,d);dp[i+1][lt][rt]=(dp[i+1][lt][rt]+dp[i][j][k])%mod;}}}}// rep(j,3)rep(k,3){// cout<<j<<", "<<k<<" : "<<dp[m][j][k]<<endl;// }ll ret=0;if(n&1){int mid=s[n/2]-'0';rep(i,10){if(i<mid){// ret+= dp[m][0][0]+dp[m][0][0]+dp[m][0][0]+// dp[m][1][0]+dp[m][1][1]+dp[m][1][2]+// dp[m][2][0]+dp[m][2][1]+dp[m][2][2];ret+=dp[m][1][0]+dp[m][1][1]+dp[m][1][2]+dp[m][2][0]+dp[m][2][1]+dp[m][2][2];ret%=mod;}else if(i==mid){ret+=dp[m][1][1]+dp[m][1][2]+dp[m][2][0]+dp[m][2][1]+dp[m][2][2];ret%=mod;}else{ret+=dp[m][2][0]+dp[m][2][1]+dp[m][2][2];ret%=mod;}}}else{ret+=dp[m][1][1]+dp[m][1][2]+dp[m][2][0]+dp[m][2][1]+dp[m][2][2];ret%=mod;}return ret;}void unk(string s,ll mod){int n=(int)s.size();if(n==1){cout<<(s[0]-'0')<<endl;return ;}ll res=0;for(int i=1;i<n;i++){int fr=i/2+i%2-1;ll add=pow_mod(10,fr,mod)*9%mod;res=(res+add)%mod;}res+=solve(s,mod);res%=mod;cout<<res<<endl;}signed main(){string s;cin>>s;unk(s,1000000000);unk(s,1000000007);return 0;}