hdu 4747 线段树/DP
先是线段树
可以知道mex(i,i),mex(i,i+1)到mex(i,n)是递增的。
首先很容易求得mex(1,1),mex(1,2)......mex(1,n)
因为上述n个数是递增的。
然后使用线段树维护,需要不断删除前面的数。
比如删掉第一个数a[1]. 那么在下一个a[1]出现前的 大于a[1]的mex值都要变成a[1]
因为是单调递增的,所以找到第一个 mex > a[1]的位置,到下一个a[1]出现位置,这个区间的值变成a[1].
然后需要线段树实现区间修改和区间求和
#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
#include <time.h>
using namespace std;
const int MAXN = ;
struct Node
{
int l, r;
long long sum;//区间和
int mx;//最大值
int lazy;//懒惰标记,表示赋值为相同的
}segTree[MAXN * ];
void push_up(int i)
{
if (segTree[i].l == segTree[i].r)
{
return;
}
segTree[i].sum = segTree[i << ].sum + segTree[(i << ) | ].sum;
segTree[i].mx = max(segTree[i << ].mx, segTree[(i << ) | ].mx);
}
void Update_Same(int i, int v)
{
segTree[i].sum = (long long)v * (segTree[i].r - segTree[i].l + );
segTree[i].mx = v;
segTree[i].lazy = ;
}
void push_down(int i)
{
if (segTree[i].l == segTree[i].r)
{
return;
}
if (segTree[i].lazy)
{
Update_Same(i << , segTree[i].mx);
Update_Same((i << ) | , segTree[i].mx);
segTree[i].lazy = ;
}
}
int mex[MAXN];
void Build(int i, int l, int r)
{
segTree[i].l = l;
segTree[i].r = r;
segTree[i].lazy = ;
if (l == r)
{
segTree[i].mx = mex[l];
segTree[i].sum = mex[l];
return;
}
int mid = (l + r) >> ;
Build(i << , l, mid);
Build((i << ) | , mid + , r);
push_up(i);
}
//将区间[l,r]的数都修改为v
void Update(int i, int l, int r, int v)
{
if (segTree[i].l == l && segTree[i].r == r)
{
Update_Same(i, v);
return;
}
push_down(i);
int mid = (segTree[i].l + segTree[i].r) >> ;
if (r <= mid)
{
Update(i << , l, r, v);
}
else if (l > mid)
{
Update((i << ) | , l, r, v);
}
else
{
Update(i << , l, mid, v);
Update((i << ) | , mid + , r, v);
}
push_up(i);
}
//得到值>= v的最左边位置!!!!!!!!!!!!!!!!!!!重要
int Get(int i, int v)
{
if (segTree[i].l == segTree[i].r)
{
return segTree[i].l;
}
push_down(i);
if (segTree[i << ].mx > v)
{
return Get(i << , v);
}
else
{
return Get((i << ) | , v);
}
}
int a[MAXN];
map<int, int>mp;
int nextt[MAXN];
int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int n;
while (~scanf("%d", &n) && n)
{
for (int i = ; i <= n; i++)
{
scanf("%d", &a[i]);
}
mp.clear();
int tmp = ;
for (int i = ; i <= n; i++) //先扫一遍得出1-N的MEX 因为是递增的所以tmp初始化一次就行
{
mp[a[i]] = ;
while (mp.find(tmp) != mp.end())
{
tmp++;
}
mex[i] = tmp;
cout << tmp << " ";
}
cout<<endl;
mp.clear();
for (int i = n; i >= ; i--)
{
if (mp.find(a[i]) == mp.end()) //如果找不到后面存在过的
{
nextt[i] = n + ;
}
else
{
nextt[i] = mp[a[i]];
}
mp[a[i]] = i;
}
for(int i=;i<=n;i++)
cout<<nextt[i]<<" ";
cout<<endl;
Build(, , n);
long long sum = ;
for (int i = ; i <= n; i++)
{
sum += segTree[].sum;
if (segTree[].mx > a[i])
{
int l = Get(, a[i]);
int r = nextt[i];
if (l < r)
{
Update(, l, r - , a[i]);//根据分析 l~r-1(下个a[i]出现之前)都要变成a[i];
}
}
Update(, i, i, );
}
printf("%I64d\n", sum);
}
return ;
}
然后是DP!!
首先要明白,以i结束的所有区间的值的和记为f[i]肯定不超过以i+1结束的所有区间的值的和记为f[i+1]。
所以可以根据f[i]间接推出f[i+1],记第i个数为sa[i],显然只用考虑大于等于sa[i]的数j对f[i]=f[i-1]+?的影响,。如果j出现在1~i-1区间中,比较j最晚出现的位置与覆盖完全的1~j-1的最小位置的较小位置k,那么区间j的前一次出现的位置到k位置这个区间内所有点到i位置的值都+1.
这样逐次累加,直到不影响为止。
#include<iostream>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<stack>
#include<list>
#include<queue>
#include<ctime>
#define eps 1e-6
#define INF 0x3fffffff
#define PI acos(-1.0)
#define ll __int64
#define lson l,m,(rt<<1)
#define rson m+1,r,(rt<<1)|1
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std; #define Maxn 210000 int sa[Maxn],pos[Maxn],full[Maxn]; int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
int n; while(scanf("%d",&n)&&n)
{
for(int i=;i<=n;i++)
scanf("%d",&sa[i]);
memset(pos,,sizeof(pos));
memset(full,,sizeof(full));
int last;
ll tt=,ans=; for(int i=;i<=n;i++)
{
if(sa[i]<n)//
{
last=pos[sa[i]];//前一个sa[i]的最晚位置
pos[sa[i]]=i; //最晚位置
for(int j=sa[i];j<n;j++)
{
if(j) //考虑j对前面区间的影响
full[j]=min(full[j-],pos[j]); //
else
full[j]=i;
if(full[j]>last)
tt+=full[j]-last; //last+1到full[j]区间内所有点到i的值+1,逐次累加
else
break;
}
}
printf("i:%d %I64d\n",i,tt);
ans+=tt;
}
printf("%I64d\n",ans);
}
return ;
}
hdu 4747 线段树/DP的更多相关文章
- hdu 4747 线段树
Mex Time Limit: 15000/5000 MS (Java/Others) Memory Limit: 65535/65535 K (Java/Others)Total Submis ...
- HDU 5726 线段树+dp
题意:给出一个序列,后q次询问,求给定区间gcd及整个序列有多少个序列的gcd和这个值相同 首先线段树维护区间gcd建树,之后预处理出每个gcd有多少个子序列,这时需要dp, dp[i][tmp]表示 ...
- HDU 3607 线段树+DP+离散化
题意:从左往右跳箱子,每个箱子有金币数量,只能从矮处向高处跳,求最大可获得金币数,数据规模1<=n<=1e5. 显然是一个dp的问题,不难得出dp[ i ] = max(dp[j] )+v ...
- HDU 3016 Man Down (线段树+dp)
HDU 3016 Man Down (线段树+dp) Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Ja ...
- Tsinsen A1219. 采矿(陈许旻) (树链剖分,线段树 + DP)
[题目链接] http://www.tsinsen.com/A1219 [题意] 给定一棵树,a[u][i]代表u结点分配i人的收益,可以随时改变a[u],查询(u,v)代表在u子树的所有节点,在u- ...
- hdu 5877 线段树(2016 ACM/ICPC Asia Regional Dalian Online)
Weak Pair Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 262144/262144 K (Java/Others)Total ...
- hdu 3974 线段树 将树弄到区间上
Assign the task Time Limit: 15000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) ...
- hdu 3436 线段树 一顿操作
Queue-jumpers Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) To ...
- hdu 3397 线段树双标记
Sequence operation Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Othe ...
随机推荐
- Linux系统格式化命令mke2fs命令简析
1.mke2fs配置文件: # vim /etc/mke2fs.conf [defaults] base_features = sparse_super,filetype,resize_inode,d ...
- python分布式进程
分布式进程指的是将Process进程分布到多台机器上,充分利用多态机器的性能完成复杂的任务 分布式进程在python 中依然要用到multiprocessing 模块.multiprocessing模 ...
- 基于OpenCV的三维数据点的曲面重构_MySurefaceReconstruction
在Opencv中有个Viz模块,可以显示三维物体,还可以实现三维动画,本来是很好的东东,但是里面的函数.类的说明太过简单,始终不得要领.不过其中一个扩展功能非常好,就是你可以在vtk中设计自己的模型类 ...
- jvm的学习笔记:二、类的初始化,代码实战(3)
首次主动此用导致类的初始化 MyParent4 myParent4 = new MyParent4(); MyParent4 myParent5 = new MyParent4(); 输出: MyPa ...
- docker遇到防火墙报错问题解决方法
-- 报错信息[root@localhost docker]# docker run -d -p 5000:5000 training/webapp python app.pycc61442060cb ...
- C++学习笔记-namespace
指标识符的各种可见范围.C++标准程序库中的所有标识符都被定义于一个名为std的namespace中 关于iostream <iostream>和<iostream.h>格式不 ...
- limit和offset、切片操作
#encoding: utf-8 from sqlalchemy import create_engine,Column,Integer,String,Float,func,and_,or_,\ Da ...
- 海量数据和高并发下的 Redis 业务优化实践
本文内容是我在 6 月 23 日参加的深圳 GIAC 技术大会上演讲的文字稿. 观众朋友们,我是来自掌阅的工程师钱文品,掘金小册<Redis 深度历险>的作者.今天我带来的是分享主题是:R ...
- [Python3] 040 文件 一般使用
目录 文件 1. open 函数 2. with 语句 3. 先写再读 3.1 写 3.2 读 4. "位置"的查询与移动 4.1 tell() 4.2 seek(cookie, ...
- multivariate_normal 多元正态分布
多元正态分布 正态分布大家都非常熟悉了,多元正态分布就是多维数据的正态分布,其概率密度函数为 上式为 x 服从 k 元正态分布,x 为 k 维向量:|Σ| 代表协方差矩阵的行列式 二维正态分布概率密度 ...