結果

問題 No.473 和と積の和
ユーザー はまやんはまやん
提出日時 2017-03-15 15:19:51
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
RE  
実行時間 -
コード長 1,110 bytes
コンパイル時間 1,497 ms
コンパイル使用メモリ 174,572 KB
実行使用メモリ 7,424 KB
最終ジャッジ日時 2024-07-03 23:54:42
合計ジャッジ時間 9,705 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample RE * 4
other RE * 43
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:43:70: warning: iteration 101 invokes undefined behavior [-Waggressive-loop-optimizations]
   43 |         rep(i, 0, 2020) rep(j, 0, 101) rep(k, 0, 2020) memo[i][j][k] = -1;
      |                                                        ~~~~~~~~~~~~~~^~~~
main.cpp:3:33: note: within this loop
    3 | #define rep(i,a,b) for(int i=a;i<b;i++)
      |                                 ^
main.cpp:43:40: note: in expansion of macro 'rep'
   43 |         rep(i, 0, 2020) rep(j, 0, 101) rep(k, 0, 2020) memo[i][j][k] = -1;
      |                                        ^~~
main.cpp:43:70: warning: iteration 101 invokes undefined behavior [-Waggressive-loop-optimizations]
   43 |         rep(i, 0, 2020) rep(j, 0, 101) rep(k, 0, 2020) memo[i][j][k] = -1;
      |                                                        ~~~~~~~~~~~~~~^~~~
main.cpp:3:33: note: within this loop
    3 | #define rep(i,a,b) for(int i=a;i<b;i++)
      |                                 ^
main.cpp:43:9: note: in expansion of macro 'rep'
   43 |         rep(i, 0, 2020) rep(j, 0, 101) rep(k, 0, 2020) memo[i][j][k] = -1;
      |         ^~~

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<b;i++)




int N, X;
vector<int> d;
vector<int> enumdiv(int n) {
	vector<int> S;
	for (int i = 1; i*i <= n; i++) if (n%i == 0) { S.push_back(i); if (i*i != n) S.push_back(n / i); }
	sort(S.begin(), S.end());
	return S;
}
//-----------------------------------------------------------------
int memo[101][101][101];
int f(int n, int m, int k) { // d[m]以上の数をk個使って、積をnとする組合せ
	if (k == 0 && n == 0) return 1;
	if (k == 0) return 0;
	if (n == 0) return 0;
	if (m == d.size()) return 0;

	if (0 <= memo[m][k][n]) return memo[m][k][n];

	int ret = f(n, m + 1, k);
	int x = d[n], kk = k;
	while (x % d[m] == 0) {
		x /= d[m];
		kk--;
		int nn = lower_bound(d.begin(), d.end(), x) - d.begin();
		ret += f(nn, m + 1, kk);
		if (kk == 0) break;
	}

	return memo[m][k][n] = ret;
}
//-----------------------------------------------------------------
int main() {
	cin >> N >> X;

	d = enumdiv(X + 1);
	rep(i, 0, 2020) rep(j, 0, 101) rep(k, 0, 2020) memo[i][j][k] = -1;
	cout << f(d.size() - 1, 1, N) << endl;
}
0