欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

2021.03.06【NOIP提高B组】模拟 总结

程序员文章站 2024-03-18 21:58:52
...

T1

看起来十分复杂,打表后发现答案是 \(n*m\mod p\)

具体的证明。。。

原式的物理意义,就是从坐标原点(0,0),用每一种合法的斜率,

穿过坐标[1 ~ n , 1 ~ m]的方阵中的整点的个数,总数即 n*m。

T2

考试时没想到 \(O(n^2)\) 的做法

\[\begin{aligned} A_{k+i-1}-B_i &=A_{k+j-1}-B_j \\ A_{k+i-1}-A_{k+j-1} &= B_i-B_j \\ A_i-A_j &= B_i-B_j \end{aligned} \]

所以只要先差分,再跑一边 \(\text{KMP}\) 或字符串哈希即可

字符串哈希

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1000005;
const LL M1=998244353,M2=1004535809,B=233;
int n,m,ans;
LL X1[N],X2[N],P1[N],P2[N],Y1,Y2,x[N],y[N];
inline bool check(int i) {
	return
	(Y1+1LL*P1[m]*X1[i-1]%M1)%M1==X1[i+m-1] &&
	(Y2+1LL*P2[m]*X2[i-1]%M2)%M2==X2[i+m-1];
}
int main() {
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)scanf("%lld",&x[i]);
	for(int i=1;i<=m;i++)scanf("%lld",&y[i]);
	for(int i=1;i<n;i++)x[i]=x[i+1]-x[i],x[i]+=2e9;
	for(int i=1;i<m;i++)y[i]=y[i+1]-y[i],y[i]+=2e9;
	--n,--m;
	P1[0]=P2[0]=1;
	for(int i=1;i<=n;i++) {
		X1[i]=(1LL*X1[i-1]*B%M1+1LL*x[i])%M1;
		X2[i]=(1LL*X2[i-1]*B%M2+1LL*x[i])%M2;
		P1[i]=(1LL*P1[i-1]*B)%M1;
		P2[i]=(1LL*P2[i-1]*B)%M2;
	}
	for(int i=1;i<=m;i++) {
		Y1=(1LL*Y1*B%M1+1LL*y[i])%M1;
		Y2=(1LL*Y2*B%M2+1LL*y[i])%M2;
	}
	for(int i=1;i<=n-m+1;i++)
		if(check(i))++ans;
	printf("%d",ans);
}

\(\text{KMP}\)

#include<bits/stdc++.h>
using namespace std;
const int N=1000005;
int n,m,ans,x[N],y[N],nxt[N];
int main() {
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)scanf("%d",&x[i]);
	for(int i=1;i<=m;i++)scanf("%d",&y[i]);
	for(int i=1;i<n;i++)x[i]=x[i+1]-x[i];
	for(int i=1;i<m;i++)y[i]=y[i+1]-y[i];
	--n,--m;
	for(int i=2,j=0;i<=m;i++) {
		while(j && y[i]!=y[j+1])j=nxt[j];
		if(y[i]==y[j+1])++j;
		nxt[i]=j;
	}
	for(int i=0,j=0;i<=n;i++) {
		while(j && x[i]!=y[j+1])j=nxt[j];
		if(x[i]==y[j+1])++j;
		if(j==m)++ans,j=nxt[j];
	}
	printf("%d",ans);
}

T3

考试时:不考虑 \(k\) 直接树形 DP ,奇妙地拿了 30

赛后:由于原题意思就是不能直接地祖孙关系,因为一个子树地 DFS 序是连续的

只要选的区间不相交即可,用 DP 解决

\(dp_{i,j}\) 为第 \(j\) 个果子选了 DFS 序为 \(i\) 的节点的最大美味值

可以从 \(dp_{i+1,j}\)\(dp_{\text{next subtree},j}\) ,和 \(dp_{\text{next subtree},j-1}+w_{\text{node of i}}\) 转移

数据十分奇妙,要开滚动数组。

#include<bits/stdc++.h>
using namespace std;
const int N=100005;
int n,k,op,dp[N][2],f[N],w[N],to[N],lst[N],nxt[N],rk[N],sz[N],odd,ans;
void Pre(int u) {
	rk[++odd]=u,sz[u]=1;
	for(int i=lst[u],v;i;i=nxt[i])
		Pre(v=to[i]),sz[u]+=sz[v];
}
int main() {
	freopen("galo.in","r",stdin);
	freopen("galo.out","w",stdout);
	scanf("%d%d",&n,&k),++k;
	for(int i=2;i<=n;i++) {
		scanf("%d%d",&f[i],&w[i]);
		to[i]=i,nxt[i]=lst[f[i]];
		lst[f[i]]=i;
	}
	Pre(1);
	for(int C=1,j=0;C<=k;C++,j^=1) {
		for(int i=n;i>=1;i--)
			dp[i][j]=max(dp[i+1][j],max(dp[i+sz[rk[i]]][j],dp[i+sz[rk[i]]][j^1]+w[rk[i]]));
		ans=max(ans,dp[1][j]);
	}
	printf("%d",ans);
}

T4

好家伙,看到子树修改,本蒟蒻直接一个树剖,好吧看错题+爆栈直接 WA 10

正解:直接 BFS 即可,设 \(F_i\) 为从根到这个点有多少个打了标记,若 \(F_i+Color_i\) 是偶数就不用翻转

#include<bits/stdc++.h>
using namespace std;
inline int Rd() {
	register int x=0;
	char C=getchar();
	for(;C<'0'||C>'9';C=getchar()) ;
	for(;C>'/'&&C<':';C=getchar()) x=(x<<1)+(x<<3)+(C^48);
	return x;
}
const int N=500005;
int n,x[N],vis[N],lst[N],nxt[N],to[N],f[N],q[N],head,tail=1;
int main() {
	freopen("tree.in","r",stdin);
	freopen("tree.out","w",stdout);
	n=Rd();
	for(int i=1;i<=n;i++)x[i]=Rd();
	for(int i=1,fr;i<n;i++) {
		fr=Rd(),to[i]=Rd();
		nxt[i]=lst[fr],lst[fr]=i;		
	}
	q[1]=1;
	for(int u,flg;head<tail;) {
		u=q[++head],flg=(f[u]+x[u])&1;
		if(flg)vis[u]=1;
		for(int i=lst[u],v;i;i=nxt[i])
			f[v=to[i]]=f[u]+flg,q[++tail]=v;
	}
	for(int i=1;i<=n;i++)
		if(vis[i])printf("%d ",i);
}

总结

T1:学会打表

T2:多去把式子转换一下

T3:子树修改就想想 DFS 序,多考虑转换成 DP

T4:不一味的想数据结构,考虑简化