結果

問題 No.1532 Different Products
コンテスト
ユーザー 0w1
提出日時 2021-06-28 01:31:01
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 983 ms / 4,000 ms
コード長 3,397 bytes
コンパイル時間 2,000 ms
コンパイル使用メモリ 197,344 KB
最終ジャッジ日時 2025-01-22 14:45:16
ジャッジサーバーID
（参考情報） 		judge3 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 62
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

/*^ debug  */

template <typename A, typename B> string to_string(pair<A, B> p);
template <typename A, typename B, typename C> string to_string(tuple<A, B, C> p);
template <typename A, typename B, typename C, typename D> string to_string(tuple<A, B, C, D> p);
string to_string(const string& s) { return '"' + s + '"'; }
string to_string(const char* s) { return to_string((string) s); }
string to_string(bool b) { return (b ? "true" : "false"); }
string to_string(vector<bool> v) {
  bool first = true;
  string res = "{";
  for (int i = 0; i < static_cast<int>(v.size()); i++) {
    if (!first) { res += ", "; }
    first = false;
    res += to_string(v[i]);
  }
  res += "}";
  return res;
}
template <size_t N>
string to_string(bitset<N> v) {
  string res = "";
  for (size_t i = 0; i < N; i++) { res += static_cast<char>('0' + v[i]); }
  return res;
}
template <typename A>
string to_string(A v) {
  bool first = true;
  string res = "{";
  for (const auto &x : v) {
    if (!first) { res += ", "; }
    first = false;
    res += to_string(x);
  }
  res += "}";
  return res;
}
template <typename A, typename B>
string to_string(pair<A, B> p) { return "(" + to_string(p.first) + ", " + to_string(p.second) + ")"; }
template <typename A, typename B, typename C>
string to_string(tuple<A, B, C> p) { return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ")"; }
template <typename A, typename B, typename C, typename D>
string to_string(tuple<A, B, C, D> p) { return "(" + to_string(get<0>(p)) + ", " + to_string(get<1>(p)) + ", " + to_string(get<2>(p)) + ", " + to_string(get<3>(p)) + ")"; }
void debug_out() { cerr << endl; }
template <typename Head, typename... Tail>
void debug_out(Head H, Tail... T) { cerr << " " << to_string(H); debug_out(T...); }
#ifdef LOCAL
#define debug(...) cerr << "[" << #__VA_ARGS__ << "]:", debug_out(__VA_ARGS__)
#else
#define debug(...) 42
#endif

/*  debug $*/

/*^ generic definitions  */

template<typename F>
struct _RecurFun : F {
  _RecurFun(F&& f) : F(forward<F>(f)) {}
  template<typename... Args>
  decltype(auto) operator()(Args&&... args) const { return F::operator()(*this, forward<Args>(args)...); }
};
template<typename F>
decltype(auto) RecurFun(F&& f) { return _RecurFun<F> { forward<F>(f) }; }

template<typename D, typename R>
vector<R> maps(const vector<D> &&d, function<R(D)> f) {
  vector<R> r(d.size());
  for (int i = 0; i < (int) d.size(); ++i) r[i] = f(d[i]);
  return r;
}

/*  generic definitions $*/

int main() {
  ios::sync_with_stdio(false);

  int N;
  int64_t K; {
    cin >> N >> K;
  }

  if (K == 1) {
    cout << 1 << endl;
    return 0;
  }

  const int L = 1<<17;

  vector<uint64_t> small(L+1);
  vector<uint64_t> large(L+1); {
    for (int i = 1; i <= N; ++i) {
      vector<uint64_t> psmall(L+1); {
        partial_sum(small.begin(), small.end(), psmall.begin());
      }
      for (int j = 1; j <= L; ++j) {
        int64_t x = i*j;
        if (x <= L) {
          large[j] += (K/j>=i) + large[x];
        } else {
          assert(K/x <= L);
          large[j] += (K/j>=i) + psmall[K/x];
        }
      }
      for (int j = L; j; --j) {
        if (j%i == 0) small[j] += small[j/i];
      }
      small[i] += 1;
      /* debug(i); */
      /* debug(small); */
      /* debug(large); */
    }
  }

  cout << large[1] << endl;
}
0