noip模擬34 solutions
我從來不為失敗找藉口，因為敗了就是敗了，沒人聽你訴說任何事情
今天很傷感，以來考試沒考好，二來改題改半天也改不出來
這次算是炸出來了我經常範的一些錯誤，比如除以0
一共只有25pts，是我模擬考試以來最最最慘的成績了吧
T1 Merchant
好像考場上的時候沒有太想好思路就開始打了。。。。。。
二分這個不用說，一眼就是二分。。。。
一開始碼了一個01揹包用來判斷，複雜度\(\mathcal{O(n^2logn)}\)
後來發現不對，我好像直接取前m個最大的就好了，複雜度\(\mathcal{O(nlog^2n)}\)
這個就是我考場上的22pts程式碼，為啥，因為我除以0，直接段錯誤了
考場下，我直接想起來 nth_element ，這玩意可以\(\mathcal{O(n)}\)求前m大的數啊
然後我直接把sort換成它，WA 44pts，為啥，因為我沒有將負數賦值為0
如果前m個數中有負數，我完全可以不加他，所以要和0取max
AC_code
#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=1e6+5;
ll n,m,s,maxn;
ll k[N],b[N],no[N];
ll f[N],bb[N];
bool check(ll x){
for(re i=1;i<=n;i++){
no[i]=k[i]*x+b[i];
no[i]=max(no[i],0ll);
}
nth_element(no+1,no+n-m+1,no+n+1);
ll tmp=0;
for(re i=n;i>=n-m+1;i--){
tmp+=no[i];
if(tmp>=s)return true;
}
return false;
}
signed main(){
scanf("%lld%lld%lld",&n,&m,&s);
bool flag=false;
for(re i=1;i<=n;i++){
scanf("%lld%lld",&k[i],&b[i]);
bb[i]=b[i];
}
sort(bb+1,bb+n+1);
ll tmp=0;
for(re i=n;i>=n-m+1;i--){
tmp+=bb[i];
if(tmp>=s)flag=true;
}
ll l=0,r=1000000000ll,mid;
while(l<r){
mid=l+r>>1;
if(check(mid))r=mid;
else l=mid+1;
}
printf("%lld",r);
}
T2 Equation
這個吧我考場上想到了一種線段樹+樹剖的做法。
這個做法吧是通過兩個式子的加減實現的，記錄非常的麻煩，而且還帶兩個log
碼長270，為了祭奠如此之長的程式碼，我特地把他粘在這裡
0pts
#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=1e6+5;
int n,q;
int to[N],nxt[N],head[N],rp;
void add_edg(int x,int y){
to[++rp]=y;
nxt[rp]=head[x];
head[x]=rp;
}
ll val[N],len[N];
int siz[N],son[N],fa[N],dep[N],top[N];
int dfn[N],cnt,idf[N];
void dfs_fi(int x){
siz[x]=1;son[x]=0;
for(re i=head[x];i;i=nxt[i]){
int y=to[i];
dep[y]=dep[x]+1;
dfs_fi(y);
siz[x]+=siz[y];
if(!son[x]||siz[y]>siz[son[x]])son[x]=y;
}
}
void dfs_se(int x,int f){
top[x]=f;dfn[x]=++cnt;idf[cnt]=x;
if(son[x])dfs_se(son[x],f);
for(re i=head[x];i;i=nxt[i]){
int y=to[i];
if(y==son[x])continue;
dfs_se(y,y);
}
}
int LCA(int x,int y){
while(top[x]!=top[y]){
if(dep[top[x]]<dep[top[y]])swap(x,y);
x=fa[top[x]];
}
return dep[x]<dep[y]?x:y;
}
struct XDS{
#define ls x<<1
#define rs x<<1|1
ll tr[N];
bool typ[N],tyq[N];
bool tp[N],tq[N];
bool rep,req;
void pushup(int x){
if(tyq[ls]&&typ[rs]){
typ[x]=typ[ls];
tyq[x]=!tyq[rs];
tr[x]=tr[ls]-tr[rs];
return ;
}
if(!tyq[ls]&&!typ[rs]){
typ[x]=!typ[ls];
tyq[x]=tyq[rs];
tr[x]=tr[rs]-tr[ls];
return ;
}
if(tyq[ls]&&!typ[rs]){
typ[x]=typ[ls];
tyq[x]=tyq[rs];
tr[x]=tr[ls]+tr[rs];
return ;
}
if(!tyq[ls]&&typ[rs]){
typ[x]=typ[ls];
tyq[x]=tyq[rs];
tr[x]=tr[ls]+tr[rs];
return ;
}
}
void build(int x,int l,int r){
if(l==r){
typ[x]=tyq[x]=true;
tr[x]=val[idf[l]];
return ;
}
int mid=l+r>>1;
build(ls,l,mid);
build(rs,mid+1,r);
pushup(x);
return ;
}
void ins(int x,int l,int r,int pos,ll v){
if(l==r){
typ[x]=tyq[x]=true;
tr[x]=v;
return ;
}
int mid=l+r>>1;
if(pos<=mid)ins(ls,l,mid,pos,v);
else ins(rs,mid+1,r,pos,v);
pushup(x);
return ;
}
ll query(int x,int l,int r,int ql,int qr){
if(ql>qr)return 0;
if(ql<=l&&r<=qr){
tp[x]=typ[x];
tq[x]=tyq[x];
return tr[x];
}
int mid=l+r>>1;
ll tmp1=0,tmp2=0;
if(ql<=mid)tmp1=query(ls,l,mid,ql,qr);
if(qr>mid)tmp2=query(rs,mid+1,r,ql,qr);
if(!tmp1){
tp[x]=tp[rs];tq[x]=tq[rs];
return tmp2;
}
if(!tmp2){
tp[x]=tp[ls];tq[x]=tq[ls];
return tmp1;
}
if(tq[ls]&&tp[rs]){
tp[x]=tp[ls];tq[x]=!tq[rs];
return tmp1-tmp2;
}
if(!tq[ls]&&!tp[rs]){
tp[x]=!tp[ls];tq[x]=tq[rs];
return tmp2-tmp1;
}
tp[x]=tp[ls];tq[x]=tq[rs];
return tmp1+tmp2;
}
#undef ls
#undef rs
}xds;
bool rxp,rxq,ryp,ryq;
ll get_this(int x,int y){
rxp=false;rxq=true;
ryp=false;ryq=true;
ll rex=0,rey=0;
while(top[x]!=top[y]){
if(dep[top[x]]>dep[top[y]]){
int tmp=xds.query(1,1,n,dfn[top[x]],dfn[x]);
if(xds.tq[1]&&rxp){
rxp=xds.tp[1];rxq=!rxq;
rex=tmp-rex;
}
else if(!xds.tq[1]&&!rxp){
rxp=!xds.tp[1];rxq=rxq;
rex=rex-tmp;
}
else {
rxp=xds.tp[1];rxq=rxq;
rex+=tmp;
}
x=fa[top[x]];
}
else{
int tmp=xds.query(1,1,n,dfn[top[y]],dfn[y]);
if(xds.tq[1]&&ryp){
ryp=xds.tp[1];ryq=!ryq;
rey=tmp-rey;
}
else if(!xds.tq[1]&&!ryp){
ryp=!xds.tp[1];ryq=ryq;
rey=rey-tmp;
}
else {
ryp=xds.tp[1];ryq=ryq;
rey+=tmp;
}
y=fa[top[y]];
}
}
if(dep[x]>dep[y]){
int tmp=xds.query(1,1,n,dfn[y]+1,dfn[x]);
if(xds.tq[1]&&rxp){
rxp=xds.tp[1];rxq=!rxq;
rex=tmp-rex;
}
else if(!xds.tq[1]&&!rxp){
rxp=!xds.tp[1];rxq=rxq;
rex=rex-tmp;
}
else {
rxp=xds.tp[1];rxq=rxq;
rex+=tmp;
}
}
else{
int tmp=xds.query(1,1,n,dfn[x]+1,dfn[y]);
//cout<<tmp<<endl;
if(xds.tq[1]&&ryp){
ryp=xds.tp[1];ryq=!ryq;
rey=tmp-rey;
}
else if(!xds.tq[1]&&!ryp){
ryp=!xds.tp[1];ryq=ryq;
rey=rey-tmp;
}
else {
ryp=xds.tp[1];ryq=ryq;
rey+=tmp;
}
}
if(rxp&&ryp){
ryq=!ryq;
return rex-rey;
}
if(!rxp&&!ryp){
rxq=!rxq;
return rey-rex;
}
return rex+rey;
}
signed main(){
scanf("%d%d",&n,&q);
for(re i=2;i<=n;i++){
int x;scanf("%d%lld",&x,&val[i]);
add_edg(x,i);fa[i]=x;
}
dfs_fi(1);dfs_se(1,1);
xds.build(1,1,n);
while(q--){
int ty;
scanf("%d",&ty);
if(ty==1){
int u,v;ll s;
scanf("%d%d%lld",&u,&v,&s);
if(u==v){
ll bi=s/2;
ll now=get_this(1,v);
if(1==v)printf("none\n");
else if(rxq&&ryq)printf("%lld\n",now-bi);
else if(!rxq&&!ryq)printf("%lld\n",-now-bi);
else if(!rxq&&ryq)printf("%lld\n",now+bi);
else printf("%lld\n",-now+bi);
continue;
}
ll tmp=get_this(u,v);
//cout<<dfn[u]<<" "<<dfn[v]<<endl;
//cout<<rxq<<" "<<ryq<<" "<<tmp<<endl;
if(rxq&&ryq&&s==tmp)
printf("inf\n");
else if(rxq&&ryq&&s!=tmp)
printf("none\n");
else {
ll now;
if(rxq&&!ryq){
ll bi=(s+tmp)/2;
now=get_this(1,u);
if(rxq&&ryq)printf("%lld\n",now-bi);
else if(!rxq&&!ryq)printf("%lld\n",-now-bi);
else if(!rxq&&ryq)printf("%lld\n",now+bi);
else printf("%lld\n",-now+bi);
}
else{
ll bi=(s+tmp)/2;
now=get_this(1,v);
if(rxq&&ryq)printf("%lld\n",now-bi);
else if(!rxq&&!ryq)printf("%lld\n",-now-bi);
else if(!rxq&&ryq)printf("%lld\n",now+bi);
else printf("%lld\n",-now+bi);
}
}
}
else{
int u;ll w;
scanf("%d%lld",&u,&w);
xds.ins(1,1,n,dfn[u],w);
}
}
}
但是後來發現其實沒有那麼複雜，直接利用邊權和x1來表示每個點的點權
我們只需要通過奇偶來判斷應該加還是減
直接樹狀陣列維護就好了
AC_code
#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=1e6+5;
int n,q,fa[N];
ll val[N];
int to[N],nxt[N],head[N],rp;
void add_edg(int x,int y){
to[++rp]=y;
nxt[rp]=head[x];
head[x]=rp;
}
int dfn[N],dfm[N],cnt,idf[N];
int dep[N];
void dfs(int x){
dfn[x]=++cnt;
idf[cnt]=x;
for(re i=head[x];i;i=nxt[i]){
int y=to[i];
dep[y]=dep[x]+1;
dfs(y);
}
dfm[x]=cnt;
}
ll tr[N];
int lb(int x){return x&(-x);}
void ins(int x,ll v){
for(re i=x;i<=n;i+=lb(i))
tr[i]+=v;
}
ll query(int x){
ll ret=0;
for(re i=x;i;i-=lb(i))
ret+=tr[i];
return ret;
}
signed main(){
scanf("%d%d",&n,&q);
for(re i=2;i<=n;i++){
scanf("%d%lld",&fa[i],&val[i]);
add_edg(fa[i],i);
}
dfs(1);
for(re i=2;i<=n;i++){
int tmp;
if(dep[i]%2)tmp=-1;
else tmp=1;
ins(dfn[i],val[i]*tmp);
ins(dfm[i]+1,-val[i]*tmp);
}
while(q--){
int typ;
scanf("%d",&typ);
if(typ==1){
int u,v;
ll s;
scanf("%d%d%lld",&u,&v,&s);
ll tmpu=query(dfn[u]);
ll tmpv=query(dfn[v]);
if(dep[u]%2)tmpu=-tmpu;
if(dep[v]%2)tmpv=-tmpv;
if(dep[u]%2==dep[v]%2){
ll tmp=tmpu+tmpv-s;
if(tmp%2)printf("none\n");
else if(dep[u]%2)printf("%lld\n",tmp/2);
else if(dep[v]%2==0)printf("%lld\n",-tmp/2);
}
else{
ll tmp=tmpu+tmpv;
if(tmp==s)printf("inf\n");
else printf("none\n");
}
}
else{
int u,tmp;
ll w;
scanf("%d%lld",&u,&w);
if(dep[u]%2)tmp=-1;
else tmp=1;
ins(dfn[u],-val[u]*tmp);
ins(dfm[u]+1,val[u]*tmp);
val[u]=w;
ins(dfn[u],val[u]*tmp);
ins(dfm[u]+1,-val[u]*tmp);
}
}
}
T3 Rectangle
所以我因為這個題錯過了A層？？？？
其實挺傷心的，自己改題總是特別慢，而且思路還總是有偏差
就這個題一開始我想偏了好多，
我用當前的樹狀陣列去更新前面的，而題解是直接將前面的賦值過來
我也是吐了啊，真的很是無奈，最後迫不得已只能去看程式碼
總而言之就是程式碼能力過於弱了
這個就是一個用樹狀陣列優化的列舉題；
你發現只有在邊界上有點的時候這個正方形才是合法的，
那麼我們假裝每一列只有一個點，那麼我們就可以固定這個點
從這個點向前列舉前面每一列的點，我們如果想用這兩列作為矩形的邊界
這兩個點必須要選上，我們只需要找到那些y特別大特別小的就好了
直接把每個矩形都組合出來
那麼如果擴充套件到多個點的時候，我們只需要根據這兩列的點吧整個序列分成幾個塊就好了
AC_code
#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=1e4+5;
const int M=2505;
const ll mod=1e9+7;
int n;
int sca[M][M];
bool vis[M][M];
ll ans;
struct BIT{
ll tr[M];
int lb(int x){return x&(-x);}
void ins(int x,ll v){
for(re i=x;i<=2500;i+=lb(i))
tr[i]+=v;
}
ll query(int x){
ll ret=0;
for(re i=x;i;i-=lb(i))
ret+=tr[i];
return ret;
}
}bit[M],num[M];
signed main(){
scanf("%d",&n);
for(re i=1;i<=n;i++){
int x,y;
scanf("%d%d",&x,&y);
sca[x][++sca[x][0]]=y;
}
for(re i=1;i<=2500;i++){
sca[i][sca[i][0]+1]=2501;
sort(sca[i]+1,sca[i]+sca[i][0]+1);
}
for(re i=1;i<=2500;i++){
if(!sca[i][0])continue;
for(re j=1;j<=sca[i][0];j++){
vis[i][sca[i][j]]=true;
num[i].ins(sca[i][j],1);
bit[i].ins(sca[i][j],sca[i][j]);
}
for(re j=i-1;j>=1;j--){
if(!sca[j][0])continue;
for(re k=1;k<=sca[j][0];k++){
if(vis[i][sca[j][k]])continue;
vis[i][sca[j][k]]=true;
num[i].ins(sca[j][k],1);
bit[i].ins(sca[j][k],sca[j][k]);
}
int noi=1,noj=1;
int upy=max(sca[i][1],sca[j][1]);
while(sca[i][noi+1]<=upy)noi++;
while(sca[j][noj+1]<=upy)noj++;
while(noi<=sca[i][0]&&noj<=sca[j][0]){
int mx=min(sca[i][noi+1],sca[j][noj+1]);
int mn=min(sca[i][noi],sca[j][noj]);
ll uy=bit[i].query(mx-1)+mod-bit[i].query(upy-1);
ll dy=bit[i].query(mn);
ll us=num[i].query(mx-1)+mod-num[i].query(upy-1);
ll ds=num[i].query(mn);
ans=(ans+(uy*ds%mod+mod-dy*us%mod)%mod*(i-j)%mod)%mod;
upy=mx;
if(sca[i][noi+1]<=mx)noi++;
if(sca[j][noj+1]<=mx)noj++;
}
}
}
printf("%lld",ans);
}