結果

問題 No.473 和と積の和
ユーザー char134217728char134217728
提出日時 2018-01-01 22:33:07
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 253 ms / 3,000 ms
コード長 4,022 bytes
コンパイル時間 2,969 ms
コンパイル使用メモリ 217,192 KB
実行使用メモリ 225,344 KB
最終ジャッジ日時 2023-08-24 15:47:55
合計ジャッジ時間 5,982 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
7,520 KB
testcase_01 AC 4 ms
9,708 KB
testcase_02 AC 4 ms
7,612 KB
testcase_03 AC 148 ms
91,160 KB
testcase_04 AC 5 ms
9,800 KB
testcase_05 AC 4 ms
7,524 KB
testcase_06 AC 4 ms
4,380 KB
testcase_07 AC 4 ms
7,484 KB
testcase_08 AC 5 ms
7,760 KB
testcase_09 AC 5 ms
11,768 KB
testcase_10 AC 4 ms
4,384 KB
testcase_11 AC 4 ms
4,380 KB
testcase_12 AC 4 ms
7,688 KB
testcase_13 AC 5 ms
11,772 KB
testcase_14 AC 5 ms
13,888 KB
testcase_15 AC 6 ms
15,828 KB
testcase_16 AC 6 ms
15,872 KB
testcase_17 AC 6 ms
17,916 KB
testcase_18 AC 4 ms
7,740 KB
testcase_19 AC 4 ms
9,600 KB
testcase_20 AC 4 ms
4,380 KB
testcase_21 AC 4 ms
4,380 KB
testcase_22 AC 3 ms
4,380 KB
testcase_23 AC 8 ms
24,604 KB
testcase_24 AC 11 ms
38,984 KB
testcase_25 AC 253 ms
225,344 KB
testcase_26 AC 104 ms
86,400 KB
testcase_27 AC 104 ms
110,920 KB
testcase_28 AC 3 ms
4,384 KB
testcase_29 AC 8 ms
28,248 KB
testcase_30 AC 8 ms
26,100 KB
testcase_31 AC 7 ms
24,056 KB
testcase_32 AC 7 ms
21,968 KB
testcase_33 AC 6 ms
13,988 KB
testcase_34 AC 125 ms
164,180 KB
testcase_35 AC 181 ms
142,676 KB
testcase_36 AC 16 ms
63,736 KB
testcase_37 AC 6 ms
15,816 KB
testcase_38 AC 221 ms
131,208 KB
testcase_39 AC 132 ms
90,884 KB
testcase_40 AC 124 ms
96,540 KB
testcase_41 AC 4 ms
4,384 KB
testcase_42 AC 20 ms
67,136 KB
testcase_43 AC 34 ms
84,976 KB
testcase_44 AC 131 ms
102,892 KB
testcase_45 AC 5 ms
11,712 KB
testcase_46 AC 4 ms
7,752 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp:155:1: 警告: ISO C++ では型の無い ‘main’ の宣言を禁止しています [-Wreturn-type]
  155 | main(){
      | ^~~~

ソースコード

diff # 

#include <bits/stdc++.h>
#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define FORR(i,a,b) for (int i=(a);i>=(b);i--)
#define pb push_back
#define pcnt __builtin_popcount
#define show(x) cout<<#x<<" = "<<x<<endl;
#define maxs(x,y) x = max(x,y)
#define mins(x,y) x = min(x,y)
#define fi first
#define se second
#define rng(a) a.begin(),a.end()
#define each(it,c) for(__typeof((c).begin()) it=(c).begin();it!=(c).end();it++)
#define sz(x) (int)(x).size()
#define mp make_pair

using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<pii> vpii;
typedef set<int> si;
typedef pair<ll,ll> pll;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<pll> vpll;
typedef set<ll> sl;
typedef __int128_t lll;
typedef pair<lll,lll> plll;
typedef vector<lll> vlll;
template<typename T>string join(vector<T>&v)
{stringstream s;FOR(i,0,sz(v))s<<' '<<v[i];return s.str().substr(1);}
ll gcd(ll a,ll b){if(a>b)swap(a,b);for(;a>0;b%=a,swap(a,b));return b;}
int modpow(ll a,ll n,int m){if(a==0)return a;ll p=1;for(;n>0;n/=2,a=a*a%m)if(n&1)p=p*a%m;return(int)p;}
void dout(double d){printf("%.12f\n",d);}

const int iinf = 1e9;
const ll linf = 1e18;
const int mod = 1e9+7;
const double pi = acos(-1);
const double eps = 1e-10;

struct Prime{
  vl pl;
  random_device rnd;
  mt19937_64 mt;
  Prime(){
    mt = mt19937_64(rnd());
  }
  __int128_t mp(ll _a,ll n,ll m){
    if(_a==0)return 0;
    __int128_t a=_a, p=1;
    for(;n>0;n/=2,a=a*a%m)if(n&1)p=p*a%m;
    return p;
  }
  bool miller_rabin(ll _n){
    __int128_t n, t, y;
    ll a, d;
    _n = abs(_n);
    if(_n==2) return true;
    if(_n==1||(_n&1)==0) return false;
    uniform_int_distribution<ll> dist(1, _n-1);
    n = __int128_t(_n);
    d = _n-1;
    while((d&1)==0) d>>=1;
    FOR(i, 0, 100){
      a = dist(mt);
      t = d;
      y = mp(a,t,_n);
      while(t != _n-1 && y != 1 &&y != n-1){
        y = (y*y)%_n;
        t <<= 1;
      }
      if(y != _n-1 && (t&1) == 0) return false;
    }
    return true;
  }
  bool is_prime(ll a){
    if(a < 2)return false;
    each(itr, pl){
      if(*itr * *itr > a)return true;
      if(a%*itr == 0)return false;
    }
    return true;
  }
  void init(int a){
    pl.resize(0);
    pl.pb(2);
    pl.pb(3);
    int i = 5;
    while(i <= a){
      if(is_prime(i)) pl.pb(i);
      i += 2;
      if(is_prime(i)) pl.pb(i);
      i += 4;
    }
  }
  void prime_division(ll a, vpii&v){ // vに素因数分解の結果が入る
    v.resize(0);
    for(vl::iterator itr = pl.begin(); *itr * *itr <= a; itr++){
      if(a%*itr == 0){
        ll c = 0;
        while((a%*itr) == 0){
          a /= *itr;
          c++;
        }
        v.pb(pii(*itr, c));
      }
    }
    if(a > 1)v.pb(pii(a, 1));
  }
  void divisor(ll a, vl&v){ // v に約数が入る
    v.pb(1);
    for(vl::iterator itr = pl.begin(); *itr * *itr <= a; itr++){
      if(a%*itr == 0){
        ll b = 1, vs = sz(v);
        while((a%*itr) == 0){
          a /= *itr;
          b *= *itr;
          FOR(j, 0, vs)v.pb(v[j] * b);
        }
      }
    }
    if(a>1){
      int vs = sz(v);
      FOR(j, 0, vs)v.pb(v[j] * a);
    }
  }
};
ll memo[35][2000][2000];
map<ll, int> m;
vl v;
ll N, X;
ll get(int n, ll a, ll b){
  if(memo[n][a][b] == -1){
    if(n == 1){
      memo[n][a][b] = a >= b ? 1 : 0;
    }else{
      ll su = 0, _a = v[a], _b = v[b];
      FOR(i, 0, n) _a /= _b;
      if(_a == 0){
        memo[n][a][b] = 0;
        return 0;
      }
      _a = v[a];
      FOR(i, b, sz(v)){
        if(_a % v[i] != 0) continue;
        su += get(n-1, m[_a/v[i]], i);
      }
      memo[n][a][b] = su;
    }
  }
  return memo[n][a][b];
}
main(){
  cin.tie(0);
  ios::sync_with_stdio(false);
  Prime pr;
  pr.init(34000);
  cin >> N >> X;
  if(N >= 35 || X < 3){
    cout << 0 << endl;
    return 0;
  }
  pr.divisor(X+1, v);
  sort(rng(v));
  FOR(i, 0, sz(v)) m[v[i]] = i;
  FOR(i, 0, N+1)FOR(j, 0, sz(v))FOR(k, 0, sz(v)) memo[i][j][k] = -1;
  cout << get(N, m[X+1], 1) << endl;
  return 0;
}
0