牛客网暑期ACM多校训练营（第三场

于 2018-07-30 01:35:19 发布
阅读量157 收藏
点赞数
分类专栏： 比赛总结 文章标签： 牛客 多校
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_34921856/article/details/81278408
版权

B Expected Number of Nodes
对于求期望这一类问题，我们通常可以算每个点对答案的贡献，然后累加。
假设当前要选择k个点。
对于度数小于等于2的点，它能存在的情况只有它被选进去，所以它能存在的概率即为 C(n−1,k−1)/C(n,k) C ( n − 1 , k − 1 ) / C ( n , k ) 。可以发现它们的概率的分母都为 C(n,k) C ( n , k ) ，所以我们下面只关注它的分子的值，最后乘上分母的逆元就行了。
对于度数大于2的点，当且仅当它有不少于3个孩子的子树被选了，它才会留下来。
所以这里倒过来算，所有的情况减掉它不会留下来的情况即为它的贡献。
设它有m个孩子，每个孩子的大小为sz[i]，则它不能留下来的情况为
∑i<jC(sz[i]+sz[j],k)−(m−2)∑iC(sz[i],k) ∑ i < j C ( s z [ i ] + s z [ j ] , k ) − ( m − 2 ) ∑ i C ( s z [ i ] , k )
我们在遍历j时，k个点都在i子树里面被算了m-1次，我们只需要1次，所以要减去(m-2)次。
我们可以发现，在遍历k时，组合数C(a,k)中的a不会变，所以我们只需要提前维护每个a出现的次数，而非每次都得dfs一遍。

#include<bits/stdc++.h>
using namespace std;
const int maxn=5005;
const int mod=1e9+7;
long long val[maxn];
int k;
vector<int>V[maxn];
int Size[maxn];
int cnt[maxn];
long long extgcd(long long a,long long b,long long &x,long long &y)
{
    if(a==0&&b==0)return -1;
    if(b==0){x=1;y=0;return a;}
    long long d=extgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
long long mod_reverse(long long a,long long n)
{
    long long x,y;
    long long d=extgcd(a,n,x,y);
    if(d==1)return (x%n+n)%n;
    return -1;
}
long long tab1[maxn],tab2[maxn];
void init()
{
    tab1[0]=1,tab2[0]=mod_reverse(1,mod);
    for(int i=1;i<=5000;i++)
        tab1[i]=1LL*i*tab1[i-1]%mod,tab2[i]=mod_reverse(tab1[i],mod);
}

long long C(int a,int b)
{
    if(b>a)return 0;
    return tab1[a]*tab2[b]%mod*tab2[a-b]%mod;
}
int N;
int sigle;
void dfs(int u,int pre)
{
    Size[u]=1;
    if(V[u].size()<=2)
    {
        sigle++;
        //val[u]=C(N-1,k-1);
        for(int v:V[u])if(v!=pre)
                dfs(v,u),Size[u]+=Size[v];
    }
    else
    {
        for(int v:V[u])if(v!=pre)
                dfs(v,u),Size[u]+=Size[v];
        //val[u]=C(N,k);
        for(int i=0; i<V[u].size(); i++)
        {
            int a=Size[V[u][i]];
            if(V[u][i]==pre)a=N-Size[u];
            for(int j=i+1; j<V[u].size(); j++)
            {
                int b=Size[V[u][j]];
                if(V[u][j]==pre)b=N-Size[u];
                cnt[a+b]++;
                if(cnt[a+b]>=mod)cnt[a+b]-=mod;
                //val[u]=val[u]-C(a+b,k);
                //if(val[u]<0)val[u]+=mod;
            }
            cnt[a]-=(V[u].size()-2);
            if(cnt[a]<0)cnt[a]+=mod;
            //val[u]=(val[u]+1LL*(V[u].size()-2)*C(a,k)%mod);
            //if(val[u]>mod)val[u]-=mod;
        }
    }
}


int main()
{
    init();
    scanf("%d",&N);
    int u,v;
    for(int i=1; i<N; i++)
    {
        scanf("%d %d",&u,&v);
        V[u].push_back(v);
        V[v].push_back(u);
    }
    dfs(1,0);
    for(k=1; k<=N; k++)
    {
        long long ans=sigle*C(N-1,k-1)%mod;
        ans=ans+(N-sigle)*C(N,k)%mod;
        for(int i=1;i<=N;i++)if(cnt[i])
        {
            ans=(ans-cnt[i]*C(i,k)%mod);
            if(ans<mod)ans+=mod;
        }
        long long tmp=mod_reverse(C(N,k),mod);
        ans=ans*tmp%mod;
        printf("%lld\n",ans);
    }
    return 0;
}

F Sum Of Digit
首先我们打个表可以发现sod(v)会返回(v-1)%15+1.
而我们只需要关注每段区间0到15出现的次数，这里就用线段树维护即可。
合并操作看代码，原理自己试几个数也就知道了。

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+5;
const int mod=1e9+7;
char s[maxn];
struct node
{
    int l,r;
    long long cnt[20];
    node(){memset(cnt,0,sizeof(cnt));}

}nodes[maxn<<2];
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1

int getnum(char s)
{
    if(s>='0'&&s<='9')return s-'0';
    return s-'A'+10;
}

void pushup(int rt)
{
    for(int i=0;i<16;i++)
    {
        nodes[rt].cnt[i]=(nodes[rt<<1].cnt[i]+nodes[rt<<1|1].cnt[i]);
        if(nodes[rt].cnt[i]>=mod)nodes[rt].cnt[i]-=mod;
    }
    for(int i=0;i<16;i++)
        for(int j=0;j<16;j++)
    {
        int num=i+j-1;
        if(num>=15)num-=15;
        num++;
        nodes[rt].cnt[num]=(nodes[rt].cnt[num]+nodes[rt<<1].cnt[i]*nodes[rt<<1|1].cnt[j]%mod);
        if(nodes[rt].cnt[num]>=mod)nodes[rt].cnt[num]-=mod;
    }
}

void build(int l,int r,int rt)
{
    if(l==r)
    {
        nodes[rt].l=l,nodes[rt].r=r;
        nodes[rt].cnt[getnum(s[l])]++;
        return;
    }
    nodes[rt].l=l,nodes[rt].r=r;
    int m=l+r>>1;
    build(ls);
    build(rs);
    pushup(rt);
}
int tab[30];
void init()
{
    tab[0]=1;
    tab[1]=1021;
    for(int i=2;i<=15;i++)
        tab[i]=1021LL*tab[i-1]%mod;
}

void update(int L,char c,int l,int r,int rt)
{
    if(l==r)
    {
        memset(nodes[rt].cnt,0,sizeof(nodes[rt].cnt));
        nodes[rt].cnt[getnum(c)]=1;
        return;
    }
    int m=l+r>>1;
    if(L<=m)update(L,c,ls);
    else update(L,c,rs);
    pushup(rt);
}

void Union(node &a,node b)
{
    node c=a;
    for(int i=0;i<=15;i++)
    {
        a.cnt[i]=(a.cnt[i]+b.cnt[i]);
        if(a.cnt[i]>=mod)a.cnt[i]-=mod;
    }
    for(int i=0;i<=15;i++)
        for(int j=0;j<=15;j++)
    {
        int num=(i+j-1);
        if(num>=15)num-=15;
        num++;
        a.cnt[num]=(a.cnt[num]+c.cnt[i]*b.cnt[j]%mod);
        if(a.cnt[num]>=mod)a.cnt[num]-=mod;
    }
}

node query(int L,int R,int l,int r,int rt)
{
    if(L<=l&&r<=R)
        return nodes[rt];
    int m=l+r>>1;
    node a;
    if(L<=m)
        Union(a,query(L,R,ls));
    if(R>m)
        Union(a,query(L,R,rs));
    return a;

}

int main()
{
    //freopen("in.txt","r",stdin);
    init();
    int n,q;
    scanf("%d %d",&n,&q);
    scanf("%s",s+1);
    build(1,n,1);
    int opt;
    while(q--)
    {
        scanf("%d",&opt);
        if(opt==1)
        {
            int p;
            char ch[3];
            scanf("%d %s",&p,ch);
            update(p,ch[0],1,n,1);
        }
        else
        {
            int l,r;
            scanf("%d %d",&l,&r);
            node a=query(l,r,1,n,1);
            long long sum=0;
            for(int i=0;i<=15;i++)
            {
                sum=(sum+a.cnt[i]*tab[i]%mod);
                if(sum>=mod)sum-=mod;
            }
            printf("%lld\n",sum);
        }
    }
    return 0;
}

G Coloring Tree
这题我也不知道如何解释，真的很难想。。
具体解释:https://www.nowcoder.com/discuss/88556?type=101&order=0&pos=1&page=1
放个证明给大家看。
这里写图片描述

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
vector<int>V[5005];
int d[5005][5005];
void dfs(int s,int now,int pre)
{
    for(int i=0;i<V[now].size();i++)
    {
        int v=V[now][i];
        if(v==pre)continue;
        d[s][v]=d[s][now]+1;
        dfs(s,v,now);
    }
}
int times[5005];
bool vis[5005];
int N,K,D;
long long solve(int dnum)
{
    for(int i=1;i<=N;i++)times[i]=K;
    queue<int>Q;
    memset(vis,0,sizeof(vis));
    Q.push(1);
    while(!Q.empty())
    {
        int u=Q.front();Q.pop();
        if(vis[u])continue;
        vis[u]=1;
        for(int i=1;i<=N;i++)if(vis[i]!=1&&d[u][i]<dnum)
            times[i]--;
        for(int v:V[u])
            Q.push(v);
    }
    long long res=1;
    for(int i=1;i<=N;i++)
        res=res*times[i]%mod;
    return res;
}

int main()
{

    scanf("%d %d %d",&N,&K,&D);
    int u,v;
    for(int i=1;i<N;i++)
    {
        scanf("%d %d",&u,&v);
        V[u].push_back(v);
        V[v].push_back(u);
    }
    for(int i=1;i<=N;i++)
        dfs(i,i,0);
    printf("%lld\n",((solve(D)-solve(D+1))%mod+mod)%mod);


}

H Diff-prime Pairs
卡了2个半小时才把这个题写出来。。
对于每个质数prm[i]而言，每隔N/prm[i],就会有一个数拥有这个质数，那么它提供的贡献就是i*N/prm[i] (该贡献意为它和比他小的质数产生的，每个因子包含这个质数的数可以将它的幂减1，然后将小于这个质数的质数的幂加1，产生一对符合题意的数)。
那么最后的答案再乘以2即可。

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e7+5;
bool is[maxn];
int prm[maxn], id;
int Maxprm[maxn];
void init()
{
    memset(is, 1, sizeof(is));
    is[0] = is[1] = 0;
    int k = 0;
    prm[k++] = 2;
    for (int i = 4;i < maxn;i+=2)is[i] = 0,Maxprm[i]=max(Maxprm[i],0);
    int e = (int)sqrt(0.0 + maxn) + 1;
    int i;
    for (i = 3;i < e;i += 2)if(is[i])
    {
        prm[k++] = i;
        Maxprm[i]=max(Maxprm[i],k-1);
        for (int s = i, j = i+i;j < maxn;j += s)
            is[j] = 0,Maxprm[j]=max(Maxprm[j],k-1);
    }
    for (;i < maxn;i += 2)if (is[i])
    {
        prm[k++] = i;
        //Maxprm[i]=max(Maxprm[i],k-1);
        for(int s=i,j=i;j<maxn;j+=s)
            Maxprm[j]=max(Maxprm[j],k-1);
    }
    id = k;
}
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}


int main()
{
    init();
    int N;
    scanf("%d",&N);
    long long ans=0;
    for(int i=1;i<id;i++)
    {
        if(prm[i]>N)break;
        ans+=N/prm[i]*i;
        //cout<<i<<' '<<Maxprm[i]<<endl;
        //if(Maxprm[i]==0)cout<<i<<endl;
    }
    cout<<ans*2LL<<endl;
}

I Expected Size of Random Convex Hull
这题没太多的解释，打表找规律吧。。
虽然它的规律是调和级数。。但博主不知道它的具体数学意义是啥。。

#include<bits/stdc++.h>
using namespace std;
int main()
{
    int x,y;
    for(int i=1;i<=3;i++)
        scanf("%d %d",&x,&y);
    int N;
    scanf("%d",&N);
    double ans=0;
    for(int i=1;i<N;i++)
        ans+=(1.0/i)*2;
    printf("%.12f\n",ans);


}

J Distance to Work
主要是多边形与圆交面积，而且这题卡精度卡的莫名其妙。。

#include<bits/stdc++.h>
using namespace std;
const int maxn=205;
const double eps=1e-10;
const double PI=acos(-1.0);
int sgn(double x)
{
    if(fabs(x)<eps)return 0;
    if(x>0)return 1;
    return -1;
}
struct Point
{
    double x,y;
    Point(){}
    Point(double sx,double sy):x(sx),y(sy){}
    void read()
    {
        scanf("%lf %lf",&x,&y);
    }
    Point operator-(const Point &b)const
    {
        return Point(x-b.x,y-b.y);
    }
    double operator^(const Point &b)const
    {
        return x*b.y-y*b.x;
    }
    double operator*(const Point &b)const
    {
        return x*b.x+y*b.y;
    }
    double getdis()
    {
        return sqrt(x*x+y*y);
    }
    Point operator+(const Point &b)const
    {
        return Point(x+b.x,y+b.y);
    }

}Poly[maxn];
double PS;
int n;

bool InCircle(Point a,double r)
{
    if(sgn(a.getdis()-r)<=0)return 1;
    return 0;
}

typedef const Point CP;

Point intersection(CP &u1,CP &u2,CP &v1,CP &v2)
{
    Point ret=u1;
    double t=((u1.x-v1.x)*(v1.y-v2.y)-(u1.y-v1.y)*(v1.x-v2.x))
            /((u1.x-u2.x)*(v1.y-v2.y)-(u1.y-u2.y)*(v1.x-v2.x));
    ret.x+=(u2.x-u1.x)*t;
    ret.y+=(u2.y-u1.y)*t;
    return ret;
}

double dis(CP &a,CP &b)
{
    double x=a.x-b.x;
    double y=a.y-b.y;
    return sqrt(x*x+y*y);
}

void intersection_line_circle(CP &c,double r,CP &l1,CP &l2,Point &p1,Point &p2)
{
    Point p=c;
    double t;
    p.x+=l1.y-l2.y;
    p.y+=l2.x-l1.x;
    p=intersection(p,c,l1,l2);
    t=sqrt(r*r-dis(p,c)*dis(p,c))/dis(l1,l2);
    p1.x=p.x+(l2.x-l1.x)*t;
    p1.y=p.y+(l2.y-l1.y)*t;
    p2.x=p.x-(l2.x-l1.x)*t;
    p2.y=p.y-(l2.y-l1.y)*t;
}

bool Inmid(Point a,Point b,Point c)
{
    if(sgn((c-a).getdis()+(c-b).getdis()-(a-b).getdis())<=0)return 1;
    return 0;
}

double getang(Point a,Point b)
{
    double ang=atan2(b.y,b.x)-atan2(a.y,a.x);
    if(sgn(ang+PI)<0)ang+=2*PI;
    if(sgn(ang-PI)>0)ang-=2*PI;
    return ang;


}

double xmult(Point a,Point b,Point c)
{
    return (a-c)^(b-c);
}

double disptoline(CP &p,CP &l1,CP &l2)
{
    return fabs(xmult(p,l1,l2))/dis(l1,l2);
}
int intersect_seg_circle(CP &c,double r,CP &l1,CP &l2)
{
    double t1=dis(c,l1)-r,t2=dis(c,l2)-r;
    Point t=c;
    if(t1<eps||t2<eps)
        return t1>-eps||t2>-eps;
    t.x+=l1.y-l2.y;
    t.y+=l2.x-l1.x;
    return xmult(l1,c,t)*xmult(l2,c,t)<eps&&disptoline(c,l1,l2)-r<eps;
}
double sqr(double x){return x*x;}
double dis2(Point a,Point b)
{
    double x=a.x-b.x;
    double y=a.y-b.y;
    return x*x+y*y;
}

double area_triangle(CP &p1,CP &p2,CP &p3)
{
    return fabs(xmult(p1,p2,p3))/2;
}

double area_circle_angle(CP &p1,CP &p2,CP &c,double r)
{
    double alpha=fabs(atan2(p1.y-c.y,p1.x-c.x)-atan2(p2.y-c.y,p2.x-c.x));
    if(alpha>PI)alpha=2*PI-alpha;
    return alpha/2*r*r;
}
double dmult(CP &p1,CP &p2,CP &p0)
{
    return (p1-p0)*(p2-p0);
}

double disptoseg(CP &p,CP &l1,CP &l2)
{
    Point t=p;
    t.x+=l1.y-l2.y,t.y+=l2.x-l1.x;
    if(xmult(l1,t,p)*xmult(l2,t,p)>eps)
        return dis(p,l1)<dis(p,l2)?dis(p,l1):dis(p,l2);
    return fabs(xmult(p,l1,l2))/dis(l1,l2);
}

double area_triangle_circle(CP &c,double r,CP &p1,CP &p2)
{
    double x=xmult(p2,c,p1);
    int flag=sgn(x);
    if(flag==0)return 0;
    double r2=sqr(r);
    double s=0,l1=dis2(p1,c),l2=dis2(p2,c);
    if((l1<=r2)&&(l2<=r2))
        return area_triangle(p2,c,p1)*flag;
    if((l1>r2)&&(l2>r2))
    {
        Point p3,p4;
        s=area_circle_angle(p2,p1,c,r);
        if(disptoseg(c,p1,p2)<r)
        {
            intersection_line_circle(c,r,p1,p2,p3,p4);
            if(dis2(p3,p1)>dis2(p4,p1))
                swap(p3,p4);
            s-=area_circle_angle(p3,p4,c,r)-area_triangle(p3,c,p4);
        }
        return s*flag;
    }
    if(l1<l2)
    {
        Point p3,p4;
        intersection_line_circle(c,r,p1,p2,p3,p4);
        if(dmult(p3,p2,p1)<=0)p3=p4;
        s=area_triangle(p1,p3,c)+area_circle_angle(p3,p2,c,r);
        return s*flag;
    }
    else
    {
        Point p3,p4;
        intersection_line_circle(c,r,p1,p2,p3,p4);
        if(dmult(p3,p1,p2)<=0)p3=p4;
        s=area_triangle(p2,p3,c)+area_circle_angle(p3,p1,c,r);
        return s*flag;
    }
}


double getS(double r)
{
   double ans=0;
   for(int i=1;i<=n;i++)
   {
       ans+=area_triangle_circle(Point(0,0),r,Poly[i],Poly[i%n+1]);
   }
   return fabs(ans);
}

bool check(double r,double haveS)
{
    double thes=getS(r);
    if(sgn(thes-haveS)<=0)
        return 1;
    return 0;
}

void solve()
{
    double x,y,P,Q;
    scanf("%lf %lf %lf %lf",&x,&y,&P,&Q);
    for(int i=1;i<=n;i++)Poly[i]=Poly[i]-Point(x,y);
    P=Q-P;
    //double haveS=PS*p;
    double l=0,r=1e7;
    for(int t=1;t<=100;t++)
    {
        double mid=(l+r)/2;
        double a=getS(mid);
        if(a*Q<P*PS)
            l=mid;
        else r=mid;
    }
    for(int i=1;i<=n;i++)Poly[i]=Poly[i]+Point(x,y);
    printf("%.12f\n",l+eps);

}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        Poly[i].read();
    for(int i=2;i<n;i++)
        PS+=((Poly[i]-Poly[1])^(Poly[i+1]-Poly[1]));
    PS*=0.5;
    if(PS<0)
    {
        reverse(Poly+1,Poly+1+n);
        PS*=-1;
    }
    int Q;
    scanf("%d",&Q);
    for(int i=1;i<=Q;i++)
        solve();
    return 0;
}

比赛的时候迷迷糊糊的。。会的题都要写好久。。看来比赛前一晚上不能熬夜了。

大学要有梦想
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值