Luogu CF1633D，DP

Luogu CF1633D

题目传送门

绿色的01背包模板题~~（可是我花了一个多小时才做对awa）~~

---

最开始我以为对于一个数$i$，如果$1$变化到$i$需要的次数不止一次那么就不能选，但过了好久一直没写对awa…再看看题？发现自己的想法和题意背道而驰
原来…$1$可以经过$x$次操作（$x<=k$）变化到$i$

---

那这不就是一道01背包的板子吗？

对于一个数$i$，$1$变化到$i$需要的次数就是$i$这个“物品”的“体积”，而$c[i]$即为$i$的“价值”，$k$即为背包容量；每个$i$最多只能被选择$1$次。

所以我们定义$1$个$num$数组，$num[i]$即为$1$变化到$i$所需的次数，也就是$i$的“体积”

void judge(){
    for(int i = 1; i <= 1000; i++)
        for(int j = 1; j <= i; j++)
            if(!num[i+i/j])num[i+i/j] = num[i]+1;
}

首先，我们看到题目的操作是$a[i]=a[i]+\lfloor \frac{a[i]}{x}\rfloor$，因此我们在模拟中用$j$模拟$x$，将小于$i$的正整数分别代入进行操作。同时，因为先达成的值会比后达成的值小，因此不能被覆盖。
后来我们发现，如果只有循环操作是不够的，会出现比实际情况更大的错误值
经过模拟，我们发现，在一些情况下，后达成的值在进一步变化成另一值时，所用总次数会小于先达成的值变化为同一值所用的总次数。比如，$7$->$9$需要1->2->3->6->7->9共$5$次，而$8$->$9$则只需1->2->4->8->9共$4$次。因此我们现将$num$数组除$num[1]$和$num[0]$之外的其他位置都初始化为极大值，再在循环操作中不断取最小变化次数

void judge(){
    //将num赋极大值
    memset(num, 0x7f7f7f7f, sizeof(num));
    num[0] = num[1] = 0;//0不参与，1变到1是0次
    for(int i = 1; i <= 1000; i++){
        for(int j = 1; j <= i; j++)
            num[i+i/j] = minn(num[i+i/j], num[i] + 1);//按照题面进行一次操作
    }
}

这样就完美解决了现有问题。

接下来就是$01$背包的板子啦

for(int j = 1; j <= n; j++){
    for(int g = k; g >= num[b[j]]; g--){
        dp[g] = maxx(dp[g], dp[g - num[b[j]]] + c[j]);
    }
}
cout << dp[k] << endl;

---

但是这种方法非常浪费时间，时间复杂度大约为$O(nk)$，也就是$10^9$，远超$2s$的时间限制。因此我们需要对时间进行优化

输出$num$数组后我们发现，$1$变化到$1000$只需要$12$次，因此对于$n$个数，最多只需要$12n$次，也就是说$k$不需要大于$12000$。也就是说$10^6$的范围完全是在浪费时间，白白循环。同时这也意味着，如果$k>=12n$，那么所有的$a[i]$都可以变得与$b[i]$相等，所有的$c[i]$都可以被计入总收益！因此我们在$01$背包前加入了一个判断$k$大小的操作，如果$k>=12n$，那么总收益就是$c$数组中所有元素之和。输出之后就可以直接进行下一组数据的处理了。

if(k >= 12 * n){
    int sum = 0;
    for(int j = 1; j <= n; j++){
        sum += c[j];
    }
    cout << sum << endl;
    continue;
}

---

AC代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <iomanip>

using namespace std;
int t;//t组数据
int n, k;//数组长度和最多操作次数
int b[1001], c[1001], dp[1000001];
int num[2005];//num[i]是指1变到i所需的次数

int maxx(int x, int y){
    return x > y ? x : y;
}
int minn(int x, int y){
    return x < y ? x : y;
}

void judge(){
    //将num赋极大值
    memset(num, 0x7f7f7f7f, sizeof(num));
    num[0] = num[1] = 0;//0不参与，1变到1是0次
    for(int i = 1; i <= 1000; i++){
        for(int j = 1; j <= i; j++)
            num[i+i/j] = minn(num[i+i/j], num[i] + 1);//按照题面进行一次操作
    }
}
void pre(){
    cin >> n >> k;
    memset(b, 0, sizeof(b));
    memset(c, 0, sizeof(c));
    memset(dp, 0, sizeof(dp));
    for(int i = 1; i <= n; i++)
        cin >> b[i];
    for(int i = 1; i <= n; i++)
        cin >> c[i];
}

int main(){

    //freopen("code.in", "r", stdin);
    //freopen("code.out", "w", stdout);

    judge();
    cin >> t;
    for(int i = 1; i <= t; i++){
        pre();
        if(k >= 12 * n){
            int sum = 0;
            for(int j = 1; j <= n; j++){
                sum += c[j];
            }
            cout << sum << endl;
            continue;
        }
        //display();
        //01背包的模板
        for(int j = 1; j <= n; j++){
            for(int g = k; g >= num[b[j]]; g--){
                dp[g] = maxx(dp[g], dp[g - num[b[j]]] + c[j]);
            }
        }
        cout << dp[k] << endl;
    }

    return 0;
}