作为一个生活散漫的人,小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的袜子 莫队+分块的更多相关文章

  1. [日常摸鱼]bzoj2038[2009国家集训队]小Z的袜子-莫队算法

    今天来学了下莫队-这题应该就是这个算法的出处了 一篇别人的blog:https://www.cnblogs.com/Paul-Guderian/p/6933799.html 题意:一个序列,$m$次询 ...

  2. BZOJ 2038 [2009国家集训队]小Z的袜子 莫队

    2038: [2009国家集训队]小Z的袜子(hose) 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=2038 Descriptionw ...

  3. BZOJ2038 2009国家集训队 小Z的袜子(hose) 【莫队】

    BZOJ2038 2009国家集训队 小Z的袜子(hose) Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼 ...

  4. BZOJ2038: [2009国家集训队]小Z的袜子(hose) -- 莫队算法 ,,分块

    2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 3577  Solved: 1652[Subm ...

  5. [BZOJ2038] [2009国家集训队]小Z的袜子(hose) 莫队算法练习

    2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec  Memory Limit: 259 MBSubmit: 10299  Solved: 4685[Sub ...

  6. P1494 [国家集训队]小Z的袜子/莫队学习笔记(误

    P1494 [国家集训队]小Z的袜子 题目描述 作为一个生活散漫的人,小\(Z\)每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小\(Z\)再也无法忍受这恼人的找袜子过程,于是他 ...

  7. [bzoj2038][2009国家集训队]小Z的袜子(hose)_莫队

    小Z的袜子 hose 2009-国家集训队 bzoj-2038 题目大意:给定一个n个袜子的序列,每个袜子有一个颜色.m次询问:每次询问一段区间中每种颜色袜子个数的平方和. 注释:$1\le n,m\ ...

  8. BZOJ2038: [2009国家集训队]小Z的袜子(hose)

    Time Limit: 20 Sec  Memory Limit: 259 MB Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天, ...

  9. bzoj2038: [2009国家集训队]小Z的袜子(hose) [莫队]

    Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……具体来说,小Z把这N只袜 ...

随机推荐

  1. Hive记录-Impala jdbc连接hive和kudu参考

    1.配置环境Eclipse和JDK 2.加载hive jar包或者impala jar包 备注:从CDH集群里面拷贝出来 下载地址:https://www.cloudera.com/downloads ...

  2. golang 常量

    介绍 常量使用const修改 常量在定义时必须初始化 常量定义后不能修改 常量只能修饰bool.数值类型(int.float).string类型 golang中没有硬性规定常量必使用大写字母,但仍然通 ...

  3. Golang入门教程(三)beego 框架安装

    beego 是一个快速开发 Go 应用的 HTTP 框架,他可以用来快速开发 API.Web 及后端服务等各种应用,是一个 RESTful 的框架,主要设计灵感来源于 tornado.sinatra ...

  4. 跨域请求传递Cookie问题

    问题描述 前后端完全分离的项目,前端使用Vue + axios,后端使用SpringMVC,容器为Tomcat. 使用CORS协议解决跨域访问数据限制的问题,但是发现客户端的Ajax请求不会自动带上服 ...

  5. 让overflow:auto页面滚动条出现时不跳动

    今天看到一篇张鑫旭的文章,转载过来.https://www.zhangxinxu.com/wordpress/2015/01/css-page-scrollbar-toggle-center-no-j ...

  6. C#子类重写父类函数的两种方法

    (1)使用Virtual关键字Override从写 父类子类代码如下,不能修改public 为其它权限 public virtual void Clear() { UpdateView(); } pu ...

  7. 10.tesseract

    1.Tesseract-OCR简介  一个Google支持的开源的OCR图文识别开源项目.支持多种语言(我使用的是3.02 版本,支持包括英文,简体中文,繁体中文),支持Windows,Linux,M ...

  8. 信号量Semaphore

    信号量说简单点就是为了线程同步,或者说是为了限制线程能运行的数量. 那它又是怎么限制线程的数量的哩?是因为它内部有个计数器,比如你想限制最多5个线程运行,那么这个计数器的值就会被设置成5,如果一个线程 ...

  9. 虚拟机中Linux安装(转)

    地址:http://blog.csdn.net/u013142781/article/details/50529030 不是每一个程序员都必须玩过linux,只是博主觉得现在的很多服务器都是linux ...

  10. IIS服务器的安全保护措施

    转载自:https://www.williamlong.info/archives/118.html 部分内容做了修改. 通过标注Web服务器安全级别以及可用性的安全策略,网络管理员将能够从容地在不同 ...