原文地址:http://www.cnblogs.com/GXZlegend/p/6803860.html


题目描述

作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……
具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择。

输入

输入文件第一行包含两个正整数N和M。N为袜子的数量,M为小Z所提的询问的数量。接下来一行包含N个正整数Ci,其中Ci表示第i只袜子的颜色,相同的颜色用相同的数字表示。再接下来M行,每行两个正整数L,R表示一个询问。

输出

包含M行,对于每个询问在一行中输出分数A/B表示从该询问的区间[L,R]中随机抽出两只袜子颜色相同的概率。若该概率为0则输出0/1,否则输出的A/B必须为最简分数。(详见样例)

样例输入

6 4
1 2 3 3 3 2
2 6
1 3
3 5
1 6

样例输出

2/5
0/1
1/1
4/15


题解

莫队算法

若长度为l的区间内含有某数i个,那么随机选择两个数都为i的概率为(i*(i-1))/(l*(l-1))

所以一段区间内两数相等的概率为∑(cnt[i]*(cnt[i]-1))/(l*(l-1))

先不考虑l*(l-1),设原有某数n个,加1对答案的贡献为(n+1)*n-n*(n-1)=2*n,减1对答案的贡献为(n-1)*(n-2)-n*(n-1)=-(2*n-2)。

然后区间平移得到∑(cnt[i]*(cnt[i]-1)),最后和l*(l-1)求一下gcd,变成分数输出即可。

#include <cstdio>
#include <cmath>
#include <algorithm>
#define N 50010
using namespace std;
typedef long long ll;
struct data
{int l , r , b , p;
}a[N];
int c[N];
ll ans1[N] , ans2[N] , cnt[N];
ll gcd(ll a , ll b)
{return b ? gcd(b , a % b) : a;
}
bool cmp(data x , data y)
{return x.b == y.b ? x.r < y.r : x.b < y.b;
}
int main()
{int n , m , i , si , lp = 1 , rp = 0 , now = 0;ll tmp;scanf("%d%d" , &n , &m);si = (int)sqrt(n);for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &c[i]);for(i = 1 ; i <= m ; i ++ ) scanf("%d%d" , &a[i].l , &a[i].r) , a[i].b = (a[i].l - 1) / si , a[i].p = i;sort(a + 1 , a + m + 1 , cmp);for(i = 1 ; i <= m ; i ++ ){while(lp < a[i].l) now -= 2 * cnt[c[lp]] - 2 , cnt[c[lp]] -- , lp ++ ;while(rp > a[i].r) now -= 2 * cnt[c[rp]] - 2 , cnt[c[rp]] -- , rp -- ;while(lp > a[i].l) lp -- , now += 2 * cnt[c[lp]] , cnt[c[lp]] ++ ;while(rp < a[i].r) rp ++ , now += 2 * cnt[c[rp]] , cnt[c[rp]] ++ ;ans1[a[i].p] = now , ans2[a[i].p] = (ll)(a[i].r - a[i].l + 1) * (a[i].r - a[i].l);}for(i = 1 ; i <= m ; i ++ ) tmp = gcd(ans1[i] , ans2[i]) , printf("%lld/%lld\n" , ans1[i] / tmp , ans2[i] / tmp);return 0;
}

转载于:https://www.cnblogs.com/GXZlegend/p/6803860.html

【bzoj2038】[2009国家集训队]小Z的袜子(hose) 莫队算法相关推荐

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

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

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

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

  3. 2038: [2009国家集训队]小Z的袜子(hose)+莫队入门

    题目链接:2038: [2009国家集训队]小Z的袜子(hose) 题目: Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再 ...

  4. BZOJ 2038: [2009国家集训队]小Z的袜子(莫队算法例题)

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

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

    bzoj2038[2009国家集训队]小Z的袜子(hose) 题意: 把N只袜子从1到N编号,每次求从编号为L到R的袜子中抽两只,有多大的概率抽到颜色相同的袜子. 题解: 不知道要用什么数据结构,但是 ...

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

    2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec Memory Limit: 259 MB Submit: 19269 Solved: 8851 [Sub ...

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

    [传送门:BZOJ2038] 简要题意: 给出n只袜子,每只袜子都有颜色 有多个询问,每次询问一个区间L,R,求出在这个区间内选出两只相同颜色袜子的概率,以最简分数形式输出(不用化成整数,如果概率为0 ...

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

    题目传送门 强啊学了一发莫队. 解法: 刚看到这题没啥思路就只有一个n方的垃圾.. %题解说是莫队,学了一发. 莫队大概就是用来优化暴力的一个东西. 这道题的数学解法相必大家都会.. 在n个袜子里面选 ...

  9. BZOJ 2038: [2009国家集训队]小Z的袜子(hose)【莫队算法裸题学习笔记】

    2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec  Memory Limit: 259 MB Submit: 9894  Solved: 4561 [Su ...

最新文章

  1. linux作为办公机器,Ubuntu能否做日常办公使用?
  2. asp.net ajax删除数据,Asp.net MVC 2 使用Ajax删除数据
  3. Vim as Python IDE on windows(转)
  4. 分布式 | Dubbo 架构设计详解
  5. Java的几何布朗运动
  6. 计算机网络【三】物理层数据通信
  7. 基于JAVA+SpringBoot+Vue+Mybatis+MYSQL的图书馆管理系统
  8. 网关是什么?有什么作用?
  9. 使用SubSonic生成数据访问层步骤
  10. netty权威指南学习
  11. BiomaRt 将小鼠的ENTREZID转化为人类的ENTREZID(同源ENTREZID转换)
  12. 力扣1823题:找出游戏获胜者(约瑟夫环)
  13. import()函数和import语句
  14. Javaweb后端阶段
  15. python以读写方式打开文件_Python读写文件
  16. mldonkey 的使用
  17. mysql 的delete from 子查询限制
  18. Jove:@AlanCui 在手机上用QQ或GTalk来贴Fanfou吧, 不用敲空格或省略号. 好像你输入多少个中文字符就要补多少个空格...
  19. 华为快应用IDE:微信小程序转换成快应用功能火热公测啦!
  20. 操作系统不等于 Linux,六问操作系统新时代!

热门文章

  1. AES加密,解决了同步问题,和随机密钥和固定密钥,多端通信加密不一致解决办法...
  2. 【BZOJ3489】A simple rmq problem(树套树)
  3. JS之模板技术(aui / artTemplate)
  4. [Android]Notification汇总
  5. C#中获取程序集版本号的方法
  6. JAVA(时间对比排序程序)
  7. 003 - JQuery APP viewport rem
  8. js中短路运算符 ||
  9. 24期分期免息可以提前还吗?
  10. 你是怎么发现你的同事很有钱的?