NOIP模拟21:「Median·Game·Park」

T1:Median

  线性筛+桶+随机化(??什么鬼?)。
  首先,题解一句话秀到了我:

考虑输入如此诡异,其实可以看作随机数据

  随机数据??
  这就意味着分布均匀。。
  又考虑到w<=k<=n
  可以用桶了。
  中位数暴力算的话是排序后取中间。
  但是时间明显不允许。只能\(O(n)\)过掉。所以要维护两个中位数指针(k%2==1当然就是一个了)。
  由于数据随机,分布均匀,所以可以直接跳桶。
  笑死,我当时不信还把数据输了出去,发现有的相邻数据差了几百,就这还能直接跳。。。。。好吧,我肤浅了。。。。
  差点没调出来的代码:

200行的煌煌大作QWQ
#include<bits/stdc++.h>
using namespace std;
namespace STD
{
	#define ll long long
	#define rr register
	const int N=1.1e7;
	const int MAXN=1.8e8+3;
	int n,k;
	int w;
	ll cnt,prime[N];
	bool notprime[MAXN];
	double ans;
	int s2[N];
	int ton[N<<1];
	inline void Prime()
	{
		for(rr int i=2;i<MAXN;i++)
		{
			if(!notprime[i]) prime[++cnt]=i;
			for(rr int j=1;j<=cnt&&i*prime[j]<MAXN;j++)
			{
				notprime[i*prime[j]]=1;
				if(!(i%prime[j]))
					break;
			}
		}
	}
	inline int read()
	{
		rr int x_read=0,y_read=1;	
		rr char c_read=getchar();
		while(c_read<'0'||c_read>'9')
		{
			if(c_read=='-') y_read=-1;
			c_read=getchar();
		}	
		while(c_read<='9'&&c_read>='0')
		{
			x_read=(x_read*10)+(c_read^48);
			c_read=getchar();
		}
		return x_read*y_read;
	}
};
using namespace STD;
int main()
{
	Prime();
	n=read(),k=read(),w=read();
	for(rr int i=1;i<=n;i++)
	{
		prime[i]=prime[i]*i%w;
		s2[i]=prime[i]+prime[i/10+1];
	}
	if(k&1)
	{
		int *p=ton+s2[1],l=1,sum=1;
		ton[s2[1]]++;
		for(rr int i=2;i<=k;i++)
		{
			sum++;
			ton[s2[i]]++;
			if(s2[i]<=(p-ton)) l++;
			if(l<((sum>>1)+1))
				while(l<((sum>>1)+1))
				{
					p++;
					while(!(*p)) p++;
					l+=(*p);
				}
			else
				while(l-(*p)>=((sum>>1)+1))
				{
					l-=(*p);
					p--;
					while(!(*p)) p--;
				}	
		}
		ans+=(p-ton);
		for(rr int i=k+1;i<=n;i++)
		{
			ton[s2[i]]++;
			if(s2[i]<=p-ton) l++;
			ton[s2[i-k]]--;
			if(s2[i-k]<=p-ton) l--;
			if(l<((k>>1)+1))
				while(l<((k>>1)+1))
				{
					p++;
					while(!(*p)) p++;
					l+=(*p);
				}
			else
				while(l-(*p)>=((k>>1)+1))
				{
					l-=(*p);
					p--;
					while(!(*p)) p--;
				}
			ans+=(p-ton);
		}
		printf("%.1lf\n",ans);
	}
	else
	{
		int *p1=ton+min(s2[1],s2[2]),*p2=ton+max(s2[1],s2[2]),l2=2,l1=1,sum=2;
		ton[s2[1]]++;
		ton[s2[2]]++;
		for(rr int i=3;i<=k;i++)
		{
			sum++;
			ton[s2[i]]++;
			if(s2[i]<=(p1-ton)) l1++;
			if(s2[i]<=(p2-ton)) l2++;
			if(l2<((sum>>1)+1))
				while(l2<((sum>>1)+1))
				{
					p2++;
					while(!(*p2)) p2++;
					l2+=(*p2);
				}
			else
				while(l2-(*p2)>=((sum>>1)+1))
				{
					l2-=(*p2);
					p2--;
					while(!(*p2)) p2--;
				}
			if(l1<(sum>>1))
				while(l1<(sum>>1))
				{
					p1++;
					while(!(*p1)) p1++;
					l1+=(*p1);
				}
			else
				while(l1-(*p1)>=(sum>>1))
				{
					l1-=(*p1);
					p1--;
					while(!(*p1)) p1--;
				}
		}
		double temp=((p1-ton)+(p2-ton));
		ans+=(temp/2.00);
		for(rr int i=k+1;i<=n;i++)
		{
			ton[s2[i]]++;
			if(s2[i]<=(p1-ton)) l1++;
			if(s2[i]<=(p2-ton)) l2++;
			ton[s2[i-k]]--;
			if(s2[i-k]<=p1-ton) l1--;
			if(s2[i-k]<=p2-ton) l2--;
			if(l2<((k>>1)+1))
				while(l2<((k>>1)+1))
				{
					p2++;
					while(!(*p2)) p2++;
					l2+=(*p2);
				}
			else
				while(l2-(*p2)>=((k>>1)+1))
				{
					l2-=(*p2);
					p2--;
					while(!(*p2)) p2--;
				}
			if(l1<(k>>1))
				while(l1<(k>>1))
				{
					p1++;
					while(!(*p1)) p1++;
					l1+=(*p1);
				}
			else
				while(l1-(*p1)>=(k>>1))
				{
					l1-=(*p1);
					p1--;
					while(!(*p1)) p1--;
				}
			temp=(p1-ton)+(p2-ton);
			ans+=temp/2.00;
		}
		printf("%.1lf\n",ans);
	}
}

  其他的事情
  在讨论时土哥提到了一个叫”对顶堆”的东西来维护中位数,当然,\(O(nlogn)\)会TLE。
  这个东西其实就是维护两个堆,一个大根堆,一个小根堆。
  小根堆里的数全部大于大根堆里的数,这样就有单调性了,相当于排了个序。
  但比直接用排序算法在N上少了个指数。
  当有数进来时,先与两个堆顶比较如果大于大根堆顶就进小根堆,否则进大根堆。
  然后比较两个堆的大小,然后将多的数放进另一个堆即可。
  中位数就是堆顶之和除以2
  至于说k%2==1的情况,你就不要把中位数往堆里放即可。

T2:Game

  考场上一眼看出来就是贪心,直接放了个堆上去,还纳闷为啥这么简单呢。。
  然后T了。
  正解是\(O(nk)\)
  还是桶。一场考试三道题,两道考桶。。。。
  这题有个看起来很明显但是你往往会忽视的性质:

如果你拿进序列的数比当前序列里的最大值还要大,那么它下一轮一会被拿走

  简单到无需证明。。
  但是他会决定你是\(A\)还是\(T\)。
  记得开longlong。

#include<bits/stdc++.h>
using namespace std;
namespace STD
{
	#define ll long long
	#define rr register 	
	#define inf INT_MAX
	const int N=100004;
	const int K=2004;
	int n,k,p;
	int *po;
	ll score[2];
	ll a[N];
	int b[N];
	int read()
	{
		rr int x_read=0,y_read=1;
		rr char c_read=getchar();
		while(c_read<'0'||c_read>'9')
		{
			if(c_read=='-') y_read=-1;
			c_read=getchar();
		}
		while(c_read<='9'&&c_read>='0')
		{
			x_read=(x_read*10)+(c_read^48);
			c_read=getchar();
		}
		return x_read*y_read;
	}
};
using namespace STD;
int main()
{
	n=read(),k=read();
	for(rr int i=1;i<=n;i++)
		a[i]=read();
	while(k--)
	{
		int p=read();
		ll temp=-inf;
		int roun=1;
		int now=p;
		for(rr int i=1;i<=p;i++)
		{
			temp=max(temp,a[i]);
			b[a[i]]++;
		}
		po=b+temp;	
		temp=-inf;
		while(roun<=n)
		{
			if(temp>(po-b))
			{
				score[roun&1]+=temp;
				temp=-inf;
			}
			else
			{
				score[roun&1]+=(po-b);
				int x=*po;
				x--;
				*po=x;
			}
			now++;
			if(now<=n)
			{
				if(a[now]>(po-b))
					temp=a[now];
				else
					b[a[now]]++;
			}
			while(((*po)==0)&&(po>b))
				po--;
			roun++;	
		}
		printf("%lld\n",score[1]-score[0]);
		score[0]=score[1]=0;
	}
}

T3:Park

  还在推方程,先鸽掉好了。

上一篇:快速排序实现


下一篇:boost::gil::median_filter用法的测试程序