結果

問題 No.1973 Divisor Sequence
ユーザー renren_utsrenren_uts
提出日時 2022-06-16 07:21:35
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 727 ms / 2,000 ms
コード長 1,613 bytes
コンパイル時間 1,067 ms
コンパイル使用メモリ 112,376 KB
実行使用メモリ 179,700 KB
最終ジャッジ日時 2024-10-05 18:49:25
合計ジャッジ時間 7,156 ms
ジャッジサーバーID
(参考情報)
judge4 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 147 ms
6,816 KB
testcase_01 AC 2 ms
6,816 KB
testcase_02 AC 186 ms
30,552 KB
testcase_03 AC 158 ms
6,820 KB
testcase_04 AC 171 ms
19,320 KB
testcase_05 AC 165 ms
17,224 KB
testcase_06 AC 181 ms
25,828 KB
testcase_07 AC 154 ms
7,496 KB
testcase_08 AC 163 ms
15,424 KB
testcase_09 AC 165 ms
17,416 KB
testcase_10 AC 195 ms
20,812 KB
testcase_11 AC 150 ms
6,816 KB
testcase_12 AC 178 ms
21,488 KB
testcase_13 AC 172 ms
8,780 KB
testcase_14 AC 189 ms
18,264 KB
testcase_15 AC 171 ms
22,332 KB
testcase_16 AC 180 ms
26,640 KB
testcase_17 AC 185 ms
18,092 KB
testcase_18 AC 148 ms
6,820 KB
testcase_19 AC 202 ms
22,496 KB
testcase_20 AC 158 ms
6,864 KB
testcase_21 AC 187 ms
30,948 KB
testcase_22 AC 173 ms
22,508 KB
testcase_23 AC 727 ms
179,700 KB
testcase_24 AC 331 ms
46,768 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:67:12: warning: structured bindings only available with '-std=c++17' or '-std=gnu++17' [-Wc++17-extensions]
   67 |   for(auto [i,vaku]:mp){
      |            ^

ソースコード

diff #

#include <iostream>
#include <stdio.h>
#include <cmath>
#include <vector>
#include <tuple>
#include <bitset>
#include <map>
#include <queue>
#include <time.h>
#include <utility>
#include <numeric>
#include <deque>
#include <cassert>
#include <sstream>
#include <cctype>
#include <unordered_map>
#include <algorithm>
#include <chrono>
#define ll long long
#include <time.h>
#define mod % 1000000007
#define rep(i,n) for(long long i = 0;i < n;i ++)
ll MAX = 5'000'000'000'000'000'000LL;
using namespace std;

 vector<ll> p;

void MAKE_P(){
  p = {2,3};
  for(int i = 5;i <= 1000000;i ++){
    bool z = true;
    for(ll j:p){
      if(i%j==0){
        z = false;
        break;
      }
      if(j*j>i)break;
    }
    if(z)p.push_back(i);
  }
}

int main(){
  ll n,m;
  cin >> n >> m;
  if(m==1){
    cout << 1 << endl;
    return 0;
  }
  MAKE_P();
  vector<ll> q;
  rep(i,int(p.size())){
    if(m%p[i]==0){
      q.push_back(0);
      while(m%p[i] == 0){
        q[int(q.size())-1] ++;
        m = m/p[i];
      }
    }
  }
  if(m>1)q.push_back(1);
  map<ll,ll> mp;
  for(int i:q){
    mp[i] = 1;
  }
  vector<ll> ans(50);
  for(auto [i,vaku]:mp){
    if(i==0){
      ans[0] = 1;
      continue;
    }
    vector<vector<ll>> dp(n+1,vector<ll>(i+1,0));
    dp[0][0] = 1;
    rep(j,n){
      rep(k,i+1){
        dp[j][k] = dp[j][k] mod;
        rep(l,i+1){
          if(l+k>i)break;
          dp[j+1][l] += dp[j][k];
        }
      }
    }
    rep(j,i+1)ans[i] += dp[n][j];
    ans[i] = ans[i] mod;
  }
  ll sum = 1;
  rep(i,int(q.size())){
    sum = (sum*ans[q[i]]) mod;
  }
  cout << sum << endl;
  
  return 0;
}
0