1103. Integer Factorization (30)

来源:互联网 发布:cad编程入门 编辑:程序博客网 时间:2024/06/06 13:12

1.题面

https://www.patest.cn/contests/pat-a-practise/1103

2.题意

将一个数字分解成若干正整数幂之和,比如

169 = 6^2 + 6^2 + 6^2 + 6^2 + 5^2
如果存在多种分解,要求输出正整数之和最小的,如果还有多解要求输出字典序最小的

3.思路

貌似暴力搜索加剪枝也可以得到答案,但是搜索保存答案似乎更快,十几毫秒就过了

4.代码

/*****************************************************************    > File Name: cpp_acm.cpp    > Author: Uncle_Sugar    > Mail: uncle_sugar@qq.com    > Created Time: Thu 16 Feb 2017 02:21:19 CST*****************************************************************/# include <cstring># include <iostream># include <iomanip># include <vector>using namespace std;const int size  = 10 + 1000 ; int n, k, p;vector<int>  vct;int vis[size][size];int fa_vis[size][size];int sum[size][size];int iniv[size];int ipow(int n, int k){    int ret = 1;    for (; k; k>>=1, n *= n){        if (k&1) ret *= n;    }    return ret;}int getAns(int _n, int _k){    if (vis[_n][_k] != -1) return vis[_n][_k];    int ret = 0;    for (int k = vct.size()-1; k >= 0; k--){        int i = vct[k];        if (_n - i >= 0 && getAns(_n - i, _k - 1)){            if (ret == 0){                fa_vis[_n][_k] = _n-i;                sum[_n][_k] = sum[_n-i][_k-1] + iniv[i];                ret = 1;            }else if (sum[_n][_k] < sum[_n-i][_k-1] + iniv[i]){                fa_vis[_n][_k] = _n-i;                sum[_n][_k] = sum[_n-i][_k-1] + iniv[i];            }        }    }    return vis[_n][_k] = ret;}int main(){    std::ios::sync_with_stdio(false);cin.tie(0);    memset(vis, -1, sizeof(vis));    fa_vis[0][0] = -1;    cin >> n >> k >> p;    for (int i = 1; i <= n; i++){        int t = ipow(i, p);        iniv[t] = i;        if (t > n) break;        vct.push_back(t);    }    for (int i = 0; i <= n; i++) vis[i][0] = 0;    vis[0][0] = 1;    if (getAns(n, k)){        cout << n;        char flag = 1;        while (fa_vis[n][k] != -1){            int t = fa_vis[n][k];            if (flag) cout << " =", flag = 0;            else cout << " +";            cout << " " << iniv[n-t] << "^" << p;            n = fa_vis[n][k]; --k;        }        cout << endl;    }else {        cout << "Impossible" << endl;    }    return 0;}



0 0