2018.10.17多校联测测试总结

2018.10.17 多校联测测试总结

题目链接

T1梦境

解题思路:

题目意思是将点和区间作最大匹配,
先将区间按左端点和点一起从小到大排序
一路扫过去,类似于二维偏序,用堆维护插入点的右端点的最小值,然后查询

#pragma GCC optimize(2)
#include<queue>
#include<cstdio>
#include<algorithm>
using namespace std;
int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') ch=='-'&&(f=-1),ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int N=2e5+10;
int n,m,ans,tot;
priority_queue<int,vector<int>,greater<int> > que;
struct AC{int x,y;}a[N<<1];
bool cmp(AC x,AC y){return x.x==y.x?x.y>y.y:x.x<y.x;}
int main(){
    n=read(),m=read();
    for (int i=0;i<n;++i)a[++tot]=(AC){read(),read()};
    for (int i=0;i<m;++i)a[++tot]=(AC){read(),0};
    sort(a+1,a+n+m+1,cmp);
    for (int i=1,q=n+m;i<=q;++i){
        if (a[i].y) que.push(a[i].y);
        else {
            while ((!que.empty())&&a[i].x>que.top())que.pop();
            if (!que.empty()) que.pop(),++ans;
        }
    }
    printf("%d\n",ans);
    return 0;
}

T2 玩具

解题思路:

预处理\(dp[i][j]\)示i个点的森林,有j个点在第⼀棵树的概率,转移考虑第i个点是否在第⼀棵子树中,我们有状态转移方程
\[dp[i][j] = dp[i − 1][j − 1] ∗ (j − 1) ∗ inv[i] + dp[i − 1][j] ∗ (i − j) ∗ inv[i]\]
考虑修改算法三中状态的含义,令\(f[i][j]\)表示有i个点的树,深度不超过j的概率,\(g[i][j]\)表示有i个点的森林,深度不超过j的概率,\(f[i][j]\)直接从\(g[i-1][j-1]\)转移来;\(g[i][j]\)考虑枚举第⼀棵树的大小k,从⼀棵树和⼀个森林转移来,同时还要乘上第⼀棵子树大小为k的概率,我们有状态转移方程:
\[g[i][j] = ∑ f[k][j] ∗ g[i − k][j] ∗ dp[i][k]\]
k=1具体的细节可以见标程。最后只要用\(f[n][j]-f[n][j-1]\)就可以得到深度为j的树的概率
时间复杂度: \(O(n^3)\)

#pragma GCC optimize(2)
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
typedef long long ll;
const int N=210;
int n,P,ans;
int inv[N],dp[N][N],f[N][N],g[N][N],get[N];
int ksm(int x,int k){
    int su;
    for (su=1;k;x=1ll*x*x%P,k>>=1) 
        if (k&1) su=1ll*su*x%P;
    return su;
}
int solve1(int i,int j);
int solve2(int i,int j){
    if (j<0) return 0;
    if (!i) return 1;
    if (~g[i][j]) return g[i][j];
    g[i][j]=0;
    for (int k=1;k<=i;++k)
        (g[i][j]+=1ll*solve1(k,j)*solve2(i-k,j)%P*dp[i][k]%P)%=P;
    return g[i][j];
}
int solve1(int i,int j){
    if (~f[i][j]) return f[i][j];
    return f[i][j]=solve2(i-1,j-1);
}
int main(){
    cin>>n>>P;
    for (int i=1;i<=n;++i) inv[i]=ksm(i,P-2);
    dp[1][1]=1;
    for (int i=2;i<=n;++i)
        for (int j=1;j<=i;++j)
            dp[i][j]=(1ll*dp[i-1][j-1]*(j-1)%P*inv[i]+1ll*dp[i-1][j]*(i-j)%P*inv[i])%P;
    memset(f,-1,sizeof f);
    memset(g,-1,sizeof g);
    for (int i=2;i<=n;++i)get[i]=solve1(n,i),ans=(ans+1ll*(get[i]-get[i-1]+P)*i)%P;
    cout<<(ans-1+P)%P<<endl;
    return 0;
}

T3(飘雪圣域)

解题思路:

对于每次询问的\(x_i,y_i\),所有左右端点位于这个区间的边都会将两个联通块连在一起(因为这是一棵树,加一条边少一个联通块),所以此题就是对于每个询问作二维偏序,此题没有强制在线,所以主席树和树状数组都可写。

主席树:

#pragma GCC optimize(2)
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') ch=='-'&&(f=-1),ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int N=2e5+10,M=4e6;
int n,m,tot;
int rt[N],ls[M],rs[M],t[M];
struct AC{int x,y;}a[N];
bool cmp(AC x,AC y){return x.x<y.x;}
void updata(int p){t[p]=t[ls[p]]+t[rs[p]];}
void change(int&p,int l,int r,int x,int las){
    if (p==las) p=++tot,t[p]=t[las],ls[p]=ls[las],rs[p]=rs[las];
    if (l==r) return ++t[p],void();
    int mid=(l+r)>>1;
    if (x<=mid) change(ls[p],l,mid,x,ls[las]);
    else change(rs[p],mid+1,r,x,rs[las]);
    updata(p);
}
int query(int p1,int p2,int l,int r,int x,int y){
    if (x<=l&&r<=y)return t[p2]-t[p1];
    int mid=(l+r)>>1;
    return (x<=mid?query(ls[p1],ls[p2],l,mid,x,y):0)+(y>mid?query(rs[p1],rs[p2],mid+1,r,x,y):0);
}
int main(){
    n=read(),m=read();
    for (int i=1;i<n;++i)
        a[i]=(AC){read(),read()};
    sort(a+1,a+n,cmp);
    for (int i=1,j=1;i<=n;++i){
        rt[i]=rt[i-1];
        while (a[j].x==i&&j<n) change(rt[i],1,n,a[j].y,rt[i-1]),++j;//printf("...%d %d\n",a[j].x,a[j].y),++j;
    }
    for (int i=0,x,y;i<m;++i){
        x=read(),y=read();
        printf("%d\n",y-x+1-query(rt[x-1],rt[y],1,n,x,y));
    }
    return 0;
}

树状数组

#pragma GCC optimize(2)
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define low(x) (x&(-x))
using namespace std;
typedef long long ll;
int read(){
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9') ch=='-'&&(f=-1),ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
const int N=2e5+10;
int n,m,tot;
int c[N],ans[N];
struct AC{int x,y,z,id;}a[N<<1];
bool cmp(AC x,AC y){return x.x==y.x?x.z<y.z:x.x>y.x;}
void add(int x,int y){while (x<=n) c[x]+=y,x+=low(x);}
int query(int x){
    int su=0;
    while (x) su+=c[x],x-=low(x);
    return su;
}
int main(){
    n=read(),m=read();
    for (int i=1,x,y;i<n;++i)x=read(),y=read(),x>y?swap(x,y):void(),a[++tot]=(AC){x,y,0,i};
    for (int i=0,x,y;i<m;++i)x=read(),y=read(),x>y?swap(x,y):void(),a[++tot]=(AC){x,y,1,i};
    sort(a+1,a+n+m,cmp);
    for (int i=1,q=n+m;i<q;++i){
        if (a[i].z) ans[a[i].id]=a[i].y-a[i].x+1-query(a[i].y);
        else add(a[i].y,1);
    }
    for (int i=0;i<m;++i)printf("%d\n",ans[i]);
    return 0;
}
posted @ 2018-10-19 17:28 蒟蒻czx 阅读(...) 评论(...) 编辑 收藏