2018.8.9提高B组模拟考试

Posted water-radish

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了2018.8.9提高B组模拟考试相关的知识,希望对你有一定的参考价值。

今天连续3道题都出锅...无F♂A可说

T1 题意简述:jzoj5775

   解题思路:还算简单...

             考虑先暴力算出(1,1)的魔音值,然后递推。

             发现 (x,y)->(x+1,y) 魔音值增加了1~x行中起始点的数量,减少了x+1~n行中起始点的数量。

                  (x,y)->(x-1,y) 魔音值增加了x~n行中起始点的数量,减少了1~x-1行中起始点的数量。

                  (x,y)->(x,y+1) 魔音值增加了1~y列中起始点的数量,减少了y+1~n列中起始点的数量。

                  (x,y)->(x,y-1) 魔音值增加了y~n列中起始点的数量,减少了1~y-1列中起始点的数量。

             发现可以求每行和每列的前缀和,可以把时间复杂度压到O(2n)。

             有细心的观众可能会说:z呢?

             emmm...由于出题人乱出数据,导致所有数据点的z都默认为1。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
ll n,m,k,sumx[100001],sumy[100001],sum;
int main()
{
    freopen("shuru.in","r",stdin);
    freopen("shuru.out","w",stdout);
    scanf("%lld%lld%lld",&n,&m,&k);
    for(ll i=1;i<=m;i++)
    {
        ll u,v,w;
        scanf("%lld%lld%lld",&u,&v,&w);
        sumx[u]++;
        sumy[v]++;
        sum+=w;
    }
    for(ll i=1;i<=n;i++)
        sum+=sumx[i]*(i-1)+sumy[i]*(i-1);
    for(ll i=2;i<=n;i++)
        sumx[i]+=sumx[i-1],sumy[i]+=sumy[i-1];
    ll xsum=sum,ysum=sum,mnx=sum,ansx=1,mny=sum,ansy=1;
    for(ll i=2;i<=n;i++)
    {
        xsum+=2*sumx[i-1]-sumx[n];
        if(xsum<mnx) mnx=xsum,ansx=i;
        ysum+=2*sumy[i-1]-sumy[n];
        if(ysum<mny) mny=ysum,ansy=i;
    }
    printf("%lld
%lld %lld
",mnx+mny-sum,ansx,ansy);
    return 0;
}

 


 

T2 题意简述:jzoj5776

   解题思路:和第一题基本相同。

             先dfs求出一个点的结果,同时把每个点的子树大小算出。

             然后dfs,每次查询到一个点时修改它与它父亲的子树大小及边权,然后算出这个点的结果。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<queue>
#define INF 0x3f3f3f3f
#define ll long long
using namespace std;
ll n,mn=INF,pnt,cnt,head[700001],rdc[700001],siz[700001],ans[700001];
struct uio{
    ll nxt,to,val;
}edge[1400001];
void add(ll x,ll y,ll z)
{
    edge[++cnt].nxt=head[x];
    edge[cnt].to=y;
    edge[cnt].val=z;
    head[x]=cnt;
}
ll dfs(ll x,ll fa)
{
    siz[x]=1;
    ll tmp=0;
    for(ll i=head[x];i;i=edge[i].nxt)
    {
        ll y=edge[i].to;
        if(y==fa) continue;
        tmp+=dfs(y,x)+edge[i].val*siz[y];
        siz[x]+=siz[y];
    }
    return tmp;
}
void dt(ll x,ll fa)
{
    for(ll i=head[x];i;i=edge[i].nxt)
    {
        ll y=edge[i].to;
        if(y==fa) continue;
        ans[y]=ans[x]-edge[i].val*siz[y];
        siz[x]-=siz[y],siz[y]+=siz[x];
        edge[i].val+=rdc[x]-rdc[y];
        ans[y]+=edge[i].val*siz[x];
        if(ans[y]<mn) mn=ans[y],pnt=y;
        if(ans[y]==mn&&y<pnt) pnt=y;
        dt(y,x);
        siz[y]-=siz[x],siz[x]+=siz[y];
        edge[i].val+=rdc[y]-rdc[x];
    }
}
int main()
{
    freopen("yggdrasil.in","r",stdin);
    freopen("yggdrasil.out","w",stdout);
    scanf("%lld",&n);
    for(ll i=1;i<=n;i++)
        scanf("%lld",&rdc[i]);
    for(ll i=1;i<n;i++)
    {
        ll u,v,w;
        scanf("%lld%lld%lld",&u,&v,&w);
        add(u,v,w-rdc[u]);
        add(v,u,w-rdc[v]);
    }
    ans[1]=dfs(1,0),mn=ans[1],pnt=1;
    dt(1,0);
    printf("%lld
%lld
",pnt,mn);
    return 0;
}

 


 

T3 题意简述:jzoj5786

   解题思路:很...麻...烦...

             倍增LCA+线段树维护 <-- 这是出题人的思路。

             由Menteur_Hxy大佬亲身试验,倍增是会T的

             由ErkkiErkko大佬亲身试验,set也是会T的

             由本蒟蒻亲身试验,就算用树剖+线段树,也是会RE的。

             ......

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
#include<set>
#define INF 0x3f3f3f3f
using namespace std;
int n,m,cnt,tot,head[800001];
int hvyson[800001],fa[800001],siz[800001],dep[800001],top[800001],dfsno[800001],vis[800001],mp[800001];
int num[3200001];
struct uio{
    int nxt,to;
}edge[1600001];
//set<int> st;
//set<int>::iterator it;
void add(int x,int y)
{
    edge[++cnt].nxt=head[x];
    edge[cnt].to=y;
    head[x]=cnt;
}
void dfs1(int x,int f,int depth)
{
    dep[x]=depth;
    fa[x]=f;
    siz[x]=1;
    int maxson=-1;
    for(int i=head[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(y==f)
            continue;
        dfs1(y,x,depth+1);
        siz[x]+=siz[y];
        if(siz[y]>maxson)
        {
            hvyson[x]=y;
            maxson=siz[y];
        }
    }
}
void dfs2(int x,int topf)
{
    dfsno[x]=++tot;
    mp[tot]=x;
    top[x]=topf;
    if(!hvyson[x])
        return;
    dfs2(hvyson[x],topf);
    for(int i=head[x];i;i=edge[i].nxt)
    {
        int y=edge[i].to;
        if(y==fa[x]||y==hvyson[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;
}
void revise(int now,int l,int r,int x,int y)
{
    if(l==r)
    {
        num[now]+=y;
        return;
    }
    int mid=(l+r)>>1;
    if(x<=mid) revise(now<<1,l,mid,x,y);
    else revise(now<<1|1,mid+1,r,x,y);
    num[now]=num[now<<1]+num[now<<1|1];
}
int query(int now,int l,int r,int x)
{
    if(l==r&&num[now]) return 1;
    else if(l==r&&!num[now]) return 0;
    int mid=(l+r)>>1;
    if(x<=mid) return query(now<<1,l,mid,x);
    else return query(now<<1|1,mid+1,r,x);
}
int query1(int now,int l,int r,int L,int R)
{
    if(!num[now]) return INF;
    if(l==r) return l;
    if(L<=l&&r<=R)
    {
        int mid=(l+r)>>1;
        if(num[now<<1|1]) return query1(now<<1|1,mid+1,r,L,R);
        else return query1(now<<1,l,mid,L,R);
    }
    int mid=(l+r)>>1;
    if(mid>=R) return query1(now<<1,l,mid,L,R);
    else
    {
        int tmp=INF;
        if(num[now<<1|1]) tmp=query1(now<<1|1,mid+1,r,L,R);
        if(tmp!=INF) return tmp;
        else return query1(now<<1,l,mid,L,R);
    }
}
int query2(int now,int l,int r,int L,int R)
{
    if(!num[now]) return INF;
    if(l==r) return l;
    if(L<=l&&r<=R)
    {
        int mid=(l+r)>>1;
        if(num[now<<1]) return query2(now<<1,l,mid,L,R);
        else return query2(now<<1|1,mid+1,r,L,R);
    }
    int mid=(l+r)>>1;
    if(mid+1<=L) return query2(now<<1|1,mid+1,r,L,R);
    else
    {
        int tmp=INF;
        if(num[now<<1]) tmp=query2(now<<1,l,mid,L,R);
        if(tmp!=INF) return tmp;
        else return query2(now<<1|1,mid+1,r,L,R);
    }
}
int main()
{
    freopen("watch.in","r",stdin);
    freopen("watch.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
    {
        int u;
        scanf("%d",&u);
        add(i+1,u),add(u,i+1);
    }
    dfs1(1,0,1),dfs2(1,1);
    for(int i=1;i<=m;i++)
    {
        int u;
        scanf("%d",&u);
        if(u>0&&!vis[u]) {revise(1,1,n,dfsno[u],1);vis[u]=1;continue;}
        else if(u>0&&vis[u]) {revise(1,1,n,dfsno[u],-1);vis[u]=0;continue;}
        if(!num[1]) {printf("0
");continue;}
        u=-u;
        int c=query(1,1,n,dfsno[u]);
        if(c) {printf("%d
",u);continue;}
        int a=query1(1,1,n,1,dfsno[u]-1);
        int b=query2(1,1,n,dfsno[u]+1,n);
        if(a!=INF) a=lca(u,mp[a]);
        if(b!=INF) b=lca(u,mp[b]);
        if(b==INF) printf("%d
",a);
        else if(a==INF) printf("%d
",b);
        else printf("%d
",(dep[a]>dep[b]? a:b));
//        if(u>0&&!vis[u]) {st.insert(dfsno[u]);vis[u]=1;continue;}
//        else if(u>0&&vis[u])  {st.erase(dfsno[u]);vis[u]=0;continue;}
//        if(st.empty()) {printf("0
");continue;}
//        u=-u;
//        if(st.count(u)) {printf("%d
",u);continue;}
//        it=st.upper_bound(dfsno[u]);
//        if(it!=st.end())
//        {
//            a=lca(u,mp[*it]);
//            if(it==st.begin())
//                {printf("%d
",a);continue;}
//            it--;
//            b=lca(u,mp[*it]);
//            printf("%d
",(dep[a]>dep[b]? a:b));
//            continue;
//        }
//        it--;
//        b=lca(u,mp[*it]);
//        printf("%d
",b);
    }
    return 0;
}

 

以上是关于2018.8.9提高B组模拟考试的主要内容,如果未能解决你的问题,请参考以下文章

2018.8.10提高B组模拟考试

2018.8.6提高A组模拟考试

2021.7.21提高B组模拟8T1 好数(模拟)

2021.8.13提高B组模拟5T1 Brothers(暴力)

2021.8.10提高B组模拟2T1 单峰(快速幂)

2021.8.10提高B组模拟2T1 单峰(快速幂)