BZOJ 2038: [2009国家集训队]小Z的袜子(hose)
生活随笔
收集整理的這篇文章主要介紹了
BZOJ 2038: [2009国家集训队]小Z的袜子(hose)
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
題目描述
與 小B的詢問 類似
對于一段區(qū)間[l,r],任選兩個數(shù)構(gòu)成排列的方案數(shù)A(r-l+1,2)=(r-l+1)*(r-l)
設(shè)num[i]為i在區(qū)間[l,r]內(nèi)重復(fù)的次數(shù),那么選出的兩個數(shù)都是i的方案有num[i]*(num[i]-1)=num[i]2-num[i]種
那么在[l,r]選出相同的數(shù)的方案就有
sigma(num[i]2-num[i])=sigma(num[i]2)-(r-l+1)?
//因為只有在[l,r]出現(xiàn)過的數(shù)字i的num[i]才不為零,所以sigma(num[i])就是區(qū)間長度?
?
#include<complex> #include<cstdio> #include<algorithm> using namespace std; const int N=5e4+7; struct node{int l,r,id; }q[N]; int n,m; int a[N],pos[N],num[N]; long long ansx[N],ansy[N]; int qread() {int x=0;char ch=getchar();while(ch<'0' || ch>'9')ch=getchar();while(ch>='0' && ch<='9'){x=x*10+ch-'0';ch=getchar();}return x; } bool cmp(const node &a,const node &b) {if(pos[a.l]==pos[b.l])return a.r<b.r;return pos[a.l]<pos[b.l]; } long long Gcd(long long a,long long b) {if(!b)return a;return Gcd(b,a%b); } void change(long long &tot,int id,int add) {tot+=2*add*num[a[id]]+1;num[a[id]]+=add; } int main() {scanf("%d%d",&n,&m);for(int i=1;i<=n;i++)a[i]=qread();for(int i=1;i<=m;i++){q[i].l=qread();q[i].r=qread();q[i].id=i;}long long tmp=sqrt(n);for(int i=1;i<=n;i++)pos[i]=(i-1)/tmp+1;sort(q+1,q+m+1,cmp);int a=0,b=0,l=1,r=0;long long tot=0;for(int i=1;i<=m;i++){while(l<q[i].l)change(tot,l++,-1);while(l>q[i].l)change(tot,--l,1);while(r<q[i].r)change(tot,++r,1);while(r>q[i].r)change(tot,r--,-1);if(q[i].l==q[i].r){ansx[q[i].id]=0;ansy[q[i].id]=1;continue;}ansx[q[i].id]=tot-(r-l+1);ansy[q[i].id]=1ll*(q[i].r-q[i].l+1)*(q[i].r-q[i].l);tmp=Gcd(ansx[q[i].id],ansy[q[i].id]);ansx[q[i].id]/=tmp;ansy[q[i].id]/=tmp;}for(int i=1;i<=m;i++)printf("%lld/%lld\n",ansx[i],ansy[i]);return 0; }?
轉(zhuǎn)載于:https://www.cnblogs.com/LeTri/p/8570154.html
總結(jié)
以上是生活随笔為你收集整理的BZOJ 2038: [2009国家集训队]小Z的袜子(hose)的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: “===” 也有不靠谱的时候
- 下一篇: PHP把下划线分隔命名的字符串 转换成驼