2021.03.06【NOIP提高B组】模拟 总结
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:不一味的想数据结构，考虑简化
