\(noip模擬27\;solutions\)

這次吧,我本來以為我能切掉兩個題,結果呢??只切掉了一個

不過,隔壁Varuxn也以為能切兩個,可惜了,他一個都沒切。。。。。。

確實他分比我高一點,但是吧,這個人就改題非常慢,所以結論就是

我牛逼,牛逼到家了

所以我應該是掛掉了100pts,下次注意,下次AK

\(T1\;妹子圖(graph)\)

這個題有一堆做法,而且旗鼓相當,複雜度都一樣,跑的也差不多快

我的做法是在無向圖上跑魔改拓撲排序+魔改DIJ,非常的噁心,但是吧,理解之後極其簡單

按照堆優化的DIJ搞了一個優先佇列,這樣我們先把半仙的家能夠連到的地方壓到佇列中

因為我們是按照邊權小的在前面,所以我們當前找到的一定是當前困難抵抗值能夠到達的最大的圖

顯然有一些困難抵抗值能夠見到的妹子種類數是一樣的,這些我們就不必儲存了

這樣的話,我們的數組裡最多隻有2e5個點,直接保證了時間和空間複雜度

我們開始彈佇列,然後更新這個陣列,這裡面存的東西就是

能夠看到當前這麼多妹子的最小困難抵抗值,和能夠看到的妹子

最後的時候直接在這個數組裡二分查詢,搞一搞邊界就好了

複雜度是\(O(nlogn+qlogn)\)

AC_code

#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=5e5+5;
int n,m,qq,hom,opt,mod;
int c[N];
ll num[N],fro[N],bas[N],cnt;
int to[N*2],nxt[N*2],val[N*2],head[N],rp;
void add_edg(int x,int y,int z){
to[++rp]=y;
val[rp]=z;
nxt[rp]=head[x];
head[x]=rp;
}
int ton[N],sum;
bool vis[N];
struct node{
int val,tag;
node(){}
node(int x,int y){
val=x;tag=y;
}
bool operator < (node x)const{
return val>x.val;
}
};
priority_queue<node> q;
signed main(){
scanf("%d%d%d%d%d",&n,&m,&qq,&hom,&opt);
if(opt)scanf("%d",&mod);
for(re i=1;i<=n;i++){
scanf("%d",&c[i]);
}
for(re i=1,x,y,z;i<=m;i++){
scanf("%d%d%d",&x,&y,&z);
add_edg(x,y,z);
add_edg(y,x,z);
}
ton[c[hom]]++;sum++;
for(re i=head[hom];i;i=nxt[i])
q.push(node(val[i],to[i]));
int cge=0;bas[++cnt]=1;num[cnt]=sum;
while(!q.empty()){
int dis=q.top().val;
int now=q.top().tag;
q.pop();
if(vis[now])continue;
vis[now]=true;
cge=max(cge,dis);
ton[c[now]]++;
if(ton[c[now]]==1){
sum++;
if(bas[cnt]!=cge)
bas[++cnt]=cge;
num[cnt]=sum;
}
for(re i=head[now];i;i=nxt[i]){
q.push(node(val[i],to[i]));
}
}
bas[++cnt]=1e9+1;
for(re i=1;i<=cnt;i++){
fro[i]=fro[i-1]+1ll*num[i-1]*(bas[i]-bas[i-1]);
//cout<<bas[i]<<" "<<num[i]<<" "<<fro[i]<<endl;
}
ll lans=0;
for(re i=1;i<=qq;i++){
ll x,y;
scanf("%lld%lld",&x,&y);
if(opt==1){
x=(x^lans)%mod+1;
y=(y^lans)%mod+1;
}
if(x>y)swap(x,y);
int ld=upper_bound(bas+1,bas+cnt+1,x)-bas;
int rd=upper_bound(bas+1,bas+cnt+1,y)-bas-1;
//cout<<ld<<" "<<rd<<" "<<num[ld]<<" "<<bas[ld]<<endl;
ll ans=fro[rd]-fro[ld];
ans+=1ll*(bas[ld]-x)*num[ld-1];
ans+=1ll*(y-bas[rd]+1)*num[rd];
printf("%lld\n",ans);
lans=ans;
}
}


\(T2\;妹子腿(tree)\)

真是沒想到啊,沒想到,這個竟然是個大暴力,主要是他還能過,氣死人,比線段樹還快

考場上沒敢寫,氣人。

就是維護被標記的點的標記時間以及和要查詢的點的距離,

這個距離就是樹剖+LCA,完事!!!

注意優化,要不然過不了,如果要標記的點已經在不斷的擴散中被標記了,就不要加進去了。

AC_code

#include<bits/stdc++.h>
using namespace std;
#define re register int
#define fi first
#define se second
#define pa pair<int,int>
#define mpa(x,y) make_pair(x,y)
const int N=1e5+5;
int n,m;
int to[N*2],nxt[N*2],head[N],rp;
void add_edg(int x,int y){
to[++rp]=y;
nxt[rp]=head[x];
head[x]=rp;
}
int fa[N],dep[N];
int siz[N],son[N],top[N];
void dfs1(int x){
siz[x]=1;son[x]=0;
for(re i=head[x];i;i=nxt[i]){
int y=to[i];
if(y==fa[x])continue;
fa[y]=x;dep[y]=dep[x]+1;
dfs1(y);
siz[x]+=siz[y];
if(!son[x]||siz[y]>siz[son[x]])son[x]=y;
}
}
void dfs2(int x,int f){
top[x]=f;
if(son[x])dfs2(son[x],f);
for(re i=head[x];i;i=nxt[i]){
int y=to[i];
if(y==fa[x]||y==son[x])continue;
dfs2(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;
}
int DIS(int x,int y){
return dep[x]+dep[y]-2*dep[LCA(x,y)];
}
pa ji[N];
int cnt;
signed main(){
scanf("%d%d",&n,&m);
for(re i=1,x,y;i<n;i++){
scanf("%d%d",&x,&y);
add_edg(x,y);add_edg(y,x);
}
dfs1(1);dfs2(1,1);
for(re i=1,typ,x;i<=m;i++){
scanf("%d%d",&typ,&x);
if(typ==1){
int flag=0;
for(re j=1;j<=cnt;j++){
if(DIS(ji[j].fi,x)<=i-ji[j].se){
flag=1;break;
}
}
if(flag==1)continue;
ji[++cnt]=mpa(x,i);
}
if(typ==2)cnt=0;
if(typ==3){
int flag=0;
for(re j=1;j<=cnt;j++){
if(DIS(ji[j].fi,x)<=i-ji[j].se){
flag=1;break;
}
}
if(flag==1)printf("wrxcsd\n");
else printf("orzFsYo\n");
}
}
}


\(T3\;妹子腰(序列)\)

這個題我一眼就是原題,而且是昨天的原題

後來發現不對,是極長上升而不是最長上升

所以我忘記了如何維護這個東西了,考完試又去回顧了一下,好難喲!!!

就是一個\(O(log^2n)\)的線段樹,注意維護一個底線,來保證是極長的。。

AC_code

#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=2e5+5;
const ll inf=0x3f3f3f3f3f3f3f3f;
const ll mod=998244353;
int n,a[N],maxd;
struct node{
#define ls x<<1
#define rs x<<1|1
ll sum[N*4];
int mxra[N*4];
ll work(int x,int l,int r,int maxn){
//cout<<x<<" "<<l<<" "<<r<<endl;
if(l==r){
if(mxra[x]>maxn)return sum[x];
return 0;
}
int mid=l+r>>1;
if(maxn>mxra[x])return 0;
else if(maxn>mxra[rs])return work(ls,l,mid,maxn);
else{
ll ret=0;
ret=(ret+work(rs,mid+1,r,maxn))%mod;
ret=(ret+work(ls,l,mid,mxra[rs]))%mod;
return ret;
}
}
void pushup(int x){
mxra[x]=max(mxra[ls],mxra[rs]);
//sum[x]=(sum[ls]+sum[rs])%mod;
return ;
}
void ins(int x,int l,int r,int pos,int mx,int v){
if(l==r){
mxra[x]=mx;
sum[x]=v;
return ;
}
int mid=l+r>>1;
if(pos<=mid)ins(ls,l,mid,pos,mx,v);
else ins(rs,mid+1,r,pos,mx,v);
pushup(x);return ;
}
ll query(int x,int l,int r,int pos){
if(pos==0)return 0;
if(r<=pos){
ll ret=work(x,l,r,maxd);
maxd=max(maxd,mxra[x]);
return ret;
}
int mid=l+r>>1;ll ret=0;
if(pos>mid)ret=(ret+query(rs,mid+1,r,pos))%mod;
ret=(ret+query(ls,l,mid,pos))%mod;
return ret;
}
#undef ls
#undef rs
}xds;
ll ans,dp[N];
signed main(){
scanf("%d",&n);
for(re i=1;i<=n;i++)scanf("%d",&a[i]);
for(re i=1;i<=n;i++){
maxd=0;
dp[i]=xds.query(1,1,n,a[i]-1);
//cout<<dp[i]<<" ";
if(dp[i]==0)dp[i]=1;
xds.ins(1,1,n,a[i],i,dp[i]);
//cout<<dp[i]<<endl;
}
int mx=0;
for(re i=n;i>=1;i--){
if(a[i]>mx)ans=(ans+dp[i])%mod,mx=a[i];
}
printf("%lld",ans);
}