(PAT 1103) Integer Factorization(가방 문제 해결을 위한 깊이 우선 순위)
3080 단어 ACM 알고리즘 문제
Input Specification:
Each input file contains one test case which gives in a line the three positive integers N (≤400), K (≤N) and P (1
Output Specification:
For each case, if the solution exists, output in the format:
N = n[1]^P + ... n[K]^P
where
n[i]
( i
= 1, ..., K
) is the i
-th factor. All the factors must be printed in non-increasing order. Note: the solution may not be unique. For example, the 5-2 factorization of 169 has 9 solutions, such as 122+42+22+22+12, or 112+62+22+22+22, or more. You must output the one with the maximum sum of the factors. If there is a tie, the largest factor sequence must be chosen -- sequence { a1,a2,⋯,aK } is said to be larger than { b1,b2,⋯,bK } if there exists 1≤L≤K such that ai=bi for ibL.
If there is no solution, simple output
Impossible
. Sample Input 1:
169 5 2
Sample Output 1:
169 = 6^2 + 6^2 + 6^2 + 6^2 + 5^2
Sample Input 2:
169 167 3
Sample Output 2:
Impossible
문제 해결 방법:
1. N을 구성할 수 있는 모든 수의 제곱과 한 그룹에 저장하기
2. DFS를 이용하여 여러 그룹을 훑어보며, 각 수에 대해 선택과 선택의 두 가지 상황이 있다
선택(중복 가능): 선택할 수 없을 때까지 계속 선택
선택 안 함: 다음 수를 계속 반복합니다.
#include
#include
#include
using namespace std;
int N, K, P, maxFacSum = -1;
const int MAXN = 100;
int FAC[MAXN];
int cindex = 0;
vector temp, ans;
int mpow(int number, int p) {
int ans = 1;
for (int i = 1; i <= p; ++i) {
ans *= number;
}
return ans;
}
void initFac() {
int ntemp = 0;
while (ntemp <= N) {
FAC[cindex] = mpow(cindex + 1, P);
ntemp = FAC[cindex];
cindex++;
}
}
void FDFS(int index, int numK, int sum, int facSum) {
if (sum == N && numK == K) {
if (facSum > maxFacSum) {
ans = temp; // ans
maxFacSum = facSum;
}
return;
}
if (numK > K || sum > N) {
return; //
}
if (index >= 0) {
temp.push_back(index + 1);
FDFS(index, numK + 1, sum + FAC[index], facSum + index); //
temp.pop_back();
FDFS(index - 1, numK, sum, facSum); //
}
}
int main() {
scanf("%d %d %d", &N, &K, &P);
initFac();
FDFS(cindex - 1, 0, 0, 0);
if (maxFacSum == -1) {
printf("Impossible
");
}
else {
printf("%d = %d^%d", N, ans[0], P);
for (int i = 1; i < ans.size(); ++i) {
printf(" + %d^%d", ans[i], P);
}
}
return 0;
}