热门标签 | HotTags
当前位置:  开发笔记 > 编程语言 > 正文

动态规划解不包含相同数字的子串个数问题

比赛描述    仙灵女巫露露,对于魔法的热忱可是超出常人,要是发现了什么上古遗留下的魔法,她总是想方设法地获得,然后研究分析。而最近,他又从邪恶小法师维嘉那里获得了一个

比赛描述

    仙灵女巫露露,对于魔法的热忱可是超出常人,要是发现了什么上古遗留下的魔法,她总是想方设法地获得,然后研究分析。而最近,他又从邪恶小法师维嘉那里获得了一个“奇怪”的魔法卷轴;

    这个魔法卷轴上有一大串数字,而且根据卷轴上的描述,这个魔法的威力指数来自于这一串数字中“魔法区间”的数量;

    所谓“魔法区间”指的是一段连续的闭区间,且这段区间上的所有数字均不相同;

    现在,露露想知道这个魔法的威力指数,你能帮帮她么?



输入

先输入一个正整数T,表示样例个数,1T10。

对于每一个样例,先输入一个正整数n,表示卷轴上的数字个数(1n106);

再输入n个整数,第i个数ai,表示卷轴上第i个数(0ai106)。

输出

对于每个样例,输出一个正整数,即威力指数。

题目保证结果在int范围内。

样例输入

1
3
1 2 3

样例输出

6

提示

  1. 读入数据请使用 scanf();

  2. 对于样例,共有{1},{2},{3},{1,2},{2,3},{1,2,3},6个魔法区间,所以威力为6。


思路:
        这道题数据范围是1000000,理论上允许的最大时间复杂度为nlog(n)以下,首先考虑用动态规划求解。
        动态规划的目标即为找到状态转移方程,假设一共n个数字在数组num[n]中,其中包含前i个数字的解 为answer[i], 那么我们就是要找到状态转移方程T,使
                                                         answer[i+1]=T(answer[i])             (1)
        那么answer[i+1]和answer[i+1]存在怎样的关系呢?

       明显的,前i+1个数字所创造的“魔法空间”肯定是包含了前i个数字创造的“魔法空间”的,并且多出了一部分包含第i个数字的魔法空间。也就是说,(1)式可
进一步写成:
                                                         answer[i+1]=answer[i]+F(i+1)       (2)
        其中F(i+1)为包含第i+1个数字的魔法空间数量。
        
        所以这道题的目标进一步转化成:当我们有i+1个数字时,包含第i+1个数字的魔法空间有多少?

        当然, 我们可以从i+1个数开始回溯,一个一个把数字添加到空集中,直到集合无法保持没有重复数字,那么这个集合的长度就是包含第i+1个数字的序列数(所有序列都依次为后缀子序列)。
        以长度为4的序列 1 2 3 1 为例
        当包含3个数字时, 含有的魔法空间数为6.即 1, 2, 3, 12,23,123 。
        当加入第四个数字后,又增加了231,31,1 (不同于第一个1), 所以1 2 3 1 的魔法空间数为9 。
优化
        求F的过程,可以尽量采取一些聪明的做法。 不断地添加数字并每一次检查序列是否合法是不明智的。
        当我们从第i个数字向前回溯时,显然到一个位置j,且num[j]==num[i]时必须停止,因为集合已经重复了。 那么我们从i+1回溯时,也不可能超过j这个位置。(因为依然包含了num[i]这个数,它与num[j]是重复的)。 我们设第i个数的前一个相同值数字所在位置为mark[i].所以在每一次回溯时,我们可以设定一个值limit(初始为0),表示我们可以达到的最前方的位置。然后在回溯的过程中,对于每一个k(klimit) . 因为我们不可能超过mark[k]这个位置。
         那么回溯将在以下条件结束:
                  达到limit 或者num[k]==num[i+1].
        回溯的距离也就是F(i+1)。 从而我们得到了完整的状态转移方程。

       由以上思路完成的代码如下:
     
#include
#include
using namespace std;
int f(int n)
{
    return (n+1)*n/2;
}
int num[1000000];
int mark[1000000];
int main()
{
    int i,j,n,T;
    scanf("%d",&T);
    //cin>>T;
    while(T--)
    {
        cin>>n;
        memset(mark,0,sizeof(mark));
        for(i=0;i>num[i];
        }
        int sum=0;
        for(i=0;i=limit;j--)
            {
                if(num[j]==num[i])
                    break;
                temp++;
                if(mark[j]>limit)
                    limit=mark[j];
            }
            mark[i]=j+1;
        //    cout<

进一步优化
           不幸的是,以上算法在南邮acm系统上发生了TLE (TIME LIMIT EXCEED)。
          回顾我们的做法和我们的初衷, 会发现如果我们面临的是n个完全不同的数字,那么每次在第i个位置回溯的距离都将是i, 算法的时间复杂度为可怕的On2 .

          看看我们哪里犯傻了呢?

          犯傻的地方就是上一段中的黑体字:
         当我们从第i个数字向前回溯时,显然到一个位置j,且num[j]==num[i]时必须停止,因为集合已经重复了。 那么我们从i+1回溯时,也不可能超过j这个位置。
         

          这个做法并非smart,实际上,我们应该保存的不应该是每一个数字前一个数字的位置,尽管这样可以帮助我们减少回溯的距离,不过仔细想想,在第i+1个位置能回溯的距离,不是仅仅和第i个位置能回溯的距离有关吗?

           我们来整理一下思路,假设我们之前有i个数字,并且已经得到了一个解(动态规划的立足点),那么当加入第i+1个数字时,得到的解可以通过求F(i+1)(包含第i+1个数字的区间数)来得到, 而 这个区间数其实就是包含第i+1个数字的最长的魔法空间的长度。 而这个长度,实际上就是包含第i个数字的最长的魔法空间的长度再加上1!
           没错,对于F来讲, 依然是一个动态规划问题,不过是一个非常简单的动态规划问题:

            F(i+1)=F(i)+1           (3)

           但是我们仍然要拿到mark[num[i+1]], 因为前一个与第i+1个数字相同的数字J可能包含于F(i)所代表的魔法区间中, 而现在因为加入了i+1个数字,所以J会“截断”F(i).

          所以对于这个问题,我们不必再去回溯(也就是说现在的时间复杂度是On ), 只需在考察第i+1个数字时,分析两种情况:

          包含第i个数字的最长空间不含有与i+1 相同的数字: 那么F(i+1)=F(i)+1 

          否则 F(i+1)= i+1 -j .


          注意文章中的变量名称与程序不是完全对应的。程序的实现手段可能更加抽象, 但是思路就是上文所述。

          代码敬上

       

#include
#include
using namespace std;
int num[1000001];
int limit[1000001];
int max_[1000001];
int main()
{
    int i,j,n,T;
    scanf("%d",&T);
//    scanf_(T);
    while(T--)
    {
        scanf("%d",&n);
    //    scanf_(n);
        memset(max_,-1,sizeof(max_));
        for(i=0;i

      






动态规划解不包含相同数字的子串个数问题


推荐阅读
  • HDU 2372 El Dorado(DP)的最长上升子序列长度求解方法
    本文介绍了解决HDU 2372 El Dorado问题的一种动态规划方法,通过循环k的方式求解最长上升子序列的长度。具体实现过程包括初始化dp数组、读取数列、计算最长上升子序列长度等步骤。 ... [详细]
  • 本文讨论了如何优化解决hdu 1003 java题目的动态规划方法,通过分析加法规则和最大和的性质,提出了一种优化的思路。具体方法是,当从1加到n为负时,即sum(1,n)sum(n,s),可以继续加法计算。同时,还考虑了两种特殊情况:都是负数的情况和有0的情况。最后,通过使用Scanner类来获取输入数据。 ... [详细]
  • 动态规划算法的基本步骤及最长递增子序列问题详解
    本文详细介绍了动态规划算法的基本步骤,包括划分阶段、选择状态、决策和状态转移方程,并以最长递增子序列问题为例进行了详细解析。动态规划算法的有效性依赖于问题本身所具有的最优子结构性质和子问题重叠性质。通过将子问题的解保存在一个表中,在以后尽可能多地利用这些子问题的解,从而提高算法的效率。 ... [详细]
  • 本文介绍了使用Java实现大数乘法的分治算法,包括输入数据的处理、普通大数乘法的结果和Karatsuba大数乘法的结果。通过改变long类型可以适应不同范围的大数乘法计算。 ... [详细]
  • 本文介绍了OC学习笔记中的@property和@synthesize,包括属性的定义和合成的使用方法。通过示例代码详细讲解了@property和@synthesize的作用和用法。 ... [详细]
  • 本文介绍了lua语言中闭包的特性及其在模式匹配、日期处理、编译和模块化等方面的应用。lua中的闭包是严格遵循词法定界的第一类值,函数可以作为变量自由传递,也可以作为参数传递给其他函数。这些特性使得lua语言具有极大的灵活性,为程序开发带来了便利。 ... [详细]
  • 本文介绍了C#中数据集DataSet对象的使用及相关方法详解,包括DataSet对象的概述、与数据关系对象的互联、Rows集合和Columns集合的组成,以及DataSet对象常用的方法之一——Merge方法的使用。通过本文的阅读,读者可以了解到DataSet对象在C#中的重要性和使用方法。 ... [详细]
  • 本文详细介绍了Linux中进程控制块PCBtask_struct结构体的结构和作用,包括进程状态、进程号、待处理信号、进程地址空间、调度标志、锁深度、基本时间片、调度策略以及内存管理信息等方面的内容。阅读本文可以更加深入地了解Linux进程管理的原理和机制。 ... [详细]
  • 1,关于死锁的理解死锁,我们可以简单的理解为是两个线程同时使用同一资源,两个线程又得不到相应的资源而造成永无相互等待的情况。 2,模拟死锁背景介绍:我们创建一个朋友 ... [详细]
  • 《数据结构》学习笔记3——串匹配算法性能评估
    本文主要讨论串匹配算法的性能评估,包括模式匹配、字符种类数量、算法复杂度等内容。通过借助C++中的头文件和库,可以实现对串的匹配操作。其中蛮力算法的复杂度为O(m*n),通过随机取出长度为m的子串作为模式P,在文本T中进行匹配,统计平均复杂度。对于成功和失败的匹配分别进行测试,分析其平均复杂度。详情请参考相关学习资源。 ... [详细]
  • 高质量SQL书写的30条建议
    本文提供了30条关于优化SQL的建议,包括避免使用select *,使用具体字段,以及使用limit 1等。这些建议是基于实际开发经验总结出来的,旨在帮助读者优化SQL查询。 ... [详细]
  • 基于layUI的图片上传前预览功能的2种实现方式
    本文介绍了基于layUI的图片上传前预览功能的两种实现方式:一种是使用blob+FileReader,另一种是使用layUI自带的参数。通过选择文件后点击文件名,在页面中间弹窗内预览图片。其中,layUI自带的参数实现了图片预览功能。该功能依赖于layUI的上传模块,并使用了blob和FileReader来读取本地文件并获取图像的base64编码。点击文件名时会执行See()函数。摘要长度为169字。 ... [详细]
  • 本文讨论了Alink回归预测的不完善问题,指出目前主要针对Python做案例,对其他语言支持不足。同时介绍了pom.xml文件的基本结构和使用方法,以及Maven的相关知识。最后,对Alink回归预测的未来发展提出了期待。 ... [详细]
  • 本文主要解析了Open judge C16H问题中涉及到的Magical Balls的快速幂和逆元算法,并给出了问题的解析和解决方法。详细介绍了问题的背景和规则,并给出了相应的算法解析和实现步骤。通过本文的解析,读者可以更好地理解和解决Open judge C16H问题中的Magical Balls部分。 ... [详细]
  • 本文介绍了P1651题目的描述和要求,以及计算能搭建的塔的最大高度的方法。通过动态规划和状压技术,将问题转化为求解差值的问题,并定义了相应的状态。最终得出了计算最大高度的解法。 ... [详细]
author-avatar
若我倆有了愛
这个家伙很懒,什么也没留下!
PHP1.CN | 中国最专业的PHP中文社区 | DevBox开发工具箱 | json解析格式化 |PHP资讯 | PHP教程 | 数据库技术 | 服务器技术 | 前端开发技术 | PHP框架 | 开发工具 | 在线工具
Copyright © 1998 - 2020 PHP1.CN. All Rights Reserved | 京公网安备 11010802041100号 | 京ICP备19059560号-4 | PHP1.CN 第一PHP社区 版权所有