洛谷NOIp热身赛题解
A 最大差值
简单树状数组,维护区间和、区间平方和,方差按照给的公式算就行了
#include
#define il inline
#define vd void
#define mod 1000000007
typedef long long ll;
namespace IO{
const int maxn=(1<<21)+1;
char ibuf[maxn],*iS,*iT,c;int f;
inline char getc(){
return iS==iT?(iT=(iS=ibuf)+fread(ibuf,1,maxn,stdin),iS==iT?EOF:*iS++):*iS++;
}
inline int gi(){
int x=0;
for(f=1,c=getc();(c<'0'||c>'9');c=getc())f=c=='-'?-1:1;
for(x=0;(c>='0'&&c<='9');c=getc())x=(x<<1)+(x<<3)+(c^48);
return x*f;
}
}
using IO::gi;
int b[100010];
int inv[100010];
int sum[100010<<4],SUM[100010<<4],N;
#define mid ((l+r)>>1)
il vd build(int n){
N=1<<(int)ceil(log(n+2)/log(2));
for(int i=1;i<=n;++i)sum[i+N]=b[i],SUM[i+N]=1ll*b[i]*b[i]%mod;
for(int i=N-1;i;--i)sum[i]=(sum[i<<1]+sum[i<<1|1])%mod,SUM[i]=(SUM[i<<1]+SUM[i<<1|1])%mod;
}
il vd update(int p,int d){
p+=N,sum[p]=d,SUM[p]=1ll*d*d%mod;
for(p>>=1;p;p>>=1)sum[p]=(sum[p<<1]+sum[p<<1|1])%mod,SUM[p]=(SUM[p<<1]+SUM[p<<1|1])%mod;
}
il std::pair query(int l,int r){
--l,++r;l+=N,r+=N;
int ret=0,RET=0;
while(l^r^1){
if(~l&1)ret=(sum[l^1]+ret)%mod,RET=(SUM[l^1]+RET)%mod;
if(r&1)ret=(sum[r^1]+ret)%mod,RET=(SUM[r^1]+RET)%mod;
l>>=1;r>>=1;
}
return std::make_pair(ret,RET);
}
int main(){
int n=gi(),m=gi();
inv[1]=1;for(int i=2;i<=n;++i)inv[i]=mod-1ll*(mod/i)*inv[mod%i]%mod;
for(int i=1;i<=n;++i)b[i]=gi();
build(n);
int o,l,r;
while(m--){
o=gi();
if(o==1)l=gi(),r=gi(),update(l,r);
else{
l=gi(),r=gi();
std::pairA=query(l,r);
int sa=A.first,SA=A.second;
int pa=1ll*sa*inv[r-l+1]%mod;
printf("%lld\n",1ll*inv[r-l+1]*(1ll*pa*pa%mod*(r-l+1)%mod+SA-1ll*2*sa*pa%mod+mod)%mod);
}
}
return 0;
}
B 攀爬者
略
C 蜈蚣
简单dp,设f[i][j]为1-j分成i段的最大收益,f[i][j]=max(f[i-1][k]+(A[k+1] xor ... xor A[j]))
#include
#define il inline
#define vd void
typedef long long ll;
il int gi(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int f[2][1010],W[1010];
il vd chkmx(int&a,int b){if(b>a)a=b;}
int main(){
int n=gi(),m=gi();
register int x=0,i,j,k;
for(i=1;i<=n;++i)W[i]=gi()^W[i-1];
memset(f[x]+1,-127,n<<2);
for(i=1;i<=m;++i){
x^=1;
memset(f[x],-127,i<<2);
for(j=i;j<=n;++j){
f[x][j]=0;
for(k=i-1;k<=j;++k)
chkmx(f[x][j],f[x^1][k-1]+(W[j]^W[k-1]));
}
}
printf("%d\n",f[x][n]);
return 0;
}
D 漂浮的鸭子
同消息传递,略
E 最大差值
略
F 随机数生成器
70分可以dp,设\(f[i]\)表示\(i\)到\(1\)的期望,那么\(f[1]=0,f[i]=1+\frac{\sum_{j=1}^{i}f[j]}{i}(i>1)\)
\(\frac{i-1}{i}f[i]=1+\frac{\sum_{j=1}^{i-1}f[j]}{i}\)
然后100分可以发现\(f[i]=1+(1/1+1/2+...+1/(i-1))\),后面的调和级数可以直接套公式
#include
#define il inline
#define vd void
typedef long long ll;
il int gi(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int main(){
int n=gi();if(n==1)return puts("0.00000"),0;
double ans=1;
if(n<=1000000)for(int i=1;i
G 大循环
题面看着很nb,实际上上下的a没有关系,F(q)是个定值,答案就是\(C_{n}^{k}\times F(q)\)
#include
#define il inline
#define vd void
#define mod 1000000007
#define int ll
typedef long long ll;
il ll gi(){
ll x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int p[500010];
il int pow(int x,int y){
int ret=1;
while(y){
if(y&1)ret=1ll*ret*x%mod;
x=1ll*x*x%mod;y>>=1;
}
return ret;
}
signed main(){
ll n=gi(),m=gi(),k=gi(),q=gi()%mod,F=0;
for(int x=1,i=0;i<=m;++i,x=1ll*x*q%mod)F=(F+1ll*x*gi()%mod)%mod;
p[0]=1;for(int i=1;i<=n;++i)p[i]=1ll*p[i-1]*i%mod;
printf("%lld\n",1ll*p[n]*pow(1ll*p[k]*p[n-k]%mod,mod-2)%mod*F%mod);
return 0;
}
H 会议座位
略
I 生日礼物
选两个数使得lcm=a,单独考虑a的每个质因数\(p^k\),那么这两个数一定有一个会有\(p^k\),另一个可以是\(p^{(0-k)}\),所以\(p^k\)有\(2k+1\)中选法。给a分解质因数,\(ans=\Pi(2k+1)\)
J HKE与他的小朋友
略
K 宝藏
二维树状数组 区间异或 区间查询异或和
和一维树状数组区间修改区间查询一样,记\(d(x,y)\)为\((x,y,n,n)\)这个子矩形的共同增量,查询\((x_0,y_0)\)只需查询\((1,1,x_0,y_0)\)的异或和。
冷静分析,\(d(x,y)\)在\((1,1,x_0,y_0)\)出现了奇数次才要计算\(d(x,y)\),条件就是\((x_0-x+1)(y_0-y+1)\)为奇数,那么\(x,x_0\)同奇偶,\(y,y_0\)同奇偶,对于每一种\(x,y\)的奇偶情况开一个树状数组即可,一共4个
#include
#define il inline
#define vd void
typedef long long ll;
il int gi(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
int n,m,N;
int t[2][2][1271][1271];
il vd update(int x,int y,int d){
int X=x&1,Y=y&1;x=(x+1)>>1,y=(y+1)>>1;
for(int i=x;i<=N;i+=i&-i)
for(int j=y;j<=N;j+=j&-j)
t[X][Y][i][j]^=d;
}
il int query(int x,int y){
int ret=0,X=x&1,Y=y&1;x=(x+1)>>1,y=(y+1)>>1;
for(int i=x;i;i-=i&-i)
for(int j=y;j;j-=j&-j)
ret^=t[X][Y][i][j];
return ret;
}
int main(){
n=gi(),m=gi();N=(n+1)>>1;
int x0,y0,x1,y1;char ch;
for(int i=1;i<=m;++i){
do ch=getchar();while(ch!='P'&&ch!='Q');
x0=gi(),y0=gi(),x1=gi(),y1=gi();
if(ch=='P'){
int t=gi(),S=0,a,b;
while(t--)a=gi(),b=gi(),S^=(b&1)<<(a-1);
update(x0,y0,S);
update(x0,y1+1,S);
update(x1+1,y0,S);
update(x1+1,y1+1,S);
}else{
int ans=query(x1,y1)^query(x1,y0-1)^query(x0-1,y1)^query(x0-1,y0-1);
for(int i=0;i<30;++i)putchar('1'+((ans>>i)&1));
puts("");
}
}
return 0;
}
L 简单的函数
回顾一下题目意思,对每个\(x\)找到最小的\(t\)满足\(t\not\mid x\),\(f(x)=f(t)+1\)
倒过来考虑,对每个\(t\)会转移到多少个\(x\)。(称为\(t\)的答案)\(x\)如果从\(t\)转移过来,那么\(lcm(1,2,...,t-1)\mid x,lcm(1,2,...,t)\not\mid x\)。
先考虑\(lcm(1,2,...,t-1)|x\),也就是说\(x\)是\(lcm(1,2,...,t-1)\)的倍数。满足这个的\(x\)有\(\lfloor\frac{n}{lcm(1,2,...,t-1)}\rfloor\)个。还要减去\(lcm(1,2,...,t)\mid x\)的情况,也就是所有\(t'>t\)的答案。
注意到\(lcm(1,...,42)\)大约是\(2*10^{17}\),那么可行的\(t\)不会超过\(42\),从\(42\)开始算即可。
#include
#define il inline
#define vd void
#define mod 1000000007
#define int ll
typedef long long ll;
il int gi(){
int x=0,f=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')f=-1;
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x*f;
}
il ll lcm(ll a,ll b){return a/std::__gcd(a,b)*b;}
ll f[]={0,0,1,2,3,2,4,2,3,2,3,2,3,2,3,2,3,2,4,2,3,2,3,2,3,2,3,2,3,2,4,2,3,2,3,2,3,2,3,2,3,2,4,2};
il int pow(int x,int y){
int ret=1;
while(y){
if(y&1)ret=1ll*ret*x%mod;
x=1ll*x*x%mod;y>>=1;
}
return ret;
}
ll L[50],a[50];
signed main(){
int n=gi();
L[1]=1;for(int i=2;i<=42;++i)L[i]=lcm(i,L[i-1]);
ll ans=1,sum=0;
for(int i=42;i;--i){
a[i]=n/L[i]-sum;
if(L[i]<=2)--a[i];
if(i==1)--a[i];
ans=1ll*ans*pow(f[i+1]+1,a[i])%mod;
sum+=a[i];
}
printf("%lld\n",ans);
return 0;
}
M 数列游戏
先dp一遍,\(f[l][r]\)表示\([l,r]\)这个区间里的数能否被消完,那么\(f[i][i+1]=[gcd(A[i],A[i+1])!=1]\)
\(f[l][r]\)有两种转移
- \(f[l+1][r-1] (gcd(A[l],A[r])!=1)\),先消完\([l+1,r-1]\),再消掉\(l,r\)
- \(f[l][k],f[k+1][r]\),消掉\([l,k]\)到\([k+1],r\)
这一段是\(O(n^3)\)的(我也不知道为啥能过)
然后就知道了哪些区间能用。剩下的看代码吧
#include
#define il inline
#define vd void
typedef long long ll;
il int gi(){
int x=0;
char ch=getchar();
while(!isdigit(ch)){
ch=getchar();
}
while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
return x;
}
ll A[1001],B[1001];
bool f[810][810];
ll g[810],sg[810];
il void chkmax(ll &a,ll b){if (a