BZOJ 3236 AHOI 2013 作业 莫队算法
题目大意:给出一些数,问在一个区间中不同的数值有多少种,和在一个区间中不同的数值有多少个。
思路:因为没有改动,所以就想到了莫队算法。然后我写了5K+的曼哈顿距离最小生成树,然后果断T了。(100s的时限啊,刷status都要刷疯了..,结果最后加了手写读入也没能A)。
后来果断放弃,写了分块版的莫队算法。
84sAC。。。这题卡的。。貌似莫队并非正解。
其有用分块来写莫队就非常easy了。仅仅须要将全部询问的区间排序。左端点所在块作为第一键值,右端点作为第二季键值排序,之后就能够转移了。理论上来时还是曼哈顿距离最小生成树比較快。可是常数实在是太大。随处可见的pair。sort。。
还不太好写。
。
唉。
。
时代的眼泪。。。
CODE(84816AC):
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define MAX 1100010
using namespace std; struct Ask{
int x,y,_id;
int block;
int a,b; bool operator <(const Ask &a)const {
if(block == a.block) return y < a.y;
return block < a.block;
}
void Read(int p) {
scanf("%d%d%d%d",&x,&y,&a,&b);
_id = p;
}
}ask[MAX]; int cnt,asks;
pair<int,int *> xx[MAX << 1];
int src[MAX]; int t;
int fenwick[MAX << 1],_fenwick[MAX];
int num[MAX];
pair<int,int> ans[MAX]; inline int GetSum(int fenwick[],int x)
{
int re = 0;
for(; x; x -= x&-x)
re += fenwick[x];
return re;
} inline void Fix(int fenwick[],int x,int c)
{
for(; x <= t; x += x&-x)
fenwick[x] += c;
} int main()
{
cin >> cnt >> asks;
int temp = 0;
for(int i = 1; i <= cnt; ++i)
scanf("%d",&xx[++temp].first),xx[temp].second = &src[i];
for(int i = 1; i <= asks; ++i) {
ask[i].Read(i);
xx[++temp] = make_pair(ask[i].a,&ask[i].a);
xx[++temp] = make_pair(ask[i].b,&ask[i].b);
}
sort(xx + 1,xx + temp + 1);
for(int i = 1; i <= temp; ++i) {
if(!t || xx[i].first != xx[i - 1].first)
++t;
*xx[i].second = t;
}
int block = static_cast<int>(sqrt(cnt));
for(int i = 1; i <= asks; ++i)
ask[i].block = ask[i].x / block;
sort(ask + 1,ask + asks + 1);
int l = 1,r = 0;
for(int x = 1; x <= asks; ++x) {
while(r < ask[x].y) {
++num[src[++r]];
Fix(fenwick,src[r],1);
if(num[src[r]] == 1) Fix(_fenwick,src[r],1);
}
while(l > ask[x].x) {
++num[src[--l]];
Fix(fenwick,src[l],1);
if(num[src[l]] == 1) Fix(_fenwick,src[l],1);
}
while(r > ask[x].y) {
--num[src[r]];
Fix(fenwick,src[r],-1);
if(!num[src[r]]) Fix(_fenwick,src[r],-1);
--r;
}
while(l < ask[x].x) {
--num[src[l]];
Fix(fenwick,src[l],-1);
if(!num[src[l]]) Fix(_fenwick,src[l],-1);
++l;
}
int p = ask[x]._id;
ans[p].first = GetSum(fenwick,ask[x].b) - GetSum(fenwick,ask[x].a - 1);
ans[p].second = GetSum(_fenwick,ask[x].b) - GetSum(_fenwick,ask[x].a - 1);
}
for(int i = 1; i <= asks; ++i)
printf("%d %d\n",ans[i].first,ans[i].second);
return 0;
}
CODE(TLE):
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define MAX 1100010
#define INF 0x3f3f3f3f
using namespace std; inline int GetInt()
{
int res=0;
char t=getchar();
for (;t>'9'||t<'0';t=getchar());
for (;t>='0'&&t<='9';t=getchar()) res=res*10+t-'0';
return res;
} struct Ask{
int x,y,_id;
int a,b; bool operator <(const Ask &a)const {
if(x == a.x) return y < a.y;
return x < a.x;
}
void Read() {
x = GetInt();
y = GetInt();
a = GetInt();
b = GetInt();
}
}ask[MAX]; struct Edge{
int x,y,length; Edge(int _,int __,int ___):x(_),y(__),length(___) {}
Edge() {}
bool operator <(const Edge &a)const {
return length < a.length;
}
}edge[MAX << 3]; int cnt,asks;
int src[MAX]; int y_x[MAX],t,gt,edges;
pair<int,int *> xx[MAX << 1];
int fenwick[MAX << 1],_fenwick[MAX]; int father[MAX]; int head[MAX],total;
int next[MAX << 1],aim[MAX << 1]; int num[MAX];
pair<int,int> ans[MAX]; inline int CalcM(const Ask &a,const Ask &b)
{
return abs(a.x - b.x) + abs(a.y - b.y);
} inline int GetPos(int x)
{
int re = 0;
for(; x <= t; x += x&-x)
if(ask[fenwick[x]].x + ask[fenwick[x]].y < ask[re].x + ask[re].y)
re = fenwick[x];
return re;
} inline void Fix(int x,int pos)
{
for(; x; x -= x&-x)
if(ask[fenwick[x]].x + ask[fenwick[x]].y > ask[pos].x + ask[pos].y)
fenwick[x] = pos;
} void MakeGraph()
{
for(int dir = 1; dir <= 4; ++dir) {
if(dir == 2 || dir == 4)
for(int i = 1; i <= asks; ++i)
swap(ask[i].x,ask[i].y);
if(dir == 3)
for(int i = 1; i <= asks; ++i)
ask[i].y *= -1;
sort(ask + 1,ask + asks + 1);
for(int i = 1; i <= asks; ++i)
xx[i] = make_pair(ask[i].y - ask[i].x,&y_x[i]);
sort(xx + 1,xx + asks + 1);
t = 0;
for(int i = 1; i <= asks; ++i) {
if(!t || xx[i].first != xx[i - 1].first)
++t;
*xx[i].second = t;
}
memset(fenwick,0,sizeof(fenwick));
for(int i = asks; i; --i) {
int temp = GetPos(y_x[i]);
if(temp)
edge[++edges] = Edge(ask[temp]._id,ask[i]._id,CalcM(ask[temp],ask[i]));
Fix(y_x[i],i);
}
}
for(int i = 1; i <= asks; ++i)
ask[i].x *= -1;
} int Find(int x)
{
if(father[x] == x) return x;
return father[x] = Find(father[x]);
} inline void Add(int x,int y)
{
next[++total] = head[x];
aim[total] = y;
head[x] = total;
} void MST()
{
sort(edge + 1,edge + edges + 1);
for(int i = 1; i <= edges; ++i) {
int fx = Find(edge[i].x);
int fy = Find(edge[i].y);
if(fx != fy) {
father[fx] = fy;
Add(fx,fy),Add(fy,fx);
}
}
} inline int GetSum(int fenwick[],int x)
{
int re = 0;
for(; x; x -= x&-x)
re += fenwick[x];
return re;
} inline void Fix(int fenwick[],int x,int c)
{
for(; x <= gt; x += x&-x)
fenwick[x] += c;
} void DFS(int x,int last)
{
static int l = 1,r = 0; while(r < ask[x].y) {
++num[src[++r]];
Fix(fenwick,src[r],1);
if(num[src[r]] == 1) Fix(_fenwick,src[r],1);
}
while(l > ask[x].x) {
++num[src[--l]];
Fix(fenwick,src[l],1);
if(num[src[l]] == 1) Fix(_fenwick,src[l],1);
}
while(r > ask[x].y) {
--num[src[r]];
Fix(fenwick,src[r],-1);
if(!num[src[r]]) Fix(_fenwick,src[r],-1);
--r;
}
while(l < ask[x].x) {
--num[src[l]];
Fix(fenwick,src[l],-1);
if(!num[src[l]]) Fix(_fenwick,src[l],-1);
++l;
} int p = ask[x]._id;
ans[p].first = GetSum(fenwick,ask[x].b) - GetSum(fenwick,ask[x].a - 1);
ans[p].second = GetSum(_fenwick,ask[x].b) - GetSum(_fenwick,ask[x].a - 1); for(int i = head[x]; i; i = next[i]) {
if(aim[i] == last) continue;
DFS(aim[i],x); while(r < ask[x].y) {
++num[src[++r]];
Fix(fenwick,src[r],1);
if(num[src[r]] == 1) Fix(_fenwick,src[r],1);
}
while(l > ask[x].x) {
++num[src[--l]];
Fix(fenwick,src[l],1);
if(num[src[l]] == 1) Fix(_fenwick,src[l],1);
}
while(r > ask[x].y) {
--num[src[r]];
Fix(fenwick,src[r],-1);
if(!num[src[r]]) Fix(_fenwick,src[r],-1);
--r;
}
while(l < ask[x].x) {
--num[src[l]];
Fix(fenwick,src[l],-1);
if(!num[src[l]]) Fix(_fenwick,src[l],-1);
++l;
}
}
} int main()
{
ask[0].x = ask[0].y = INF;
cnt = GetInt();
asks = GetInt();
int temp = 0;
for(int i = 1; i <= cnt; ++i) {
xx[++temp].first = GetInt();
xx[temp].second = &src[i];
}
for(int i = 1; i <= asks; ++i) {
ask[i].Read();
xx[++temp] = make_pair(ask[i].a,&ask[i].a);
xx[++temp] = make_pair(ask[i].b,&ask[i].b);
ask[i]._id = i;
father[i] = i;
}
sort(xx + 1,xx + temp + 1);
for(int i = 1; i <= temp; ++i) {
if(!gt || xx[i].first != xx[i - 1].first)
++gt;
*xx[i].second = gt;
}
MakeGraph();
MST();
memset(fenwick,0,sizeof(fenwick));
DFS(1,0);
for(int i = 1; i <= asks; ++i)
printf("%d %d\n",ans[i].first,ans[i].second);
return 0;
}
BZOJ 3236 AHOI 2013 作业 莫队算法的更多相关文章
- BZOJ 3236 AHOI 2013 作业 莫队+树状数组
BZOJ 3236 AHOI 2013 作业 内存限制:512 MiB 时间限制:10000 ms 标准输入输出 题目类型:传统 评测方式:文本比较 题目大意: 此时己是凌晨两点,刚刚做了Co ...
- 【bzoj3809/bzoj3236】Gty的二逼妹子序列/[Ahoi2013]作业 莫队算法+分块
原文地址:http://www.cnblogs.com/GXZlegend/p/6805252.html bzoj3809 题目描述 Autumn和Bakser又在研究Gty的妹子序列了!但他们遇到了 ...
- 【BZOJ3809/3236】Gty的二逼妹子序列 [Ahoi2013]作业 莫队算法+分块
[BZOJ3809]Gty的二逼妹子序列 Description Autumn和Bakser又在研究Gty的妹子序列了!但他们遇到了一个难题. 对于一段妹子们,他们想让你帮忙求出这之内美丽度∈[a,b ...
- bzoj 4540 [HNOI 2016] 序列 - 莫队算法 - Sparse-Table - 单调栈
题目传送门 传送点I 传送点II 题目大意 给定一个长度为$n$的序列.询问区间$[l, r]$的所有不同的子序列的最小值的和. 这里的子序列是连续的.两个子序列不同当且仅当它们的左端点或右端点不同. ...
- BZOJ 3289: Mato的文件管理[莫队算法 树状数组]
3289: Mato的文件管理 Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 2399 Solved: 988[Submit][Status][Di ...
- bzoj 3757 苹果树(树上莫队算法)
[题意] 有若干个询问,询问路径u,v上的颜色总数,另外有要求a,b,意为将a颜色看作b颜色. [思路] vfk真是神系列233. Quote: 用S(v, u)代表 v到u的路径上的结点的集合. 用 ...
- bzoj 3236: 洛谷 P4396: [AHOI2013]作业 (莫队, 分块)
题目传送门:洛谷P4396. 题意简述: 给定一个长度为\(n\)的数列.有\(m\)次询问,每次询问区间\([l,r]\)中数值在\([a,b]\)之间的数的个数,和数值在\([a,b]\)之间的不 ...
- bzoj 4358 Permu - 莫队算法 - 链表
题目传送门 需要高级权限的传送门 题目大意 给定一个全排列,询问一个区间内的值域连续的一段的长度的最大值. 考虑使用莫队算法. 每次插入一个数$x$,对值域的影响可以分成4种情况: $x - 1$, ...
- bzoj 2038 A-小Z的袜子[hose] - 莫队算法
作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命…… 具体来说,小Z把这N只袜子从1到N编号,然后从 ...
随机推荐
- 在不同版本号hdfs集群之间转移数据
在不同版本号hdfs集群之间转移数据 最简单的办法就是把src集群的数据导到本地,然后起还有一个进程将本地数据传到des集群上去. 只是这有几个问题: 效率减少 占用本地磁盘空间 不能应付实时 ...
- oracle 之 内存—鞭辟近里(三)
oracle 之 内存—鞭辟近里(三) 今天是2013-07-08,今天晚上突然接到一个电话,我的外甥问我的qq是多少,我感觉很吃惊,他长大了.在他现在这个年龄就开始接触网络,我难免有少许担心,希望他 ...
- Filezilla出现中文乱码
使用Filezilla client FTP客户端登陆某些FTP站点会出现中文乱码,原因是FTP服务器端编码与filezilla client端编码不一致造成的,解决方法如下:文件-站点管理-选中要登 ...
- Lichee (五岁以下儿童) sysconfig1.fex 配置系统
sysconfig配置系统,作为一个通用的软件平台,还希望通过它.能够适应用户不同的方案.通过给出一个相应的配置.用户的方案就能够自己主动执行,而不须要改动系统里面的代码,或者又一次给出參数. 一. ...
- SPARK如何使用AKKA实现进程、节点通信
SPARK如何使用AKKA实现进程.节点通信 <深入理解Spark:核心思想与源码分析>一书前言的内容请看链接<深入理解SPARK:核心思想与源码分析>一书正式出版上市 < ...
- AIX 中以并发模式挂载vg
要想以并发模式挂载VG 必须要有/usr/sbin/gsclvmd 这个进程,而些进程是安装HACMP 的必要的文件集bos.clvm.enh,同时gsclvmd 也是由HACMP 启动的,多个节点挂 ...
- Ice-2.1.2在RHEL Server 5.5上的安装
因为项目的需要,服务器上的程序需要使用Ice接口与其它程序通信,对方提供了一个Windows版的工程,我要把它移植到Linux服务器上,既然Ice是跨平台跨语言的中间件,想来移植不是很困难, ...
- sql小技巧 group by datetime类型字段,只取其中的日期部分
工作中经常会遇到,要在sql中查询报表,查询结果要求按照日期来罗列, 或按照天, 或按照月,年. 这个时候我们经常会苦恼,datetime是精确到毫秒的,如果单纯的group by datetime就 ...
- restrictkeyword
今天在移植ffmpeg到opencore时出现一个编译错误: /libavcodec/dsputil.c:545: error: expected ';', ',' or ')' before 'bl ...
- Html5 Device API详解
三.四月曾学习过html5相关知识,并就html5 device api做过一次讲解 课程时长一个小时,预期达到level 200目标,即知道html5 device api是什么,且知道怎么实现 面 ...