RE:从零开始的算法之路第四章

0.关于

如何将实际问题简单化,用循环和判断使机器执行就能得到答案.就是算法的思想
比如暴力,贪心,分治,动态规划等等,这章就是贪心和分治的思想介绍,这2种算法思想容易理解容易编程,
遇到问题可以优先考虑

1.贪心

介绍

把整个问题分为多个步骤在每个步骤都选取当前步骤的最优方案，直到所有步骤结束
这种走一步看一步的方案，会使得"目光短浅"，局部最优的组合不一定是全局最优
代码主要做的就是做一些规划，使得区部最优，达到全局最优
如何判断一个题目能用贪心？
1.当一个问题的最优解包含其子问题的最优解时。那么这个问题具有最优子结构性质,
也就是说，局部和全局最优解都是相同标准的
2.问题的最优解通过一系列局部最优选择得到的

例题

1.1区间覆盖

今年暑假不AC HDU - 2037
“今年暑假不AC？”
“是的。”
“那你干什么呢？”
“看世界杯呀，笨蛋！”
“@#$%^&*%…”

确实如此，世界杯来了，球迷的节日也来了，估计很多ACMer也会抛开电脑，奔向电视了。
作为球迷，一定想看尽量多的完整的比赛，当然，作为新时代的好青年，你一定还会看一些其它的节目，比如新闻联播（永远不要忘记关心国家大事）、非常6+7、超级女生，以及王小丫的《开心辞典》等等，假设你已经知道了所有你喜欢看的电视节目的转播时间表，你会合理安排吗？（目标是能看尽量多的完整节目）
Input
输入数据包含多个测试实例，每个测试实例的第一行只有一个整数n(n<=100)，表示你喜欢看的节目的总数，然后是n行数据，每行包括两个数据Ti_s,Ti_e (1<=i<=n)，分别表示第i个节目的开始和结束时间，为了简化问题，每个时间都用一个正整数表示。n=0表示输入结束，不做处理。
Output
对于每个测试实例，输出能完整看到的电视节目的个数，每个测试实例的输出占一行。
Sample Input
12
1 3
3 4
0 7
3 8
15 19
15 20
10 15
8 18
6 12
5 10
4 14
2 9
0
Sample Output
5
思路

• 贪心策略有3种:最早开始时间,用时最短,最早结束
• 1.若每次选取开始时间最早的节目，无法求出最优解
  2.若每次选取播放时间最短的节目，无法求出最优解
  3.若每次选取结束时间最短的节目，用数学归纳法可以求出最优解
• 用结构体存储每个节目的开始时间和结束时间
• 根据每个节目的结束时间将这个结构体数组从小到大排序
• 如果当前节目的开始时间大于上一个节目的结束时间，则选入改节目
  代码

#include<iostream>
#include <algorithm>
#include<string>
using namespace std;
class TV {
public:
    int start;
    int end;
};
bool px(class TV a, class TV b) {
    if (a.end != b.end) return a.end < b.end;
    return a.start < b.start;

}
int main()
{
    ios::sync_with_stdio(0); cin.tie(0);
    int n,b,k=0,m=0;
    TV s[100];
    while (cin >> n&&n) {
        k = 0;
        m = 0;
        for (int i = 0; i < n; i++)
            cin >> s[i].start >> s[i].end;
        sort(s, s + n, px);     
        for (int i = 0; i < n; i++)
        {
            if (s[i].start >= k) {
                m++;
                k = s[i].end;
            }
        }
        cout << m << endl;
    }
}

1.2最优装载

迷瘴 HDU - 2570
通过悬崖的yifenfei，又面临着幽谷的考验——
幽谷周围瘴气弥漫，静的可怕，隐约可见地上堆满了骷髅。由于此处长年不见天日，导致空气中布满了毒素，一旦吸入体内，便会全身溃烂而死。
幸好yifenfei早有防备，提前备好了解药材料（各种浓度的万能药水）。现在只需按照配置成不同比例的浓度。
现已知yifenfei随身携带有n种浓度的万能药水，体积V都相同，浓度则分别为Pi%。并且知道，针对当时幽谷的瘴气情况，只需选择部分或者全部的万能药水，然后配置出浓度不大于 W%的药水即可解毒。
现在的问题是：如何配置此药，能得到最大体积的当前可用的解药呢？
特别说明：由于幽谷内设备的限制,只允许把一种已有的药全部混入另一种之中（即：不能出现对一种药只取它的一部分这样的操作）。

Input
输入数据的第一行是一个整数C，表示测试数据的组数；
每组测试数据包含2行，首先一行给出三个正整数n,V,W(1<=n,V,W<=100)；
接着一行是n个整数，表示n种药水的浓度Pi%(1<=Pi<=100)。
Output
对于每组测试数据，请输出一个整数和一个浮点数；
其中整数表示解药的最大体积，浮点数表示解药的浓度(四舍五入保留2位小数)；
如果不能配出满足要求的的解药，则请输出0 0.00。
Sample Input
3
1 100 10
100
2 100 24
20 30
3 100 24
20 20 30
Sample Output
0 0.00
100 0.20
300 0.23
思路

• 简单的贪心算法
• 浓度越小优先度越高,不断混合知道浓度超标

代码

#define _CRT_SECURE_NO_WARNINGS
#include<bits/stdc++.h>
const int maxn = 1e5 + 10;
#define re(x) for(int i=0;i<x;++i)
#define rey(x) for(int j=0;j<x;++j)
#define Cheak(x,y) (x<W&&x>=0&&y>=0&&y<H)
typedef long long LL;
using namespace std;
int dir[4][2] = { {-1,0},{0,-1},{1,0},{0,1} };
int T, n, a[maxn], m, * p, flag = 1;
int main()
{

	cin >> T;
	while (T--) {
		int n, V, num = 0;
		double w, p, next = 0;
		cin >> n >> V >> w;
		re(n)cin >> a[i];
		sort(a, a + n);
		for (int i = 0; i < n; i++)
		{
			next += a[i];
			if (next / (i + 1) > w)break;

			num++;
			p = next / (num * 100 );

		}
		if (num == 0)cout << "0 0.00\n";
		else printf("%d %.2lf\n", num * V, p);
	}
}

1.3huffman编码

例题

Entropy POJ - 1521
输入一个字符串，分别用普通ASCII编码(每个字符8bit)和Huffman编码.
输出编码后的长度并输出压缩比
输入样例:
AAAAABCD
输出样例:
64 13 4.9
思路

• 首先要会huffman编码实现
• 然后不要被空格坑到

代码

#include <iostream>
#include <queue>
#include <string>
#include <stdio.h>
using namespace std;
string s;
int k=0;
struct node{
    int w;
    int parent;
    int key;
    friend bool operator < (const node & a,const node & b){
    if(b.w<a.w)return true;
    else return false;
    }
    }alp[60];
int value(char c){
    if(c=='_')return 26;
    else return c-'A';
    }
void huf(){
    int i,j;
    int length=s.length();
    priority_queue<node>q;
    for(i=0;i<length;i++)
    alp[value(s[i])].w++;
    for(i=0;i<=26;i++){
        if(alp[i].w!=0){
        q.push(alp[i]);
        }
    }
    if(q.size()==1)
    printf("%d %d 8.0\n",8*length,length);
    else {
         j=27;
        while(q.size()>1){
            node s1=q.top();
            q.pop();
            node s2=q.top();
            q.pop();
            alp[j].w=s1.w+s2.w;
            alp[s1.key].parent=j;
            alp[s2.key].parent=j;
            q.push(alp[j]);
            j++;
            }
        int res=0;
        for(int ii=0;ii<27;ii++){
             if(alp[ii].w!=0){
                int len=1;
                int par=alp[ii].parent;
                while(par!=j-1){
                    len++;
                    par=alp[par].parent;
 
                    }
                res+=len*alp[ii].w;
                }
            }
        double dd=8*length;
        printf("%d %d %.1lf\n",8*length,res,dd/res);
        }
    }
int main()
{
    while(cin>>s&&s!="END"){
        int i;
        for(i=0;i<60;i++){
            alp[i].w=0;
            alp[i].key=i;
            }
        huf();
        }
 
    return 0;
}
 

2.分治

介绍

人们在遇到一个大的问题时,自然而然的就会想把这个问题分解:各个击破分而治之
使用分治法要有2个前提:
1平衡子问题:子问题的规模大致相同,最好能二分
2独立子问题:子问题之间没有联系,相互独立
局部的优化扩展到全局,分解,解决,合并
分治的思想就是递归的过程
比如快排和汉诺塔

例题

Inversion HDU - 4911
bobo具有序列a1，a2，…，an。 允许他交换两个相邻的数字不超过k次。

求出交换后的最小反转次数。

注意：反转数是对（i，j）的对数，其中1≤i<j≤n且ai> aj。
输入
输入包含多个测试。 对于每个测试：

第一行包含2个整数n，k（1≤n≤105,0≤k≤109）。 第二行包含n个整数a1，a2，…，an（0≤ai≤109）。
输出
对于每个测试：

单个整数表示最小反转数。
样本输入
3 1
2 2 1
3 0
2 2 1
样本输出
1
2
思路

• 在归并排序中,2个部分合并时,左右都是有序的
• 所以一旦左边有数大于右边的,就只要移动右边第一个下标减左边下标的次数
• 在出现这种情况就把移动次数累加起来,记作m
• 当k>=m时,可以完全有序,逆序数为0
• 如果k>m,则还有k-m次就完全有序,所以逆序数为m-k
  代码

#define _CRT_SECURE_NO_WARNINGS
#include<bits/stdc++.h>
const int maxn = 1e5 + 10;
#define re(x) for(int i=0;i<x;++i)
#define rey(x) for(int j=0;j<x;++j)
#define Cheak(x,y) (x<W&&x>=0&&y>=0&&y<H)
typedef long long LL;
using namespace std;
int dir[4][2] = { {-1,0},{0,-1},{1,0},{0,1} };
LL T, n, a[maxn], m, * p, flag = 1, b[maxn];
void Merge(LL l, LL mid, LL r) {
	LL i = l, j = mid + 1, t = 0;
	while (i <= mid && j <= r)
	{
		if (a[i] > a[j])
		{
			b[t++] = a[j++];
			m += mid + 1 - i;
		}
		else b[t++] = a[i++];
	}
	while (j <= r) b[t++] = a[j++];
	while (i <= mid) b[t++] = a[i++];
	re(t)a[l + i] = b[i];

}
void Mergesort(LL l, LL r) {
	LL mid = (l + r) / 2;
	if (l < r)
	{
		Mergesort(l, mid);
		Mergesort(mid + 1, r);
		Merge(l, mid, r);

	}
}
int main()
{
	ios::sync_with_stdio(0); cin.tie(0);
	LL n, k;
	while (cin >> n >> k)
	{
		re(n)cin >> a[i];
		m = 0;
		Mergesort(0, n - 1);
		if (m <= k)cout << "0\n";
		else cout << m - k << endl;
	}
}