0x07 内容简介与例题习题

《算法竞赛进阶指南》读书笔记汇总
这里面是我在阅读《算法竞赛进阶指南》这本书时的一些思考，有兴趣可以瞧瞧！
如若发现什么问题，可以通过评论或者私信作者提出。希望各位大佬不吝赐教！

贪心的基本思想

贪心是一种在每次决策时采取当前意义下最优策略的算法。一般通过猜测得出策略之后，加以证明即可使用。通常的证明方法有以下几种：
反证法
邻项交换法（通过交换相邻两项，证明原决策比交换之后的决策更优）
决策包容性
范围扩展（证明任何对局部最优策略作用范围的扩展都不会造成整体结果变差）

那我们直接来看几道例题吧。

【例题】防晒（AcWing110）

题目链接
思路：
首先给出贪心策略：按照$minspf$递减的顺序把奶牛排序，依次考虑每只奶牛。对于每头奶牛，找到一瓶还有剩余的且$spf$值最大的防晒霜去用。
接下来我们来证明这个策略的正确性。这里使用范围扩展证明，考虑这一步策略的作用范围扩展到其他奶牛之后的影响。由于奶牛已经按照$minspf$值递减的顺序排序，所以每一瓶大于当前奶牛$minspf$值的防晒霜，都会大于之后奶牛的$minspf$值。那么也就是说， 对于当前奶牛可以使用的任意两瓶防晒霜而言，假设这两瓶防晒霜是$x,y$，如果$sp{f}_x<sp{f}_y$，那么对于之后的奶牛而言，对于这两瓶防晒霜，只有三种可能情况，$x,y$都能用、$x,y$都不能用或者$x$能用$y$不能用。如果$x,y$都不能用或者$x$能用$y$不能用，那么正确性显然。如果是$x,y$都能用的情况，那么满足当前的奶牛与满足之后的奶牛是等价的，因为每只奶牛对答案的贡献是等价的，所以这个策略不会使得结果变差。

AC代码：

#include<bits/stdc++.h>
#define N 2505

using namespace std;

int n,m;
struct node{
    int mi,mx;
    bool operator < (const node& c)const{
        return mi > c.mi;
    }
}cow[N];

struct Node{
    int spf,cover;
    bool operator < (const Node& b) const{
        return spf > b.spf;
    }
}spfs[N];

void solve(){
    cin >> n >> m;
    for(int i = 1;i <= n;i ++)
        cin >> cow[i].mi >> cow[i].mx;
    for(int i = 1;i <= m;i ++)
        cin >> spfs[i].spf >> spfs[i].cover;

    sort(cow + 1,cow + 1 + n);
    sort(spfs + 1,spfs + 1 + m);

    int ans = 0;
    for(int i = 1;i <= n;i ++){
        for(int j = 1;j <= m;j ++){
            if(spfs[j].spf <= cow[i].mx && spfs[j].spf >= cow[i].mi && spfs[j].cover){
                ans ++;
                spfs[j].cover --;
                break;
            }
        }
    }

    cout << ans << endl;
}

int main(){
    solve();
    return 0;
}

【例题】国王游戏（AcWing114）

题目链接
思路：
首先给出贪心策略，按照每个大臣左、右手上数的乘积从小到大排序，就是最优排队方案。
接下来使用邻项交换法证明这个策略的正确性。
对于任意一种排队顺序，设$n$名大臣左右手上的数字分别是$A[1]...A[n]$与$B[1]...B[n]$，国王手里的数字是$A[0]$与$B[0]$
假设我们交换顺序的两位大臣是$i$和$i+1$，那么在交换之前这两位大臣获得的奖励分别为：
$\frac{1}{B[i]}\ast \prod _{j=0}^{i-1}A[j]$与$\frac{A[i]}{B[i+1]}\ast \prod _{j=0}^{i-1}A[j]$
交换之后两位大臣获得的奖励为：
$\frac{1}{B[i+1]}\ast \prod _{j=0}^{i-1}A[j]$与$\frac{A[i+1]}{B[i]}\ast \prod _{j=0}^{i-1}A[j]$
我们要比较交换前后的最大值变换情况，其实就是比较下面两个式子：
$max(\frac{1}{B[i]},\frac{A[i]}{B[i+1]})$与$max(\frac{1}{B[i+1]},\frac{A[i+1]}{B[i]})$
我们把他们同时乘上$B[i]\ast B[i+1]$得到：
$max(B[i+1],B[i]\ast A[i])$与$max(B[i],B[i+1]\ast A[i+1])$
由于所有的数都是正整数，那么$B[i+1]\le B[i+1]\ast A[i+1]$，$B[i]\le B[i]\ast A[i]$
所以只有当$A[i]\ast B[i]\le A[i+1]\ast B[i+1]$，交换前更优。也就是说，只要存在逆序对，那么交换逆序对可以使得情况变优。当无法交换也即按顺序时，为最优局面。得证。
使用这个贪心策略之后，用高精度写一下就好了。（高精好麻烦呜呜呜

AC代码：

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

typedef pair<int, int> PII;

const int N = 1005, base = 10000;

PII nums[N];

void output(vector<int>& a)
{ 
    printf("%d", a.back());
    for (int i = a.size() - 2; i >= 0; i -- )
        printf("%04d", a[i]);
    puts("");
}

vector<int> mul(vector<int> a, int b)
{
    int t = 0;
    for (int i = 0; i < a.size(); i ++ )
    {
        a[i] = a[i] * b + t;
        t = a[i] / base;
        a[i] %= base;
    }
    while(t)
    {
        a.push_back(t % base);
        t /= base;
    }
        
    return a;
}

vector<int> div(vector<int> a, int b)
{
    int t = 0;
    for (int i = a.size() - 1; i >= 0; i -- )
    {
        a[i] += t * base;
        t = a[i] % b;
        a[i] /= b;
    }
    while(a.size() && a.back() == 0) a.pop_back();
    if (a.size() == 0) a.push_back(0);
    return a;
}

vector<int> max_vec(vector<int> a, vector<int> b)
{
    if (a.size() == b.size())
    {
        for (int i = a.size() - 1; i >= 0; i -- )
            if (a[i] > b[i]) return a;
            else if(b[i] > a[i]) return b;
        return a;
    }
    else if(a.size() > b.size()) return a;
    else return b;
}

int main()
{
    int n;
    scanf("%d", &n); 
    for (int i = 0; i <= n; i ++ )
    {
        int a, b;
        scanf("%d%d", &a, &b);
        nums[i] = {a * b, b};
    }
    sort(nums + 1, nums + 1 + n);
    
    vector<int> res(1, 1);
    vector<int> ans(1, 0);
    for (int i = 0; i <= n; i ++ )
    {
        if (i)
    	{
    		ans = max_vec(ans, div(res, nums[i].second));
		}
        res = mul(res, nums[i].first / nums[i].second);
	}
    
    output(ans);
    
    return 0;
}

未完成题目

【例题】给树染色（AcWing115）

【练习】任务（AcWing127）