bzoj5277: [Usaco2018 Open]Out of Sorts
Posted akcqhzdy
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj5277: [Usaco2018 Open]Out of Sorts相关的知识,希望对你有一定的参考价值。
被tkj大爷艹爆了5555整套模拟赛都是神仙思路题
那么这题题解
还有一个神仙做法,zory巨神在考场上找规律AC,自己都不会证。。我证明了一下(然而这货还是不认可自己的做法)
按照分割点的思路,我们for循环一次,每次找到比当前点小且最远的点,ans+=j-i+1。毫无疑问,当前点的位移到该点之后停止,分割点的产生时间也就是距离。对于一个数,停下当且仅当它撞到了比他大的数,考虑计数和实际走的差别,对于两个数x,y,x在前且<y,设a是x的目标点,b是y的目标点,必有a<b,可以发现对于计数,第一次加上的是x~y y~a,第二次加上的是y~a a~b的贡献,而实际上,第一次加的是x~y y~a a~b,第二次加的是y~a,而这两个值是相等的。
#include<cstdio> #include<iostream> #include<cstring> #include<cstdlib> #include<algorithm> #include<cmath> using namespace std; typedef long long LL; const LL mod=1e9+7; int read() { int x=0,f=1;char ch=getchar(); while(ch<‘0‘||‘9‘<ch){if(ch==‘-‘)f=-1;ch=getchar();} while(‘0‘<=ch&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x*f; } LL a[1100000],up[1100000],lw[1100000]; LL l[1100000],suml[1100000],sumli[1100000]; int main() { int n,x,y; n=read(); for(int i=1;i<=n;i++) x=read(), y=read(), a[x]=y; //----------------------------------------------------- up[n-1]=a[n-1];for(int i=n-2;i>=0;i--)up[i]=max(up[i+1],a[i]); lw[0]=a[0]; for(int i=1;i<n;i++) lw[i]=min(lw[i-1],a[i]); int now=n; for(int i=0;i<n;i++) while(now>lw[i])l[--now]=i; while(now>0)l[--now]=n-1; for(int i=n-1;i>=0;i--) { suml[i]=(suml[i+1]+l[i])%mod; sumli[i]=(sumli[i+1]+l[i]*(n-i)%mod)%mod; } //----------------------------------------------------- LL ans=0,d; for(int i=1;i<=n;i++) { d=(i*(up[i]-lw[i])*(up[i]+lw[i]+1)/2)%mod; ans=(ans+d)%mod; d=i*lw[i]%mod*(up[i]-lw[i])%mod; ans=((ans-d)%mod+mod)%mod; d=(sumli[lw[i]]-sumli[up[i]])%mod-(suml[lw[i]]-suml[up[i]])*(n-up[i])%mod; ans=((ans-d)%mod+mod)%mod; } printf("%lld ",ans); return 0; }
以上是关于bzoj5277: [Usaco2018 Open]Out of Sorts的主要内容,如果未能解决你的问题,请参考以下文章
[BZOJ3535][Usaco2014 Open]Fair Photography
[BZOJ] 3126: [Usaco2013 Open]Photo
bzoj 3540: [Usaco2014 Open]Fair Photography