Public Round #1 题解

A

题意

分析

首先发现每次操作只和左右两个人的值有关，所以考虑差分
差分后的数组为b，那么每次的操作相当于合并两个相等的$b_i,b_{i+1}$为$2b_i$
发现可以通过分组来计算，正数和负数肯定不能合并到一起，然后就是$\frac{b_i}{lowbit(b_i)}$相同的才可以一起合并
每段我们通过dp计算最小的长度，设$f_i$表示到 i 位的长度最小为多少
设$g_{i,j}$表示第 i 位为右端点，合并到$2^j$的左端点

总复杂度为$O(nlogv)$

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=3e5+5;
const int inf=0x3f3f3f3f;
int a[maxn],b[maxn],op[maxn];
int T,n,ans;
int lowbit(int x)
{
    return x&(-x);
}
int calc(int x)
{
    if(!x) return x;
    return x/lowbit(x);
}
int c[maxn],g[maxn][30],f[maxn];
void solve(int l,int r)
{
    if(!b[l]) ans++;
    else
    {
        for(int i=l;i<=r;i++) c[i]=int(log((double)(b[i]&-b[i]))/log(2.0)+0.5);
        for(int i=l;i<=r;i++) g[i][c[i]]=i;
        for(int j=1;j<=30;j++)
            for(int i=l;i<=r;i++)
                if(c[i]<j && g[i][j-1]!=inf && g[i][j-1]>l) g[i][j]=g[g[i][j-1]-1][j-1];
        f[l-1]=0;
        for(int i=l;i<=r;i++)
        {
            f[i]=i-l+1;
            for(int j=0;j<=30;j++)
                if(g[i][j]!=inf) f[i]=min(f[i],f[g[i][j]-1]+1);
        }
        ans+=f[r];
    }
}
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n); ans=0;
        for(int i=0;i<=n;i++)
        {
            f[i]=inf;
            for(int j=0;j<=30;j++) g[i][j]=inf;
        }
        for(int i=1;i<=n;i++) scanf("%d",&a[i]),b[i-1]=a[i]-a[i-1];
        op[1]=b[1]<0?-1:1;
        b[1]=abs(b[1]);
        for(int i=2,j=1;i<=n;i++)
        {
            op[i]=b[i]<0?-1:1;
            b[i]=abs(b[i]);
            if(i==n || op[i]!=op[i-1] || calc(b[i])!=calc(b[i-1]))
                solve(j,i-1),j=i;
        }
        printf("%d\n",ans+1);
    }
    return 0;
}

B

题面

分析

考虑我们对所有连通块求出最小生成树肯定会更优
那么我们就要对于这个森林求出可以删除的边权的最大和

可以使用最大费用最大流来解决
S->i (1,0) 表示从起点到每个守卫的边，每个守卫向其可以被派遣去的点连容量为 1 费用为 0 的边。
每个点拆成入点和出点，所有连向它的边连到入点，所有连出去的边连到出点，入点向出点连容量为 1 费用为 0 的边。
每个点的出点向汇点连容量为 1 费用为“它到父亲的边权”的边。
特别地，根向汇点连容量为 1 费用为 ∞ 的边。

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int maxn=30005;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
int head[maxn],tot=1;
int S,T,n,m,k;
struct edge
{
    int to,nxt,v;
    ll c;
}e[maxn<<3];
void add(int x,int y,int z,ll c)
{
    e[++tot].to=y; e[tot].nxt=head[x]; e[tot].v=z; e[tot].c=c;  head[x]=tot;
    e[++tot].to=x; e[tot].nxt=head[y]; e[tot].v=0; e[tot].c=-c; head[y]=tot;
}
struct edges
{
    int x,y,w;
}E[300*305];
bool cmp(edges x,edges y)
{
    return x.w<y.w;
}
int f[maxn];
int find(int x)
{
    if(x==f[x]) return x;
    return f[x]=find(f[x]);
}
vector <PII> G[maxn];
int block;
void dfs(int u,int fa=0)
{
    add(u,u+n,1,0);
    if(!fa) add(u+n,T,1,1e9);
    else add(u+n,fa,1,0);
    for(auto to:G[u])
    {
        int v=to.first;
        if(v==fa) continue;
        add(v+n,T,1,to.second);
        dfs(v,u);
    }   
}
int maxflow;
ll cost,dis[maxn];
int pre[maxn],vis[maxn];
bool spfa()
{
    for(int i=S;i<=T;i++) vis[i]=0,pre[i]=-1,dis[i]=-INF;
    dis[S]=0; queue<int> q;
    q.push(S); vis[S]=1;
    while(!q.empty())
    {
        int u=q.front(); q.pop();
        vis[u]=0;
        for(int i=head[u];i;i=e[i].nxt)
        {
            int to=e[i].to;
            if(e[i].v && dis[to]<dis[u]+e[i].c)
            {
                dis[to]=dis[u]+e[i].c;
                pre[to]=i;
                if(!vis[to])
                {
                    vis[to]=1;
                    q.push(to);
                }
            }
        }
    }
    return pre[T]!=-1;
}
void mcmf()
{
    while(spfa())
    {
        int flow=inf;       
        for(int i=T;i!=S;i=e[pre[i]^1].to) flow=min(flow,e[pre[i]].v);
        for(int i=T;i!=S;i=e[pre[i]^1].to) e[pre[i]].v-=flow,e[pre[i]^1].v+=flow;
        maxflow+=flow; cost+=1ll*flow*dis[T];
    }
}
int main()
{
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=m;i++) scanf("%d%d%d",&E[i].x,&E[i].y,&E[i].w);
    sort(E+1,E+m+1,cmp);
    for(int i=1;i<=n;i++) f[i]=i;
    ll sum=0;
    for(int i=1;i<=m;i++)
    {
        int fx=find(E[i].x);
        int fy=find(E[i].y);
        if(fx==fy) continue;
        f[fx]=fy; sum+=E[i].w;
        G[E[i].x].push_back(make_pair(E[i].y,E[i].w));
        G[E[i].y].push_back(make_pair(E[i].x,E[i].w));
    }
    S=0; T=2*n+k+1;
    for(int i=1;i<=n;i++) if(f[i]==i) block++,dfs(i);
    for(int i=1;i<=k;i++)
    {
        int gs,x; scanf("%d",&gs);
        add(S,n+n+i,1,0);
        for(int j=1;j<=gs;j++)
        {
            scanf("%d",&x);
            add(n+n+i,x,1,0);
        }
    }
    mcmf();
    if(maxflow!=k || cost<1ll*block*1e9)
    {
        printf("-1\n");
        return 0;
    }
    printf("%lld\n",block*1000000000ll+sum-cost);
    return 0;
}

C

题面

分析