2019.9.22
summary
- 回到了上半年的做题状态 知道难了 还在瓜想正解
- 该拿的分未拿到
- 日常瓜想
- 计算空间出错
==欧阳说各种情况都要经历一下 心态炸麻木了就好了
足球比赛
脑壳痛
(ans=sum C_n^i((pa*(1-qb))^i*(pb*qb)^{n*P-2})
还卡卡空间== 所以就预处理出fac阶乘的逆
#include<bits/stdc++.h>
#define ll long long
#define rg register
const int N=1e7+5,M=100+5,P=1e9+7,INF=1e9+7,inf=0x3f3f3f3f;
int n,p,q,pa,pb,qa,qb,facn,sum=0,ans=0,ifac[N];
template <class t>void rd(t &x){
x=0;int w=0;char ch=0;
while(!isdigit(ch)) w|=ch=='-',ch=getchar();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
x=w?-x:x;
}
int qpow(int a,int b){
int ret=1;
while(b){
if(b&1) ret=(ll)ret*a%P;
a=(ll)a*a%P,b>>=1;
}
return ret;
}
int main(){
#ifndef ONLINE_JUDGE
freopen("T1.txt","r",stdin);
//freopen("attack.out","w",stdout);
#endif
rd(n),rd(pa),rd(pb),rd(qa),rd(qb);
if(!pa||qa==qb) return puts("0"),0;
p=(ll)pa*qpow(pb,P-2)%P,q=(ll)qa*qpow(qb,P-2)%P;
facn=1;
for(rg int i=2;i<=n;++i) facn=(ll)facn*i%P;
ifac[n]=qpow(facn,P-2);
for(rg int i=n;i;--i) ifac[i-1]=(ll)ifac[i]*i%P;//i!的逆
int nwp=qpow((P+1-p)%P,n),nwq=qpow((P+1-q)%P,n);
int inv1=(ll)p*qpow((P+1-p)%P,P-2)%P;
int inv2=(ll)q*qpow((P+1-q)%P,P-2)%P;
for(int i=0;i<=(n>>1);++i)
ifac[n-i]=(ll)ifac[i]*ifac[n-i]%P,ifac[i]=ifac[n-i];
for(int i=1;i<=n;++i){
sum=(sum+(ll)nwq*ifac[i-1])%P,nwp=(ll)nwp*inv1%P;
ans=(ans+(ll)nwp*sum%P*ifac[i])%P,nwq=(ll)nwq*inv2%P;
}
ans=(ll)ans*facn%P;
if(p==1) ans=sum;
ans=(ll)ans*facn%P;
printf("%d
",ans);
return 0;
}
文明
其实很简单 考场上想出来了找两个点的中点来跳 可是并没有想到树剖然后线段树来维护其和
考场上打的暴力换根 再加上暴力上跳 结果没有想到跳到中点来 然后导致减重复了 其间还有很多低级错误(QAQ)
看到题解那个上跳到中点的操作之后就知道该怎么做了 一下子就打出来了
#include<bits/stdc++.h>
#define ll long long
#define ls (o<<1)
#define rs (o<<1|1)
using namespace std;
const int N=500000+5,M=100+5,P=7,INF=1e9+7,inf=0x3f3f3f3f;
int n,m,ans;
template <class t>void rd(t &x){
x=0;int w=0;char ch=0;
while(!isdigit(ch)) w|=ch=='-',ch=getchar();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
x=w?-x:x;
}
int head[N],tot=0;
struct edge{int v,nxt;}e[N<<1];
void add(int u,int v){e[++tot]=(edge){v,head[u]},head[u]=tot;}
int idx=0,dfn[N],id[N],son[N],f[N],dep[N],sz[N],top[N];
void dfs1(int u,int fa){
dep[u]=dep[fa]+1,f[u]=fa,sz[u]=1;
for(int i=head[u],v,mxs=-1;i;i=e[i].nxt)
if((v=e[i].v)!=fa){
dfs1(v,u),sz[u]+=sz[v];
if(mxs<sz[v]) mxs=sz[v],son[u]=v;
}
}
void dfs2(int u,int topf){
top[u]=topf,dfn[u]=++idx,id[idx]=u;
if(!son[u]) return;
dfs2(son[u],topf);
for(int i=head[u],v;i;i=e[i].nxt)
if((v=e[i].v)!=f[u]&&v!=son[u]) dfs2(v,v);
}
struct node{int sum,tg;}t[N<<2];
void pup(int o){t[o].sum=t[ls].sum+t[rs].sum;}
void updnode(int o,int l,int r,int k){t[o].sum=(r-l+1)*k,t[o].tg=k;}
void pudw(int o,int l,int r){
if(t[o].tg==-1) return;
int mid=l+r>>1;
updnode(ls,l,mid,t[o].tg),updnode(rs,mid+1,r,t[o].tg);
t[o].tg=-1;
}
void mdf(int o,int l,int r,int x,int y,int k){
if(l>y||r<x) return;
if(x<=l&&r<=y){updnode(o,l,r,k);return;}
pudw(o,l,r);
int mid=l+r>>1;
mdf(ls,l,mid,x,y,k),mdf(rs,mid+1,r,x,y,k);
pup(o);
}
int LCA(int x,int y){
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]]) swap(x,y);
x=f[top[x]];
}
return dep[x]<dep[y]?x:y;
}
int jump(int x,int k){
while(dep[x]-dep[top[x]]+1<=k){
k-=(dep[x]-dep[top[x]]+1);
x=f[top[x]];
}
return id[dfn[x]-k];
}
int main(){
freopen("T2.txt","r",stdin);
// freopen("civilization.out","w",stdout);
rd(n),rd(m);
for(int i=1,u,v;i<n;++i) rd(u),rd(v),add(u,v),add(v,u);
dfs1(1,0),dfs2(1,1);
for(int i=1;i<=(n<<2);++i) t[i].tg=-1;
while(m--){
int tt,rt;rd(tt),rd(rt);
mdf(1,1,n,1,n,1);
for(int i=2,x,lca,L,R,mid;i<=tt;++i){
rd(x);lca=LCA(rt,x);
L=dep[rt]-dep[lca]-1,R=dep[x]-dep[lca]-1,mid=L+R+1>>1;
if(mid<=R){
mid=jump(x,mid);
mdf(1,1,idx,dfn[mid],dfn[mid]+sz[mid]-1,0);
}
else{
mid=jump(rt,L+R+1-mid);
mdf(1,1,idx,1,dfn[mid]-1,0),
mdf(1,1,idx,dfn[mid]+sz[mid],idx,0);
}
}
printf("%d
",t[1].sum);
}
return 0;
}
贪玩蓝月
过于sb 因为低级错误调了一个上午==
过程过于艰辛我真的不想说
分成两个部分 分别维护两个背包 当一边空了时 就重构一下
用的ST表 当选了(f2[tt2][i]) 即左边w值之和与P取余为i时 然后再在右边找到一个对应的 使其w值之和对P取余时值在(lsim r)之间
#include<bits/stdc++.h>
#define ll long long
#define rg register
#define Max(x,y) ((x)>(y)?(x):(y))
using namespace std;
const int N=50000+5,M=500+5,INF=1e9+7,inf=0x3f3f3f3f;
int n,P,cnt=0;
ll st[M][25],val[N<<1],f1[N][M],f2[N][M];int lg[M];
char opt[5];
template <class t>void rd(t &x){
x=0;int w=0;char ch=0;
while(!isdigit(ch)) w|=ch=='-',ch=getchar();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
x=w?-x:x;
}
int q[N<<1],l,r,mid,tt1,tt2;
void add1(int w,ll v){
++tt1;
for(int i=0;i<P;++i) f1[tt1][i]=f1[tt1-1][i];
for(int i=0;i<P;++i)
f1[tt1][(i+w)%P]=max(f1[tt1][(i+w)%P],f1[tt1-1][i]+v);
}
void add2(int w,ll v){
++tt2;
for(int i=0;i<P;++i) f2[tt2][i]=f2[tt2-1][i];
for(int i=0;i<P;++i)
f2[tt2][(i+w)%P]=max(f2[tt2][(i+w)%P],f2[tt2-1][i]+v);
}
void rebuild(){
tt1=tt2=0,mid=l+r>>1;
for(int i=mid;i>=l;--i) add1(q[i],val[i]);
for(int i=mid+1;i<=r;++i) add2(q[i],val[i]);
}
int que(int x,int y){
int s=lg[y-x+1];
return Max(st[x][s],st[y-(1<<s)+1][s]);
}
ll query(int L,int R){
ll ret=-0x7fffffff;
for(int i=0;i<=P;++i) st[i][0]=f1[tt1][i];
for(int j=1;j<=20;++j)
for(int i=0;i+(1<<j)-1<=P;++i) st[i][j]=max(st[i][j-1],st[i+(1<<(j-1))][j-1]);
for(int i=0,x,y;i<P;++i){
if(f2[tt2][i]<0) continue;
x=(L-i+P)%P,y=(R-i+P)%P;
if(x<=y) ret=max(ret,f2[tt2][i]+que(x,y));
else ret=max(ret,f2[tt2][i]+Max(que(0,y),que(x,P-1)));
}
if(ret<0) return -1;
else return ret;
}
int main(){
#ifndef ONLINE_JUDGE
freopen("T1.txt","r",stdin);
//freopen("attack.out","w",stdout);
#endif
int T;rd(T);
rd(n),rd(P),lg[0]=-1,l=50001,r=mid=50000;
memset(f1[0],~inf,sizeof(f1[0]));
memset(f2[0],~inf,sizeof(f2[0]));
f1[0][0]=f2[0][0]=0,lg[0]=-1;
for(int i=1;i<=P;++i) lg[i]=lg[i>>1]+1;//,f1[0][i]=f2[0][i]=-inf
ll v;
for(int i=1,w,L,R;i<=n;++i){
scanf("%s",opt);
if(opt[0]=='I'){
rd(w),rd(v),w%=P;
if(opt[1]=='F') q[--l]=w,val[l]=v,add1(w,v);
else q[++r]=w,val[r]=v,add2(w,v);
}
else if(opt[0]=='D'){
if(opt[1]=='F'){
++l;
if(--tt1<=0) rebuild();
}
else{
--r;
if(--tt2<=0) rebuild();
}
}
else{
rd(L),rd(R),printf("%lld
",query(L,R));
// if(++cnt==15)
}
}
return 0;
}