2016级算法第三次上机-C.AlvinZH的奇幻猜想——三次方

905 AlvinZH的奇幻猜想——三次方

思路

中等题。题意简单,题目说得简单,把一个数分成多个立方数的和,问最小立方数个数。

脑子转得快的马上想到贪心,从最近的三次方数往下减,反正有1^3在最后撑着保证减完。不好意思这是错的,因为1,27,64,125...等立方数之间并不是倍数关系,不能构成贪心策略。举个反例:96=64+8+8+8+8=64+27+1+1+1+1+1,答案明显是5,而贪心会算到7。

既然不是贪心,那就是DP了,没毛病。先讲一下常规做法吧,是这样想的:相当于把一个数化成几份,求最小划分的份数,那么把所求数看成背包的总体积,每个立方数的值看作每个物品的体积,价值都看做1,问题就转化为完全背包(因为每个立方数可以取多次)恰好装满求最小的价值,仔细想想。

所以,套一下完全背包的板子吧,但是,这里的一个问题是完全背包装满,这怎么办呢(可能由于本题可以保证装满这个问题不怎么显眼)。举一反三,假设有可能装不满,怎么办?这个问题初始化dp数组时就可以解决。以下方法对于01背包同样适用:

  • 普通01背包or完全背包:初始化为0;
  • 01背包or完全背包装满求价值最小:初始化为一个大数值如 \(INF\) (0x3f3f3f3f);
  • 01背包or完全背包装满求价值最大:初始化为一个小数值如 \(-INF\) (-0x3f3f3f3f);

为什么呢?对于本题,初始化为大数值,如果没装满,最后dp[n]会依然是INF,因为我们每次比较取的都是较小值。第二次练习赛的G题就是01背包装满求最大价值,初始化为最小值即可。(熬夜写了这么多,希望大家看得到QAQ)

这题还有另外的解法,那就是类似打表,先把所有数的最小立方个数通过迭代计算得到,然后 \(O(1)\) 时间取得答案。具体可见队列迭代的参考代码二以及for循环迭代的参考代码三。其实这里面也是有DP的思想,因为每次迭代会用到之前的结果,对于多组数据来讲,同样可以节省时间。简单易懂,可以学习一下。

分析

对于完全背包直接解法,时间复杂度为 \(O(V*∑(V/wi))\) 。

迭代的话复杂度差不了多少,由于多组数据的原因,迭代打表运行总时间显得更短些,不纠结这个。

扩展:完全背包装满问题:HDU 1114。

参考代码一:完全背包装满求最小价值

//
// Created by AlvinZH on 2017/10/24.
// Copyright (c) AlvinZH. All rights reserved.
//

#include <cstdio>
#include <cstring>
#define INF 0x3f3f3f3f
#define MaxSize 1000005

int weight[105];
int ans[MaxSize];

int main()
{
    for (int i = 1; i < 105; ++i)
        weight[i] = i * i * i;

    int n;
    while(~scanf("%d", &n))
    {
        for (int i = 0; i <= n; ++i)//初始化为最大值
            ans[i] = INF;

        ans[0] = 0;

        for (int i = 0; i < 105; ++i) {
            for (int j = weight[i]; j <= n; ++j) {
                if(ans[j] > ans[j-weight[i]] + 1)
                    ans[j] = ans[j-weight[i]] + 1;
            }
        }
        if(ans[n] == n) printf("Oh NO!\n");
        else printf("%d\n", ans[n]);
    }
}

/*
 * 完全背包恰好装满问题,求最小值。
 */

参考代码二:队列迭代

//
// Created by AlvinZH on 2017/10/24.
// Copyright (c) AlvinZH. All rights reserved.
//

#include <cstdio>
#include <cstring>
#include <queue>
#define INF 0x3f3f3f3f
#define MaxSize 1000005
using namespace std;

int weight[105];
int ans[MaxSize];
queue<int> Q;

void init()
{
    memset(ans, INF, sizeof(ans));
    for (int i = 1; i < 105; ++i)
        weight[i] = i * i * i;

    for (int i = 1; i <= 100; ++i) {
        ans[weight[i]] = 1;
        Q.push(weight[i]);
    }
    while(!Q.empty())
    {
        int w = Q.front();
        Q.pop();
        for (int i = 1; i <= 100; ++i) {
            int num = weight[i] + w;
            if(num > 1000000) break;
            if(ans[num] < INF) continue;
            ans[num] = ans[w] + 1;
            Q.push(num);
        }
    }
}

int main()
{
    init();
    int n;
    while(~scanf("%d", &n))
    {
        if(ans[n] == n) printf("Oh NO!\n");
        else printf("%d\n", ans[n]);
    }
}
/*
 * 思路:直接宽搜,把最开始的数扔进队列,反复用队列中的数去更新没更新的数就行了,注意数组别越界。
 */

参考代码三:for循环迭代

/*
 Author: 曾宥崴(13422)
 Result: AC Submission_id: 391756
 Created at: Fri Nov 10 2017 18:26:40 GMT+0800 (CST)
 Problem: 905   Time: 353   Memory: 6612
*/

#include <cstdio>
#include <algorithm>
#include <iostream>
using namespace std;

int f[1000001],n;

int main()
{
    for (int i = 1; i <= 1000000; i++) f[i] = 1000000000;

    f[0] = 0;
    for (int i = 0; i <= 1000000; i++)
        for (int k = 1; i + k * k * k <= 1000000; k++)
            f[i+k*k*k] = min(f[i+k*k*k],f[i] + 1);

    while (scanf("%d",&n) != EOF)
        if (f[n] == n) printf("Oh NO!\n");
        else printf("%d\n",f[n]);

    return 0;
}
时间: 2024-10-03 03:56:38

2016级算法第三次上机-C.AlvinZH的奇幻猜想——三次方的相关文章

2016级算法第六次上机-C.AlvinZH的学霸养成记II

1032 AlvinZH的学霸养成记II 思路 中等题,贪心. 所有课程按照DDL的大小来排序. 维护一个当前时间curTime,初始为0. 遍历课程,curTime加上此课程持续时间d,如果这时curTime大于此课程DDL,表示无法学习此课程,但是我们不减去此课程,而是减去用时最长的那门课程(优先队列队首,课时最长). 贪心: 假设当前课程为B,被替换课程为A,则有A.d≥B.d,A.e≤B.e.既然curTime+A.d≤A.e,那么curTime+B.d≤B.e绝对成立,保证了B的合法性

2016级算法第六次上机-A.Bamboo之寻找小金刚

Bamboo之寻找小金刚 分析 可以抽象为许多连续线段,分别计数左拐和右拐的个数.考察叉积的基础应用. 假设ABC三点构成一个夹角∠ABC,B就是拐点,AC是辅助形成夹角.考虑线段AB和BC形成的向量 sin∠ABC= (AB * BC)/|AB|*|BC| 两个向量的叉乘除以它们的模 所以叉乘可以判断夹角是否大于180°从而确定转向.当然叉积是有方向的,可以自己选择哪条边在前,只要标准统一即可.每三个点组成一组,遍历,分别计数左拐数和右拐数.具体叉积相关操作可以看<算法导论> 注意 常见的一

2016级算法第六次上机-E.Bamboo之吃我一拳

Bamboo之吃我一拳 分析 当两个点的距离<=d时,才可以出拳,想要使得满足出拳条件的点对最少但不为0 寻找最近点对距离,得到的最近距离能够使得可以出拳的组数最少,因为除了最近点对外其他组合均不符合条件. 在一堆点中找到两个点的距离最小,暴力的O(n^2)计算量很恐怖,可以用分治思想把问题变小: 把平面上的点分为两拨,距离最近的两个点只可能出现在:第一堆,第二堆,和两堆2中各自一个点 分解 想象一条垂直线把所给点集分成两拨:所有的点要么在直线左边,要么在其右边.按x坐标升序排列. 解决 划分后

2016级算法第四次上机-B ModricWang的序列问题

1019 ModricWang的序列问题 思路 此题题意非常清晰,给定一个序列,求出最长上升子序列的长度.从数据规模来看,需要\(O(nlogn)\) 的算法. \(O(nlongn)\) 求最长上升子序列的做法如下: 维护一个数组\(f[]\) ,其中\(f[i]\) 表示当前步骤下长度为i的上升子序列的末尾元素的最小值. 需要注意的是,\(f[i]\) 一定是单调递增的,这个结论十分显然,这里就不做证明了. 使用动态规划思想,对于原序列中的每个元素,都拿去更新一次\(f[]\) .假设当前元

2016级算法第四次上机-E.Bamboo and the Ancient Spell

Bamboo and the Ancient Spell 分析 可能英文读题难度比较大,但是只要看到全大写的 "THE LONGEST COMMON SUBSEQUENCE !"应该就清楚这是考什么的了. 最长公共子序列:可以不连续.序列长度很大时,暴力方法非常费时,这也是一道比较经典的<算法导论>上的动态规划题. 设序列X=和Y=的一个最长公共子序列Z=,则: 若xm=yn,则zk=xm=yn且Zk-1是Xm-1和Yn-1的最长公共子序列: 若xm≠yn且zk≠xm ,则

2016级算法第四次上机-G.ModricWang的序列问题 II

1021 ModricWang的序列问题II 思路 此题与上一题区别不是很大,只是增加了一个长度限制,当场通过的人数就少了很多. 大体解题过程与上一题相同.区别在于对\(f[]\) 的操作.没有长度限制的时候,\(f[]\) 的更新策略是立即更新.假设间隔为\(T\),现在由于需要考虑间隔,那么在处理第\(i\) 个元素的时候,就需要看到 前\(i -T\) 个元素生成的\(f[]\) ,而不能受到第\(i-T+1\) 到 \(i-1\) 个元素的干扰.因此,考虑如下操作:每次准备更新\(f[]

2016级算法第三次上机-G.Winter is coming

904 Winter is coming 思路 难题.首先简化问题, \(n\) 个0与 \(m\) 个1排成一列,连续的0不能超过x个,连续的1不能超过y个,求排列方法数. 显然会想到这是动态规划.最快想到的方法是 \(dp[i][j][x][y]\) 表示已经有i个北境兵j个野人参与排列,且末尾有x个连续北境士兵或y个连续野人士兵的方案数.这方法显然是正确的,但是光是 \(dp[200][200][10][10]\) 数组已经十分接近本题内存限制了,保证MLE.状态转移方法是大模拟,四层fo

2016级算法第三次上机-B.Bamboo和巧克力工厂

B Bamboo和巧克力工厂 分析 三条流水线的问题,依然是动态规划,但是涉及的切换种类比较多.比较易于拓展到n条流水线的方式是三层循环,外层是第k个机器手,里面两层代表可切换的流水线 核心dp语句:cost[i][k] = min(cost[i][k], cost[j][k-1]+t[j][i]+p[i][k]) 也可以在A题的基础上详细的列出所有可能的路线切割情况然后找到最小值. 注意本题与A题中t的含义不同. 上机时给出的伪代码 //数组从0开始 const int maxx= 510;

2016级算法第三次上机-F.ModricWang的导弹防御系统

936 ModricWang的导弹防御系统 思路 题意即为:给出一个长度为n的序列,求出其最长不降子序列. 考虑比较平凡的DP做法: 令\(nums[i]\) 表示这个序列,\(f[x]\) 表示以第\(x\)个数为结尾的最长的不降子序列的长度,状态转移方程为: \[ f[i]=(\max{f[j]}+1) \;\;\;\;\;\;\; \mbox{when $nums[i]<=nums[j]$}\\] f中的最大值即为答案. 时间复杂度\(O(n^2)\),空间复杂度\(O(n)\) 当然也可