深度优先搜索+剪枝（2）

LuoGu P2064 奇妙的汽车

题目传送门
虽然标签是DP，但是这也是一道DFS+剪枝的好题哈哈哈哈哈

---

法一：正推

最开始的思路是枚举$[1,s/3]$范围内所有整数$a$，分别搜索。函数中循环2~9，搜索当前总路程，当天路程和天数。但是…时间大到无法想象…

因此我们要剪枝，避免所有不必要的搜索。

首先，对于第一天的路程$a$来说，只有当总路程$s$能够被$a$整除时，才可能满足恰好到达B城市的基本条件。因此我们枚举$[1,\sqrt{s}]$中的所有整数，找到其中$s$的因数，分别搜索，并在这个因数!=1时搜索与之对应的较大因数。

	for(int i = 1; i * i <= s; i++){//循环较小的一半因数
        if(s % i == 0){//是s的因数
            search(i, i, 1);//搜索当前因数
            if(i != 1)search(s/i, s/i, 1);//搜索与当前因数对应的较大的因数
        }
    }

在搜索函数的循环中，每次循环求出一种可能的下一天路程$num$。因为每天路程之间存在倍数关系，因此只有当剩余路程能够被下一天路程整除时，才可能恰好到达B城市。

	for(int i = 2; i <= 9; i++){
        int num = i * x;
        //只有当剩余路程能够被下一天路程整除时，才能保证恰好到达B城市
        if((s - sum) % (num) == 0)search(sum + num, num, days + 1);
    }

返回条件是走过的总路程$sum$大于等于$s$，或者天数大于已求出的最小天数。

	if(sum == s){
        ans = minn(ans, days);
        return ;
    }
    if(days > ans || sum > s)return ;

输出时要注意，只有满足$ans$值发生变动，且$ans\ge 2$(所用天数至少$2$天)时，才输出ans。否则输出$-1$

---

法二：逆推

由题意可知，$s=a+a{k}_1+a{k}_1{k}_2+a{k}_1{k}_2{k}_3+...$，等式两边同时除以$a$，得到$\frac{s}{a}=1+k_1+k_1{k}_2+k_1{k}_2{k}_3+...$ 等式两边同时-1，得到$\frac{s}{a}-1=k_1+k_1{k}_2+k_1{k}_2{k}_3+...$

重复一遍上述操作，得到$\frac{\frac{s}{a}-1}{k_1}-1=k_2+k_2{k}_3+...$

通过上述操作，我们就可以慢慢的将所有的倍数$k$和初始$a$都消掉，同时每消掉一个意味着一天。

模拟一下：

#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;
int s;
int ans = 21;

int minn(int x, int y){
    return x < y ? x : y;
}
void search(int sum, int days){
    if(sum == 1){
        ans = minn(ans, days);
        return ;
    }
    if(days > ans)return ;

    sum -= 1;//既是这一操作的最后一步，也是下一操作的基础
    for(int i = 2; i <= 9; i++){
        if(sum % i == 0){//符合要求的k
            search(sum / i, days + 1);
        }
    }
}

int main(){

    cin >> s;

    //只有当总路程能够被a整除时，才能保证恰好到达B城市
    for(int i = 1; i * i <= s; i++){//循环较小的一半因数
        if(s % i == 0){//符合要求的a
            search(s/i, 1);//搜索s/a
            if(i != 1)search(i, 1);
        }
        
    }

    if(ans != 21 && ans >= 2)cout << ans;
    else cout << "-1";

    return 0;
}