結果
| 問題 |
No.473 和と積の和
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-05-17 11:58:20 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 1,025 ms / 3,000 ms |
| コード長 | 1,203 bytes |
| コンパイル時間 | 1,610 ms |
| コンパイル使用メモリ | 182,572 KB |
| 実行使用メモリ | 96,256 KB |
| 最終ジャッジ日時 | 2024-09-17 12:29:09 |
| 合計ジャッジ時間 | 6,846 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 4 |
| other | AC * 43 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i))
#define all(x) (x).begin(),(x).end()
#define pb push_back
#define fi first
#define se second
using pi = pair<int,int>;
using P = pair<pi,int>;
int n;
ll x;
int D;
vector<int> d;
void factor()
{
int X=x+1;
d.pb(X);
for(int i=2; i*i<=X; ++i)
{
if(X%i==0)
{
d.pb(i);
d.pb(X/i);
}
}
sort(all(d));
d.erase(unique(all(d)),d.end());
D = d.size();
}
map<P,int> dp;
int f(int rem, int didx, int val)
{
if(didx==D) return (rem==0 && val==1);
if(d[didx]>val)
{
if(rem==0) return (val==1);
return 0;
}
ll min_val=1;
rep(i,rem)
{
min_val*=d[didx];
if(min_val>val) return 0;
}
P p = P(pi(rem,didx),val);
if(dp.count(p)) return dp[p];
int ret=0;
ll div=1;
for(int i=0; i<=rem; ++i)
{
if(val%div==0) ret += f(rem-i, didx+1, val/div);
else break;
div *= d[didx];
}
dp[p] = ret;
return ret;
}
int main()
{
cin >>n >>x;
factor();
cout << f(n,0,x+1) << endl;
return 0;
}