BZOJ2038 [2009国家集训队]小Z的袜子 莫队+分块
作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……
具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,
甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择
题解:
先推公式:
对于一个区间$[l,r]$,每种袜子的个数$c_i$,其拿到一对相同袜子的种类数等于:
总的种类则为:
所以最终的概率就是:
(博客园的latex老出bug,脑壳痛)
显然,问题是如何求分子,也就是一个区间内,每种数字个数的的平方和
首先想到线段树,但是这个题显然不满足区间加法,普通的线段树无法解决
即,无法通过$[l,m]$与$[m+1,r]$的答案,直接得到$[l,r]$的答案
于是考虑采用莫队算法离线询问,因为我们发现,一旦我们得到了$[l,r]$的答案
我们就可以$O(1)$的得到$[l\pm 1,r]$与$[l,r\pm 1]$的答案,我们按区间排序询问即可
那如何最少次数的转移来得到所有答案呢?
问题就等价于平面上的曼哈顿最小生成树了
曼哈顿最小生成树的算法是$O(nlogn)$的
但在这道题上,根据证明,复杂度会退化
于是考虑使用数列分块,每次移动区间的时候,按照块移动
我们记录当前的询问的$l$和$r$
询问按照$l$所属的块进行排序,同一块内按照$r$排序
#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn=1e5+10;
const int INF=0x3f3f3f3f;
int casn,n,m,k,size;
ll num[maxn],sum[maxn],ans1[maxn],ans2[maxn];
int id[maxn];
struct node{
int l,r,id;
}cmd[maxn];
int cmp(node a,node b){
if(id[a.l]==id[b.l]) return a.r<b.r;
return id[a.l]<id[b.l];
}
void update(ll &now,int pos,int x) {
now-=sum[num[pos]]*sum[num[pos]];
sum[num[pos]]+=x;
now+=sum[num[pos]]*sum[num[pos]];
}
int main() {
scanf("%d%d",&n,&m);
size=sqrt(n);
for(int i=1;i<=n;i++) id[i]=(i-1)/size+1;
for(int i=1;i<=n;i++) scanf("%lld",num+i);
int a,b;
for(int i=1;i<=m;i++){
scanf("%d%d",&a,&b);
cmd[i]=(node){a,b,i};
}
ll now=0,gcd,c,d;
int l=1,r=0;
sort(cmd+1,cmd+1+m,cmp);
for(int i=1;i<=m;i++){
while(l<cmd[i].l)update(now,l++,-1);
while(r>cmd[i].r)update(now,r--,-1);
while(l>cmd[i].l)update(now,--l,1);
while(r<cmd[i].r)update(now,++r,1);
if(cmd[i].l==cmd[i].r){
ans2[cmd[i].id]=1;
continue;
}
c=now-(cmd[i].r-cmd[i].l+1);
d=(ll)(cmd[i].r-cmd[i].l)*(cmd[i].r-cmd[i].l+1);
gcd=__gcd(c,d);
ans1[cmd[i].id]=c/gcd;
ans2[cmd[i].id]=d/gcd;
}
for(int i=1;i<=m;i++) {
printf("%lld/%lld\n",ans1[i],ans2[i]);
}
return 0;
}
没有用读入优化,开了之后应该是不到1000ms
BZOJ2038 [2009国家集训队]小Z的袜子 莫队+分块的更多相关文章
- [日常摸鱼]bzoj2038[2009国家集训队]小Z的袜子-莫队算法
今天来学了下莫队-这题应该就是这个算法的出处了 一篇别人的blog:https://www.cnblogs.com/Paul-Guderian/p/6933799.html 题意:一个序列,$m$次询 ...
- BZOJ 2038 [2009国家集训队]小Z的袜子 莫队
2038: [2009国家集训队]小Z的袜子(hose) 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=2038 Descriptionw ...
- BZOJ2038 2009国家集训队 小Z的袜子(hose) 【莫队】
BZOJ2038 2009国家集训队 小Z的袜子(hose) Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼 ...
- BZOJ2038: [2009国家集训队]小Z的袜子(hose) -- 莫队算法 ,,分块
2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 3577 Solved: 1652[Subm ...
- [BZOJ2038] [2009国家集训队]小Z的袜子(hose) 莫队算法练习
2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec Memory Limit: 259 MBSubmit: 10299 Solved: 4685[Sub ...
- P1494 [国家集训队]小Z的袜子/莫队学习笔记(误
P1494 [国家集训队]小Z的袜子 题目描述 作为一个生活散漫的人,小\(Z\)每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小\(Z\)再也无法忍受这恼人的找袜子过程,于是他 ...
- [bzoj2038][2009国家集训队]小Z的袜子(hose)_莫队
小Z的袜子 hose 2009-国家集训队 bzoj-2038 题目大意:给定一个n个袜子的序列,每个袜子有一个颜色.m次询问:每次询问一段区间中每种颜色袜子个数的平方和. 注释:$1\le n,m\ ...
- BZOJ2038: [2009国家集训队]小Z的袜子(hose)
Time Limit: 20 Sec Memory Limit: 259 MB Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天, ...
- bzoj2038: [2009国家集训队]小Z的袜子(hose) [莫队]
Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……具体来说,小Z把这N只袜 ...
随机推荐
- Hive记录-Impala jdbc连接hive和kudu参考
1.配置环境Eclipse和JDK 2.加载hive jar包或者impala jar包 备注:从CDH集群里面拷贝出来 下载地址:https://www.cloudera.com/downloads ...
- golang 常量
介绍 常量使用const修改 常量在定义时必须初始化 常量定义后不能修改 常量只能修饰bool.数值类型(int.float).string类型 golang中没有硬性规定常量必使用大写字母,但仍然通 ...
- Golang入门教程(三)beego 框架安装
beego 是一个快速开发 Go 应用的 HTTP 框架,他可以用来快速开发 API.Web 及后端服务等各种应用,是一个 RESTful 的框架,主要设计灵感来源于 tornado.sinatra ...
- 跨域请求传递Cookie问题
问题描述 前后端完全分离的项目,前端使用Vue + axios,后端使用SpringMVC,容器为Tomcat. 使用CORS协议解决跨域访问数据限制的问题,但是发现客户端的Ajax请求不会自动带上服 ...
- 让overflow:auto页面滚动条出现时不跳动
今天看到一篇张鑫旭的文章,转载过来.https://www.zhangxinxu.com/wordpress/2015/01/css-page-scrollbar-toggle-center-no-j ...
- C#子类重写父类函数的两种方法
(1)使用Virtual关键字Override从写 父类子类代码如下,不能修改public 为其它权限 public virtual void Clear() { UpdateView(); } pu ...
- 10.tesseract
1.Tesseract-OCR简介 一个Google支持的开源的OCR图文识别开源项目.支持多种语言(我使用的是3.02 版本,支持包括英文,简体中文,繁体中文),支持Windows,Linux,M ...
- 信号量Semaphore
信号量说简单点就是为了线程同步,或者说是为了限制线程能运行的数量. 那它又是怎么限制线程的数量的哩?是因为它内部有个计数器,比如你想限制最多5个线程运行,那么这个计数器的值就会被设置成5,如果一个线程 ...
- 虚拟机中Linux安装(转)
地址:http://blog.csdn.net/u013142781/article/details/50529030 不是每一个程序员都必须玩过linux,只是博主觉得现在的很多服务器都是linux ...
- IIS服务器的安全保护措施
转载自:https://www.williamlong.info/archives/118.html 部分内容做了修改. 通过标注Web服务器安全级别以及可用性的安全策略,网络管理员将能够从容地在不同 ...