bzoj2038
2016-06-10 10:51
155 查看
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=2038
题目大意:作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……
具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择。
题解:莫队咯,实际上以前写过这道题,但是今天zfy又讲了一遍,真可谓温故而知新,对莫队的理解又更深了一些。
莫队网上博客很多,就不说这个了,但是有一个小优化,似乎网上没提到过,就是我们在排序的时候对于偶数块我们,用右边节点位置从大到小排序,基数块,则从小到大,为什么自己想想,据说会比平常块30%哦。
代码:
View Code
题目大意:作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿。终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命……
具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是完整的一双,甚至不在意两只袜子是否一左一右,他却很在意袜子的颜色,毕竟穿两只不同色的袜子会很尴尬。
你的任务便是告诉小Z,他有多大的概率抽到两只颜色相同的袜子。当然,小Z希望这个概率尽量高,所以他可能会询问多个(L,R)以方便自己选择。
题解:莫队咯,实际上以前写过这道题,但是今天zfy又讲了一遍,真可谓温故而知新,对莫队的理解又更深了一些。
莫队网上博客很多,就不说这个了,但是有一个小优化,似乎网上没提到过,就是我们在排序的时候对于偶数块我们,用右边节点位置从大到小排序,基数块,则从小到大,为什么自己想想,据说会比平常块30%哦。
代码:
#include<iostream> #include<cstring> #include<cstdio> #include<algorithm> #include<cmath> #define maxn 50001 #define ll long long using namespace std; int read() { int x=0; char ch; while (ch=(char)getchar(),ch<'0'||ch>'9'); while (x=x*10+ch-'0',ch=(char)getchar(),ch>='0'&&ch<='9'); return x; } struct data {int l,r,id; ll a,b;} a[maxn]; int n,m,ans; int c[maxn],pos[maxn]; ll s[maxn]; bool cmp(data a, data b) { if (pos[a.l]==pos[b.l]) if (pos[a.l]&1) return a.r<b.r; else return a.r>b.r; return pos[a.l]<pos[b.l]; } bool cmp_id(data a, data b) { return a.id<b.id; } ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);} ll sqr(ll x) { return (ll)x*x; } void init() { n=read(); m=read(); for (int i=1; i<=n; i++) c[i]=read(); int kk=int (sqrt(n)); for (int i=1; i<=m; i++) { a[i].l=read(); a[i].r=read(); a[i].id=i; } for (int i=1; i<=n; i++) pos[i]=(i-1) /kk+1; } void updata(int p,int add) { ans-=sqr(s[c[p]]); s[c[p]]+=add; ans+=sqr(s[c[p]]); } void work() { int r=0,l=1; for (int i=1; i<=m; i++) { for (; r<a[i].r; r++) updata(r+1,1); for (; r>a[i].r; r--) updata(r,-1); for (; l<a[i].l; l++) updata(l,-1); for (; l>a[i].l; l--) updata(l-1,1); if (a[i].l==a[i].r) { a[i].a=0; a[i].b=1; continue; } a[i].a=ans-(a[i].r-a[i].l+1); a[i].b=(ll)(a[i].r-a[i].l)*(a[i].r-a[i].l+1); ll k=gcd(a[i].a,a[i].b); a[i].a/=k; a[i].b/=k; } } int main() { init(); sort(a+1,a+m+1,cmp); work(); sort(a+1,a+m+1,cmp_id); for (int i=1; i<=m; i++) { printf("%lld/%lld\n",a[i].a,a[i].b); } }
View Code
相关文章推荐
- 智慧源自发问,问题重于回答。
- <Sicily>Funny Game
- Peaceful Commission_hdu1814_2-sat
- Linux的系统结构
- 2016 UESTC Training for Search Algorithm & String A - Xiper的奇妙历险(1) 八皇后问题、dfs
- 掌握CoordinatorLayout
- svn学习之路-概念篇
- Mozilla 开始为开源项目提供 50 万美元的安全审计基金
- ubuntu超级账户
- 227. Basic Calculator
- UVA-253-Cube painting
- <Hadoop>Pig安装及测试
- Ural 1874 Football Goal
- 顺序栈
- Word Pattern
- 【CSS笔记一】开始学习CSS,为网页添加样式
- C++第二次上机实验-学生成绩
- shell--赋值,判断,循环
- 程序员应该阅读的10本经典书籍
- 用Toolkit来获取Image图像