結果

問題 No.528 10^9と10^9+7と回文
ユーザー yukudo
提出日時 2017-06-10 00:24:28
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 95 ms / 1,000 ms
コード長 2,260 bytes
コンパイル時間 1,478 ms
コンパイル使用メモリ 167,080 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-23 18:30:55
合計ジャッジ時間 3,380 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
#define REP(i,n) for(int i=0,_n=(int)(n);i<_n;++i)
template<class T>bool chmin(T&a,const T&b){return a>b?(a=b,true):false;}
template<class T>bool chmax(T&a,const T&b){return a<b?(a=b,true):false;}

ll mod_pow(ll a, ll b, ll p) {
  ll res = 1;
  while (b > 0) {
    if (b & 1) res = (res * a) % p;
    a = (a * a) % p;
    b >>= 1;
  }
  return res;
}

ll f(int N, ll MOD) {
  return 9 * mod_pow(10, (N-1)/2, MOD) % MOD;
}

const int MAX_N = 100000 + 10;
ll dp[MAX_N][2];

ll solve(vector<int> a, ll MOD) {
  memset(dp, 0, sizeof(dp));
  int N = (int)a.size();
  dp[0][0] = 1;
  for (int i = 0; i < (N+1)/2; i++) {
    for (int k = 0; k < 2; k++) {
      for (int d = (i==0?1:0); d <= 9; d++) {
        if (k == 1 || d < a[i]) {
          (dp[i+1][1] += dp[i][k]) %= MOD;
        } else if (k == 0 && a[i] == d) {
          (dp[i+1][0] += dp[i][k]) %= MOD;
        }
      }
    }
  }
  // REP(i, N+1) {
    // cout << dp[i][0] << " " << dp[i][1] << endl;
  // }
  ll res = 0;
  res += dp[(N+1)/2][1];
  bool is = true;
  vector<int> b;
  REP(i, N/2) b.push_back(a[i]);
  if (N % 2) b.push_back(a[N/2]);
  REP(i, N/2) b.push_back(a[N/2-1-i]);
  REP(i, N) {
    if (a[i] > b[i]) break;
    if (a[i] < b[i]) is = false;

  }
  if (is) {
    res = (res + 1) % MOD;
  }

  // 桁数がN未満の回文数
  for (int i = 1; i < N; i++) {
    res += f(i, MOD) % MOD;
    res %= MOD;
  }

  return res % MOD;
}


bool isok(int x) {
  stringstream ss;
  ss << x;
  string s = ss.str();
  for (int i = 0; i < (int)s.size()/2; i++) {
    if (s[i] != s[s.size() - 1 - i])
      return false;
  }
  return true;
}

ll solveSmall(const vector<int> a, ll MOD) {
  int N = (int)a.size();
  if (N >= 6) return -1;
  int n = 0;
  REP(i, N) {
    n = n * 10 + a[i];
  }
  ll res = 0;
  for (int x = 1; x <= n; x++) {
    if (isok(x)) ++res;
  }
  return res;
}

int main2() {
  string s; cin >> s;
  vector<int> a(s.size());
  REP(i, s.size()) a[i] = s[i] - '0';
  cout << solve(a, (ll)1e9) << endl;
  cout << solve(a, (ll)1e9+7) << endl;
  // cout << solveSmall(a, (ll)1e9) << endl;
  // cout << solveSmall(a, (ll)1e9+7) << endl;
  return 0;
}

int main() {
  for (;!cin.eof();cin>>ws) main2();
  return 0;
}
0