[CSP-S模拟测试]:Read(数学)

题目描述

热爱看书的你有$N$本书,第$i$本书的种类为$A[i]$。你希望每天能够看一本书,但是不希望连续两天看种类相同的书。为了达成这个条件,你需要选择一些书不看,作为一个好学生,你希望不看的书尽可能少,求最少可以有多少书不看。


输入格式

为了避免输入文件过大,我们采取如下方式生成$A[i]$。
第一行读入两个个整数$M$,$K$。
接下来一行读入$M$个整数$count[i]$,其中$N=sum count[i]$。
接下来一行读入$M$个整数$X[i]$。
接下来一行读入$M$个整数$Y[i]$。
接下来一行读入$M$个整数$Z[i]$。
然后按照下面这段代码生成。

int N = 0, S = (1 << K) - 1;
for (int i = 1; i <= M; ++i) {
	N = N + 1;
	A[N] = X[i];
	long long last = X[i];
	for (int j = 1; j < count[i]; ++j) {
		last = (last * Y[i] + Z[i]) & S;
		N = N + 1;
		A[N] = last;
	}
}

注意:因为生成$A[i]$的方法不好用语言描述,所以用代码来表达。直接照搬代码的话后果自负。


输出格式

输出一个整数表示答案。


样例

样例输入:

4 2
1 1 1 1
1 1 1 2
0 0 0 0
0 0 0 0

样例输出:

1


数据范围与提示

对于所有数据,$Nleqslant 50,000,000,Mleqslant 1,000,Kleqslant 30,X[i],Y[i]<2^k$。
$subtask#1 10pts$:$Nleqslant 20$。
$subtask#2 20pts$:$Nleqslant 1,000,000$。
$subtask#3 20pts$:$Kleqslant 20$。
$subtask#4 50pts$:无特殊限制。


题解

为方便,我们设每种书有$sum[i]$本。

首先,读不到的书的个数为$2 imes max(sum[i])-n-1$。

先来解释一下这个式子,最多的那种书显然有$max(sum[i])$,那么剩下的书有$n-max(sum[i])$种,我们进行插空,最多能插$n-max(sum[i])+1$本,化简一下式子就好了。

那么问题就在于如何求出最多的那本书了。

注意空间只有$16MB$,我们显然不能都求出来,所以我们用两个变量$id,cnt$,$cnt$初始为$0$,然后生成每一个$A[i]$,如果$cnt=0$,那么就令$id=A[i],cnt=1$,否则如果$id=A[i]$则$cnt++$,如果不等于,$cnt--$。最后只要再扫一遍求出$id$的出现次数即可。

时间复杂度:$Theta(n)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
using namespace std;
int n,m,k;
int zls,cnt,id;
int c[1001],x[1001],y[1001],z[1001];
long long last;
int main()
{
	scanf("%d%d",&m,&k);
	for(int i=1;i<=m;i++){scanf("%d",&c[i]);n+=c[i];}
	for(int i=1;i<=m;i++)scanf("%d",&x[i]);
	for(int i=1;i<=m;i++)scanf("%d",&y[i]);
	for(int i=1;i<=m;i++)scanf("%d",&z[i]);
	int S=(1<<k)-1;
	for(int i=1;i<=m;i++)
	{
		last=x[i];
		if(!cnt||id==last){cnt++;id=last;}
		else cnt--;
		for(int j=1;j<c[i];j++)
		{
			last=(last*y[i]+z[i])&S;
			if(!cnt||id==last){cnt++;id=last;}
			else cnt--;
		}
	}
	cnt=0;
	for(int i=1;i<=m;i++)
	{
		last=x[i];
		if(last==id)cnt++;
		for(int j=1;j<c[i];j++)
		{
			last=(last*y[i]+z[i])&S;
			if(last==id)cnt++;
		}
	}
	if(n&1)zls=n>>1+1;
	else zls=n>>1;
	if(cnt>zls)printf("%d",cnt*2-n-1);
	else puts("0");
	return 0;
}

rp++

原文地址:https://www.cnblogs.com/wzc521/p/11566878.html