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+5;
int n,m,b[MAXN],ans=0;
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=1;i<=n;++i)
        scanf("%lld%lld",&a[i].l,&a[i].r);
    sort(a+1,a+n+1);
    for(int i=1;i<=m;++i){
        scanf("%lld",&b[i]);
        s.insert(b[i]);
    }
    for(int i=1;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
",ans);
    return 0;
}

玩具：神仙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]=sumlimits_{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=205;
int n,mod,inv[MAXN],ans=0,g[MAXN][MAXN],f[MAXN][MAXN];
signed main(){
    scanf("%lld%lld",&n,&mod);
    inv[0]=inv[1]=1;
    for(int i=2;i<=n;++i) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
    for(int i=0;i<=n;++i) g[0][i]=1;
    for(int i=1;i<=n;++i){
        for(int j=0;j<=n;++j){
            if(j) f[i][j]=g[i-1][j-1];
            else if(i==1)  f[i][j]=1;
            for(int k=1;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=1;i<=n;++i){
        ans=(ans+1ll*(f[n][i]-f[n][i-1]+mod)%mod*i%mod)%mod;
    }
    printf("%lld
",ans);
    return 0;
}

飘雪圣域：

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

#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
",r-l+1-query(root[l-1],root[r],1,n,l,r));
	}
	return 0;
}