csps模拟测试74梦境，玩具，飘雪圣域题解

题面：https://www.cnblogs.com/Juve/articles/11679226.html

梦境：

其实还是挺水的，排序错了过不了样例，打了个二分图匹配就跑了

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<set>
 #define int long long
 using namespace std;
 const int MAXN=2e5+;
 int n,m,b[MAXN],ans=;
 struct node{
     int l,r;
     friend bool operator < (node p,node q){
         return p.r==q.r?p.l<q.l:p.r<q.r;
     }
 }a[MAXN];
 multiset<int>s;
 signed main(){
     scanf("%lld%lld",&n,&m);
     for(int i=;i<=n;++i)
         scanf("%lld%lld",&a[i].l,&a[i].r);
     sort(a+,a+n+);
     for(int i=;i<=m;++i){
         scanf("%lld",&b[i]);
         s.insert(b[i]);
     }
     for(int i=;i<=n;++i){
         multiset<int>::iterator it=s.lower_bound(a[i].l);
         if(it!=s.end()&&(*it)<=a[i].r){
             ++ans;
             s.erase(it);
         }
     }
     printf("%lld\n",ans);
     return ;
 }

玩具：神仙dp，

颓的题解和方程

令f[i][j]表示有i个点的树,深度不超过j的概率,g[i][j]表示有i个点的森林,深度不超过j的概率,

f[i][j]=g[i-1][j-1](j!=0) f[0][1]=1;

$g[i][j]=\sum\limits_{k=1}^{i}f[k][j]*g[i-k][j]*inv[i]$

 #include<iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #define int long long
 using namespace std;
 const int MAXN=;
 int n,mod,inv[MAXN],ans=,g[MAXN][MAXN],f[MAXN][MAXN];
 signed main(){
     scanf("%lld%lld",&n,&mod);
     inv[]=inv[]=;
     for(int i=;i<=n;++i) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
     for(int i=;i<=n;++i) g[][i]=;
     for(int i=;i<=n;++i){
         for(int j=;j<=n;++j){
             if(j) f[i][j]=g[i-][j-];
             else if(i==)  f[i][j]=;
             for(int k=;k<=i;++k)
                   g[i][j]=(g[i][j]+1ll*f[k][j]*g[i-k][j]%mod*inv[i]%mod)%mod;
         }
     }
     for(int i=;i<=n;++i){
         ans=(ans+1ll*(f[n][i]-f[n][i-]+mod)%mod*i%mod)%mod;
     }
     printf("%lld\n",ans);
     return ;
 }

飘雪圣域：

主席树，树状数组，莫队都可以过

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#define re register
using namespace std;
const int MAXN=2e5+5;
int n,q;
vector<int>mp[MAXN];
int tot=0,root[MAXN];
struct node{
	int ls,rs,val;
}tr[MAXN*30];
void build(int &k,int l,int r){
	k=++tot;
	if(l==r) return ;
	int mid=(l+r)>>1;
	build(tr[k].ls,l,mid);
	build(tr[k].rs,mid+1,r);
}
void insert(int &now,int pre,int l,int r,int pos){
	now=++tot;
	tr[now]=tr[pre],++tr[now].val;
	if(l==r) return ;
	int mid=(l+r)>>1;
	if(pos<=mid) insert(tr[now].ls,tr[pre].ls,l,mid,pos);
	else insert(tr[now].rs,tr[pre].rs,mid+1,r,pos);
}
int query(int x,int y,int l,int r,int opl,int opr){
	if(opl<=l&&r<=opr) return tr[y].val-tr[x].val;
	int mid=(l+r)>>1,res=0;
	if(opl<=mid) res+=query(tr[x].ls,tr[y].ls,l,mid,opl,opr);
	if(opr>mid) res+=query(tr[x].rs,tr[y].rs,mid+1,r,opl,opr);
	return res;
}
signed main(){
	scanf("%d%d",&n,&q);
	for(int i=1,u,v;i<n;++i){
		scanf("%d%d",&u,&v);
		mp[min(u,v)].push_back(max(u,v));
	}
	for(int i=1;i<=n;++i){
		root[i]=root[i-1];
		int N=mp[i].size();
		for(int j=0;j<N;++j){
			insert(root[i],root[i],1,n,mp[i][j]);
		}
	}
	while(q--){
		int l,r;
		scanf("%d%d",&l,&r);
		printf("%d\n",r-l+1-query(root[l-1],root[r],1,n,l,r));
	}
	return 0;
}