题面:
2038: [2009国家集训队]小Z的袜子(hose)
Time Limit: 20 Sec Memory Limit: 259 MB
Submit: 13997 Solved: 6334
[Submit][Status][Discuss]
Description
作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……
具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择。
Input
输入文件第一行包含两个正整数N和M。N为袜子的数量,M为小Z所提的询问的数量。接下来一行包含N个正整数Ci,其中Ci表示第i只袜子的颜色,相同的颜色用相同的数字表示。再接下来M行,每行两个正整数L,R表示一个询问。
Output
包含M行,对于每个询问在一行中输出分数A/B表示从该询问的区间[L,R]中随机抽出两只袜子颜色相同的概率。若该概率为0则输出0/1,否则输出的A/B必须为最简分数。(详见样例)
Sample Input
6 4
1 2 3 3 3 2
2 6
1 3
3 5
1 6
Sample Output
2/5
0/1
1/1
4/15
【样例解释】
询问1:共C(5,2)=10种可能,其中抽出两个2有1种可能,抽出两个3有3种可能,概率为(1+3)/10=4/10=2/5。
询问2:共C(3,2)=3种可能,无法抽到颜色相同的袜子,概率为0/3=0/1。
询问3:共C(3,2)=3种可能,均为抽出两个3,概率为3/3=1/1。
注:上述C(a, b)表示组合数,组合数C(a, b)等价于在a个不同的物品中选取b个的选取方案数。
【数据规模和约定】
30%的数据中 N,M ≤ 5000;
60%的数据中 N,M ≤ 25000;
100%的数据中 N,M ≤ 50000,1 ≤ L < R ≤ N,Ci ≤ N。
嗯刚刚学习了分块
顺便讲了一下莫队
感觉真的很神奇啊%%%
莫队就是处理不带修改的区间查询问题的离线算法
据说是由曼哈顿距离最小生成树而得的
简单来说
就是如果可以将询问按照端点排序
并且可以由此区间端点+-1得到下一个区间
那么结合分块的单点修改和区间查询
可以将左端点排序后分块
再将块内按照右端点进行排序
那么这个复杂度应该是O(n*根号n)
在5*10^4的范围以内都可以
然后小z的袜子这道题
不仅要用莫队
还要用到数学知识···
就是一个公式
取到两只同颜色的袜子的方案数: Σ(i从1到n)cnt[i]∗(cnt[i]−1)(cnt[i]>=2) (cnt[i]表示颜色为i的袜子有几只)
取的总方案数:(r-l+1)*(r-l)
上下一除一约分就好啦
但是代码太复杂不带了自己打了QAQ
所以找了找终于找到一个不鬼畜的比较好懂的代码【虽然有点长】
然后自己加了些注释
如下:
#define ll long long
#include<algorithm>
#include<iostream>
#include<cstdio>
#define N 100010
#include<cmath>
using namespace std;
struct node
{
ll id,l,r,z,m;//z就是分子,m就是分母
}ask[N];
ll num[200010],n,block;
ll a[N],q,tot;
bool comp(node aa,node bb)
{
if(aa.l/block==bb.l/block)
{
return aa.r<bb.r;
}
return aa.l/block<bb.l/block;
}
bool cmp(node aa,node bb)
{
return aa.id<bb.id;
}
void add(ll nu)//加1
{
ll pa=num[nu];
num[nu]++;
ll no=num[nu];
tot-=pa*(pa-1)/2;//修改,减去之前加的那个方案数
tot+=no*(no-1)/2;//加上加一以后的方案数
}
void pop(ll nu)//减1
{
ll pa=num[nu];
num[nu]--;
ll no=num[nu];
tot-=pa*(pa-1)/2;
tot+=no*(no-1)/2;
}
ll gcd(ll a,ll b)
{
return b==0?a:gcd(b,a%b);
}
int main()
{
scanf("%lld",&n);
scanf("%lld",&q);
block=sqrt(n);
for(ll i=1;i<=n;++i) scanf("%lld",&a[i]);
for(ll i=1;i<=q;++i)
{
scanf("%lld%lld",&ask[i].l,&ask[i].r);
ask[i].id=i;
}
sort(ask+1,ask+q+1,comp);
for(ll i=ask[1].l;i<=ask[1].r;++i) add(a[i]);//先朴素处理第一个区间
ask[1].z=tot;
ll s=(ask[1].r-ask[1].l+1);
ask[1].m=s*(s-1)/2;
if(ask[1].z==0) ask[1].m=1;
else
{
ll it=gcd(ask[1].z,ask[1].m);
ask[1].z/=it;
ask[1].m/=it;
}
ll L=ask[1].l,R=ask[1].r;
for(ll i=2;i<=q;++i)
{
while(R<ask[i].r) add(a[++R]);
while(R>ask[i].r) pop(a[R--]);
while(L>ask[i].l) add(a[--L]);
while(L<ask[i].l) pop(a[L++]);
ask[i].z=tot;
ll s=(R-L+1);
ask[i].m=s*(s-1)/2;
if(ask[i].z==0) ask[i].m=1;
else
{
ll it=gcd(ask[i].z,ask[i].m);
ask[i].z/=it;
ask[i].m/=it;//约分
}
}
sort(ask+1,ask+q+1,cmp);
for(ll i=1;i<=q;++i)
printf("%lld/%lld\n",ask[i].z,ask[i].m);
return 0;
}