洛谷 1679.神奇的四次方数

首先就是作者的暴力硬解：

#include<iostream>
#include<stdio.h>
#include<cstring>
#include<cstdlib>
#include<cmath> 
#include<vector>
#include<algorithm>
#include<stack>
#include<queue>
#include<sstream>
#include<numeric>
#include<map>
#include<limits.h>
#include<set>
#define int long long
#define MAX 1005
#define _for(i,a,b) for(int i=a;i<(b);i++)
#define ALL(x) x.begin(),x.end()
using namespace std;
using PII=pair<int, int>;
int n, m;
int counts;
int t[MAX];
int w[MAX];
int f[MAX];
int fastpow(int x, int p) {
    int res = 1;
    for (; p != 0; p /=2) {
        if (p & 1)
            res *= x;
        x = x * x;
    }
    return res;
}
signed main() {
    ios::sync_with_stdio(false);
    cin.tie(NULL); cout.tie(NULL);
    cin >> m;
    n = 1;
    for (int i = 1; i <= m; i++) {
        if (fastpow(i,4) < m && fastpow(i + 1, 4) < m)
            continue;
        else if (fastpow(i, 4) < m && fastpow(i + 1, 4) > m) {
            n++;
            m -= fastpow(i, 4);
            i = 0;
            
        }
        else if (fastpow(i, 4) == m)
            break;
    }
    cout << n;
    

接下来就是对于上面的思路的DP优化：

乍一看题目我们可能看不出是背包问题。

这里讲一下背包问题的本质，也就是对于每个物品选和不选的选择。这里很明显就是这种问题，只不过物品换成了数字。我们可以这样转化：把每一个数的四次方看作是容量，这个时候价值是多少呢？每一个数字的价值全都是1.为什么呢？因为我们要的是最少有几个数字可以拼出这个结果，所以理所当然这个数的个数本身就是价值啊，也就是1。

由于这里求的是最少，所以我们需要对于所有f数组进行扩大初始化处理。

还有一点，也就是对于物品的个数有多少的问题，也就是数字到底要取到那儿。当然是4次根号下m了，这个时候才是极限。

代码中的ceil是向上取整的函数。

#include<iostream>
#include<stdio.h>
#include<cstring>
#include<cstdlib>
#include<cmath> 
#include<vector>
#include<algorithm>
#include<stack>
#include<queue>
#include<sstream>
#include<numeric>
#include<map>
#include<limits.h>
#include<set>
#define int long long
#define MAX 200010
#define _for(i,a,b) for(int i=a;i<(b);i++)
#define ALL(x) x.begin(),x.end()
using namespace std;
using PII=pair<int, int>;
int n, m;
int counts;
int a[MAX];
int f[MAX];

signed main() {
    ios::sync_with_stdio(false);
    cin.tie(NULL); cout.tie(NULL);
    cin >> m;
    n = ceil(sqrt(sqrt(m) + 1));
    for (int i = 1; i <= m; i++)
        f[i] = 1e8;
    for (int i = 1; i <= n; i++)
        a[i] = i * i * i * i;
    for (int i = 1; i <= n; i++) {
        for (int j = a[i]; j <= m; j++) {
            f[j] = min(f[j], f[j - a[i]] + 1);
        }
    }
    cout << f[m];
    return 0;
}