*PAT_甲级_1103 Integer Factorization (30point(s)) (C++)【DFS进阶用法】

最新推荐文章于 2021-02-02 16:34:47 发布
最新推荐文章于 2021-02-02 16:34:47 发布
阅读量228 收藏
点赞数
分类专栏： PAT甲级 文章标签： PAT 甲级 C++ 1103 DFS进阶用法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41528502/article/details/104868651
版权

目录

1，题目描述

题目大意

2，思路

数据结构

算法

3，AC代码

4，解题过程

第一搏

第二搏

第三搏

1，题目描述

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

题目大意

将一个数N分解成K个因子的P次方之和，要求当结果不唯一时，取因子之和最大的，仍不唯一，则将序列按照从大到小排序，选择较大序列（从第一个数往后，直到当前位置的数组A中的数大于数组B中的数；

不存在这样的分解方法时，输出Impossible；

 

2，思路

参考大神的解法：@日沉云起【pat甲级1103. Integer Factorization (30)】

（解题过程的第二搏，第三搏中的思路均可通过测试，这里仅详细介绍最简解法）

数据结构

  • vector<int> tem：暂时存放结果序列；
  • vector<int> ans：结果序列；
  • vector<int> powValue：存储所有不超过N的n^P中的n；
  • DFS函数参数：start指向powValue的特定位置，dep深度（个数），facSum因子之和，powSum指数之和；

算法

  1. 对不超过N的n^P进行打表，并把结果存储到powValue中，DFS遍历时针对此表进行：
  2.  DFS设计：出口（dep >= K || powSum >= N，若满足dep == K && powSum == N && facSum > maxFacSum则更新ans）；遍历（选择当前的start；不选择当前的start，而是更小的数）注意入栈/出栈的时机，以及递归的参数；
  3. 函数调用及输出：

 

3，AC代码

（即第三搏中的代码）

#include<bits/stdc++.h>
using namespace std;
int N, K, P;//N给出的数 K因子个数 P因子的指数
vector<int> tem, ans, powValue;
int maxFacSum = 0;
void initPow(){     //下标对应因子 值对应因子的P次方
    powValue.push_back(0);
    while(powValue.back() < N)
        powValue.push_back(pow(powValue.size(), P));    // !!!妙啊
}
void dfs(int start, int dep, int facSum, int powSum){   //start指向powValue的特定位置 dep深度 facSum因子之和 powSum指数之和
    if(dep >= K || powSum >= N){
        if(dep == K && powSum == N && facSum > maxFacSum){
            maxFacSum = facSum;
            ans = tem;
        }
        return;
    }
    if(start >= 1){
        tem.push_back(start);
        dfs(start, dep+1, facSum+start, powSum+powValue[start]);    //powValue[start]即start的P次方
        tem.pop_back();
        dfs(start-1, dep, facSum, powSum);              //未选择当前的start 而是选择更小的数
    }
}
int main(){
#ifdef ONLINE_JUDGE
#else
    freopen("1.txt", "r", stdin);
#endif // ONLINE_JUDGE

    scanf("%d %d %d\n", &N, &K, &P);
    initPow();
    dfs(powValue.size()-1, 0, 0, 0);                    //从大到小开始选择
    if(ans.size() == K){
        printf("%d = %d^%d", N, ans[0], P);
        for(int i = 1; i < ans.size(); i++)
            printf(" + %d^%d", ans[i], P);
    }else
        printf("Impossible");
    return 0;
}

 

4，解题过程

第一搏

想到了DFS，但是连样例都没运行出来。。。

#include<bits/stdc++.h>
using namespace std;
int N, K, P;//N给出的数 K因子个数 P因子的指数
vector<int> tem, ans;
bool cmp2(int a, int b) {return a > b;}
bool cmp1(vector<int>& A, vector<int>& B){
    int sumA = 0, sumB = 0;
    if(B.size() == 0) return true;
    for(int i = 0; i < K; i++){
        sumA += A[i];
        sumB += B[i];
    }
    if(sumA != sumB) return sumA > sumB;
    sort(A.begin(), A.end(), cmp2);
    sort(B.begin(), B.end(), cmp2);
    for(int i = 0; i < K; i++)
        if(A[i] != B[i])
            return A[i] > B[i];
}

void dfs(int start, int dep, int sum){
    tem.push_back(start);
    sum += (int)pow(start, P);
//    if(sum > N){
//        tem.pop_back();
//        return;
//    }
    if(dep >= K){
        //printf("YES1\n");
        if(sum == N){
            for(int i = 0; i < tem.size(); i++)
                printf(" %d", tem[i]);
            printf("\n");
            if(cmp1(tem, ans))
                ans = tem;
        }

        //printf("YES2\n");
        tem.pop_back();
        return;
    }
    for(int i = start; i <= 20; i++){//P最小为2 N最大为400 20^2=400
        dfs(i, dep+1, sum);
        //tem.pop_back();
    }
    tem.pop_back();
}
int main(){
#ifdef ONLINE_JUDGE
#else
    freopen("1.txt", "r", stdin);
#endif // ONLINE_JUDGE

    scanf("%d %d %d\n", &N, &K, &P);
    for(int i = 1; i <= 20; i++)
        dfs(i, 1, 0);
    //cout<<ans.size()<<endl;
//    for(int i = 0; i < ans.size(); i++)
//        printf(" %d", ans[i]);
    return 0;
}

看着自己写的代码，都感觉时间复杂度要炸开了。。。

 

第二搏

向大神求助，，，

不愧是大神，逻辑简单但十分准确到位：

#include<bits/stdc++.h>
using namespace std;
int N, K, P;//N给出的数 K因子个数 P因子的指数
vector<int> tem, ans;
int maxFacSum = 0;
void dfs(int start, int dep, int facSum, int powSum){//start起点 dep深度 facSum因子之和 powSum指数之和
    if(dep >= K || powSum >= N){
        if(dep == K && powSum == N && facSum > maxFacSum){
            maxFacSum = facSum;
            ans = tem;
        }
        return;
    }
    if(start >= 1){
        tem.push_back(start);
        dfs(start, dep+1, facSum+start, powSum+pow(start*1.0, P*1.0));//选择这个start !!!这里的pow函数不要强制转换为int 会出错
        tem.pop_back();
        dfs(start-1, dep, facSum, powSum);//未选择当前的start 而是选择更小的数
    }
}
int main(){
#ifdef ONLINE_JUDGE
#else
    freopen("1.txt", "r", stdin);
#endif // ONLINE_JUDGE

    scanf("%d %d %d\n", &N, &K, &P);
    dfs((int)pow(N*1.0, 1.0/P)+1, 0, 0, 0);//从P次根号下N从大到小进行选择; +1保证强制转换后不会漏掉数; 
    if(ans.size() == K){
        printf("%d = %d^%d", N, ans[0], P);
        for(int i = 1; i < ans.size(); i++)
            printf(" + %d^%d", ans[i], P);
    }else
        printf("Impossible");
    return 0;
}

 

对一些必要的地方进行了剪枝：

  • 对因子的最大值进行了限制：                                                                                                                    
  • 对因子的最大个数dep，和指数之和powSum进行剪枝：                                                             

精妙之处：

  • 使递归从最大的数开始，这样可以避免当因子之和相等时，一个个比较序列中的数字；
  • 适当剪枝；

有种子弹从耳朵👂边飞过的感觉。。。Σ(っ °Д °;)っ

第三搏

打表法！

对不超过N的n^P进行打表，并把结果存储到powValue中，从而避免重复计算数字i的P次方。

每次从powValue中取值：下标即因子，值即因子对应的P次方；

#include<bits/stdc++.h>
using namespace std;
int N, K, P;//N给出的数 K因子个数 P因子的指数
vector<int> tem, ans, powValue;
int maxFacSum = 0;
void initPow(){     //下标对应因子 值对应因子的P次方
    powValue.push_back(0);
    while(powValue.back() < N)
        powValue.push_back(pow(powValue.size(), P));    // !!!妙啊
}
void dfs(int start, int dep, int facSum, int powSum){   //start指向powValue的特定位置 dep深度 facSum因子之和 powSum指数之和
    if(dep >= K || powSum >= N){
        if(dep == K && powSum == N && facSum > maxFacSum){
            maxFacSum = facSum;
            ans = tem;
        }
        return;
    }
    if(start >= 1){
        tem.push_back(start);
        dfs(start, dep+1, facSum+start, powSum+powValue[start]);    //powValue[start]即start的P次方
        tem.pop_back();
        dfs(start-1, dep, facSum, powSum);              //未选择当前的start 而是选择更小的数
    }
}
int main(){
#ifdef ONLINE_JUDGE
#else
    freopen("1.txt", "r", stdin);
#endif // ONLINE_JUDGE

    scanf("%d %d %d\n", &N, &K, &P);
    initPow();
    dfs(powValue.size()-1, 0, 0, 0);                    //从大到小开始选择
    if(ans.size() == K){
        printf("%d = %d^%d", N, ans[0], P);
        for(int i = 1; i < ans.size(); i++)
            printf(" + %d^%d", ans[i], P);
    }else
        printf("Impossible");
    return 0;
}

注意：在某些编译环境中，pow函数的取值可能会产生误差。比如代码在codeBlock上运行时，针对样例输出（第二搏中的代码，也出现了这样的问题，但是把强制转换去掉后，莫名其妙就OK了）

但是，在PTA平台上提交，没有任何问题：

贼稳！o(*￣▽￣*)ブ

&再见萤火虫&
关注
专栏目录
KLU-1.1.0.tar.gz_FEM_klu_lu factorization_稀疏矩阵 分解
09-22
在科学计算领域，特别是在有限元方法（Finite Element Method，简称FEM）中，高效处理大规模稀疏矩阵是至关重要的。KLU-1.1.0是一个专门针对这类问题的开源库，它提供了C语言实现的LU分解算法，为解决线性系统方程组...
PAT（甲级）1103 Integer Factorization (30point(s))
hjyp的博客
05-10 183
题目 题目链接 思路 这道题是一道用DFS进行选择的题目； 最难的点在于想到把 fac 数组预先处理出来，这样在选择的时候其实是在对 fac 数组下标所代表的数从高往低是否选择进行决策； 假设 n^p <= test， 并且此时的 n 最大，那么我们在进行选择的时候要从 考虑 n, n-1, n-2 ……1 是否应该被选择；若index代表当前在对哪个数进行选择，则 index 从 n 开始 递减到 1； nowk是当前选择了多少个数，sum是已选的和，facsum是已选底数的和； 如果 nowk =
1103 Integer Factorization (30分)
白小飞
11-06 124
The K−P factorization of a positive integer N is to write N as the sum of the P-th power of K positive integers. You are supposed to write a program to find the K−P factorization of N for any positive integers N, K and P. Input Specification: Each input fi
PAT甲级A1103 Integer Factorization （30 分）
vir_lee的博客
02-26 249
The K−PK-PK−P factorization of a positive integer NNN is to write NNN as the sum of the PPP-th power of KKK positive integers. You are supposed to write a program to find the K−PK-PK−P factorization o...
1103 Integer Factorization（30 分）
Chuyuan_li的博客
08-08 409
The K−P factorization of a positive integer N is to write N as the sum of the P-th power of K positive integers. You are supposed to write a program to find the K−P factorization of N for any positive...
1103. Integer Factorization (30)
骑火一川明
02-24 293
The K-P factorization of a positive integer N is to write N as the sum of the P-th power of K positive integers. You are supposed to write a program to find the K-P factorization of N for any positive
Numerical-Analysis.rar_crout_lu factorization
07-15
这里提到的"Numerical-Analysis.rar_crout_lu_factorization"是一个关于矩阵分解的压缩包，包含了一些特定的算法实现，如Crout法、LU分解以及LUL分解等。这些方法都是为了解决线性代数中的基本问题，即如何高效地...
QR.zip_QR factorization_cfr_decomposition_eigenvalue
09-14
a QR decomposition (also called a QR factorization) of a matrix is a decomposition of the matrix into an orthogonal and an upper triangular matrix. QR decomposition is often used to solve the linear ...
rsa.rar_factorization_rsa_rsa 签名_rsa加密_数字签名算法
09-21
数字签名算法是一类包括RSA在内的加密技术，它们提供了一种方法来验证电子文档的完整性和发件人的身份。数字签名通常结合散列函数使用，散列函数可以把任意长度的消息转化为固定长度的散列值，这样即使微小的改动也...
LU_FACTORIZATION_METHODS.rar_数学计算_Visual_C++_
08-12
总之，"LU_FACTORIZATION_METHODS.rar"压缩包中包含的程序应该演示了如何在Visual C++环境中实现LU分解法，从而解决线性代数问题。理解并掌握这些知识点不仅有助于深入理解线性代数，还能为数值计算和编程实践提供...
【PAT甲级 DFS+剪枝 C++】1103 Integer Factorization (30分)
MYMarcoreus的博客
02-02 145
# include <iostream> # include <vector> # include <algorithm> # include <cmath> using namespace std; int N, K, P; // 分别存放指定和、底数个数和次方数 int maxSum = -1; // 存放最大底数和 vector<int> tmp, ans; // 分别存放当前的底数组合and和最大的底数组合（和相等时选择底数大的） vec
PAT 1103 Integer Factorization (30point(s))
momoomoo的博客
05-31 109
Recursion helps a lot. #include <iostream> #include <stdio.h> #include<stdlib.h> #include <math.h> #include <algorithm> #include <vector> #include <map> #include <string> #include <queue> #include <s
1103 Integer Factorization (30 分)
Star's blogs
08-21 141
The K−PK-PK−P factorization of a positive integer NNN is to write NNN as the sum of the PPP-th power of KKK positive integers. You are supposed to write a program to find the K−PK-PK−P factorizati...
1103 Integer Factorization （30 point(s)）
徐伯莱
02-10 184
题解 dfs。 #include&lt;iostream&gt; #include&lt;cmath&gt; #include&lt;vector&gt; using namespace std; int n, k, p, maxFacSum = -1; vector&lt;int&gt; ans, v, tempAns; void init(){ int temp = 0, inde...
1103 Integer Factorization (30)
ysq96的博客
08-01 322
The K-P factorization of a positive integer N is to write N as the sum of the P-th power of K positive integers. You are supposed to write a program to find the K-P factorization of N for any positive...
PAT A1103 Integer Factorization (30point(s))
皮皮瓜的博客
04-21 167
题目链接 思路：使用深度遍历算法进行递归，枚举所有可能的情况，直到找到最优解。 注意点： 1.计算一个数的p次幂时必须自己写一个函数，直接调用math.h的pow（）会有精度损失的问题，因为pow()的返回值是double，向int强制转换时会损失精度，想知道为什么可以了解int 和double 在计算机中的存储方式。 2.由于n,p是固定的，不妨开一个数组vector fac去记录满足i^p&l...
1103 Integer Factorization (30point(s)) - C语言 PAT 甲级
有三只土豆
07-07 145
* 建立不大于 N 的 i ^ K 数组，节约运算时间 * DFS 算法求解，dfs(index 当前 i 值，nowk 当前 k 个数，nowindexsum 当前所有 i 值的和，nowpowsum 当前所有 i ^ K 的和) * nowk 和 nowpowsum 判断终止递归 * nowindexsum 判断字典序列大小 * 因为需要字典序最大，所以 dfs 从最大值开始
PAT A 1103（甲级）
CrazyOnes的博客
09-06 425
1103 Integer Factorization（30 分） 作者: CHEN, Yue 单位: 浙江大学 时间限制: 1200 ms 内存限制: 64 MB 代码长度限制: 16 KB The K−P factorization of a positive integer N is to write N as the sum of the P-th power of K posi...
ioDraw-v3.1.0-win64.exe，流程图、思维导图、甘特图绘制软件，Windows64位免安装版
最新发布
11-10
1、ioDraw 是免费的在线作图软件，用于制作流程图、思维导图、甘特图、白板手绘图、mermaid流程图、海报设计、折线图、柱状图、拓扑图、组织结构图、UML、ER和时序图、平面设计图、电子电路图等，无需注册即可使用 2、替代Visio的免费软件，无使用限制 3、适合程序员、白领、系统架构师、项目经理、平面设计师、产品经理使用 4、v3.1.0变更日志： 新增将文件保存到百度网盘。 修复第三方平台文件列表页显示问题。
整数分解问题研究：IFP_factorization解析
这个问题的困难性是多种加密算法安全性的基础，包括广泛使用的RSA加密算法。 在描述整数分解问题时，通常会提到一个关键的数学概念——质因数分解。质因数分解指的是将一个合数（大于1的自然数，且不是质数的数）...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值