两侧同时换到之前的修订记录 前一修订版 后一修订版 | 前一修订版 | ||
2020-2021:teams:legal_string:组队训练比赛记录:contest15 [2021/08/22 19:08] jxm2001 |
2020-2021:teams:legal_string:组队训练比赛记录:contest15 [2021/08/27 21:42] (当前版本) jxm2001 |
||
---|---|---|---|
行 4: | 行 4: | ||
^ 题目 ^ 蒋贤蒙 ^ 王赵安 ^ 王智彪 ^ | ^ 题目 ^ 蒋贤蒙 ^ 王赵安 ^ 王智彪 ^ | ||
- | | A | 2 | 0 | 0 | | + | | A | 2 | 1 | 0 | |
- | | B | 2 | 0 | 0 | | + | | B | 2 | 1 | 0 | |
| C | 2 | 1 | 0 | | | C | 2 | 1 | 0 | | ||
| D | 0 | 0 | 0 | | | D | 0 | 0 | 0 | | ||
| F | 0 | 0 | 0 | | | F | 0 | 0 | 0 | | ||
- | | G | 2 | 0 | 0 | | + | | G | 2 | 2 | 0 | |
| I | 1 | 0 | 2 | | | I | 1 | 0 | 2 | | ||
| J | 2 | 0 | 2 | | | J | 2 | 0 | 2 | | ||
行 161: | 行 161: | ||
s=(s+1LL*ans.f[p][i]*A[i])%mod; | s=(s+1LL*ans.f[p][i]*A[i])%mod; | ||
enter(s); | enter(s); | ||
+ | return 0; | ||
+ | } | ||
+ | </code> | ||
+ | </hidden> | ||
+ | |||
+ | ===== B. Best Subgraph ===== | ||
+ | |||
+ | ==== 题意 ==== | ||
+ | |||
+ | 定义 $k-\text{degree}$ 子图为每个点在子图中度数至少为 $k$ 的连通极大子图。 | ||
+ | |||
+ | 定义每个 $G$ 的子图 $S$ 的分数为 $M\times m(S)-N\times n(S)+B\times b(S)$。 | ||
+ | |||
+ | 其中 $m(S)$ 表示 $S$ 的边数,$n(S)$ 表示 $S$ 的点数,$b(S)$ 表示 $|\{(u,v)|(u,v)\in G,u\in S,v\not\in S\}|$。 | ||
+ | |||
+ | 求分数最大的 $k-\text{degree}$ 子图,并求出分数最大的 $k-\text{degree}$ 子图中 $k$ 的最大值。 | ||
+ | |||
+ | 输入保证图连通。 | ||
+ | |||
+ | ==== 题解 ==== | ||
+ | |||
+ | 首先,定义 $V(k)$ 表示所有 $k-\text{degree}$ 子图的并集,易知 $V(k+1)\subseteq V(k)$。 | ||
+ | |||
+ | 构建分层图,其中第 $k$ 层的点集为 $V(k)-V(k+1)$,同时对每个点 $u\in V(k)-V(k+1)$,$w(u)=k$。 | ||
+ | |||
+ | 然后考虑按层加点,动态维护当前每个 $k-\text{degree}$ 子图的答案。对于每个新加入的点 $u$,首先他本身产生贡献 $-N$。 | ||
+ | |||
+ | 对于 $u$ 的所有连边 $(u,v)$,如果 $w(v)\gt w(u)$,则 $(u,v)$ 会被加入 $m(S)$,同时从 $b(S)$ 删除,产生贡献 $M-B$。 | ||
+ | |||
+ | 如果 $w(v)=w(u)$,则 $(u,v)$ 也会被加入 $m(S)$,但为了贡献被计算两次,于是每个点的贡献为 $\frac M2$。 | ||
+ | |||
+ | 如果 $w(v)\lt w(u)$,则 $(u,v)$ 会被加入 $b(S)$,产生贡献 $B$。 | ||
+ | |||
+ | 于是上述过程可以将所有贡献都转化为每个点的点权。考虑加入点时并查集维护每个连通块的点权和。 | ||
+ | |||
+ | 每层点全部加完后查询每个连通块的点权和的最大值即为所有 $k-\text{degree}$ 子图的最大答案。 | ||
+ | |||
+ | 至于如果计算每个点的 $w(u)$。首先输入保证图连通,所有图 $G$ 本身就是 $1-\text{degree}$ 子图。 | ||
+ | |||
+ | 考虑先删去所有度数为 $1$ 的点,删点后可能会导致原来度数不为 $1$ 的点度数不大于 $1$,继续删除,直到无点可删,得到 $2-\text{degree}$ 子图。 | ||
+ | |||
+ | 将删去的点的 $w(u)$ 全部设为 $1$,然后再求 $3-\text{degree}$ 子图不断重复上述过程,即可得到所有 $w(u)$。 | ||
+ | |||
+ | 时间复杂度 $O\left(\text{Na}(n+m)\right)$。 | ||
+ | |||
+ | <hidden 查看代码> | ||
+ | <code cpp> | ||
+ | const int MAXN=1e6+5; | ||
+ | const LL inf=1e18; | ||
+ | struct Edge{ | ||
+ | int to,next; | ||
+ | }edge[MAXN<<1]; | ||
+ | int head[MAXN],edge_cnt; | ||
+ | void Insert(int u,int v){ | ||
+ | edge[++edge_cnt]=Edge{v,head[u]}; | ||
+ | head[u]=edge_cnt; | ||
+ | } | ||
+ | bool vis[MAXN]; | ||
+ | vector<int> q[MAXN],vec[MAXN]; | ||
+ | int deg[MAXN],w[MAXN],p[MAXN]; | ||
+ | LL sum[MAXN]; | ||
+ | int Find(int x){ | ||
+ | return x==p[x]?x:p[x]=Find(p[x]); | ||
+ | } | ||
+ | int main() | ||
+ | { | ||
+ | int n=read_int(),m=read_int(),M=read_int(),N=read_int(),B=read_int(); | ||
+ | while(m--){ | ||
+ | int u=read_int(),v=read_int(); | ||
+ | Insert(u,v); | ||
+ | Insert(v,u); | ||
+ | deg[u]++; | ||
+ | deg[v]++; | ||
+ | } | ||
+ | _rep(i,1,n) | ||
+ | q[deg[i]].push_back(i); | ||
+ | _for(k,1,MAXN){ | ||
+ | _for(j,0,q[k].size()){ | ||
+ | int u=q[k][j]; | ||
+ | if(vis[u])continue; | ||
+ | vis[u]=true; | ||
+ | w[u]=k; | ||
+ | vec[k].push_back(u); | ||
+ | for(int i=head[u];i;i=edge[i].next){ | ||
+ | int v=edge[i].to; | ||
+ | deg[v]--; | ||
+ | q[deg[v]].push_back(v); | ||
+ | } | ||
+ | } | ||
+ | } | ||
+ | _rep(u,1,n){ | ||
+ | p[u]=u; | ||
+ | sum[u]=-N; | ||
+ | for(int i=head[u];i;i=edge[i].next){ | ||
+ | int v=edge[i].to; | ||
+ | if(w[u]<w[v]){ | ||
+ | sum[u]+=M; | ||
+ | sum[u]-=B; | ||
+ | } | ||
+ | else if(w[u]==w[v]&&u<v) | ||
+ | sum[u]+=M; | ||
+ | else if(w[u]>w[v]) | ||
+ | sum[u]+=B; | ||
+ | } | ||
+ | } | ||
+ | int st=MAXN-1; | ||
+ | while(vec[st].empty())st--; | ||
+ | pair<LL,int> ans=make_pair(-inf,0); | ||
+ | for(int k=st;k;k--){ | ||
+ | for(int u:vec[k]){ | ||
+ | for(int i=head[u];i;i=edge[i].next){ | ||
+ | int v=edge[i].to; | ||
+ | if(w[v]<k)continue; | ||
+ | int x=Find(u),y=Find(v); | ||
+ | if(x!=y){ | ||
+ | p[x]=y; | ||
+ | sum[y]+=sum[x]; | ||
+ | } | ||
+ | } | ||
+ | } | ||
+ | for(int u:vec[k]) | ||
+ | ans=max(ans,make_pair(sum[Find(u)],k)); | ||
+ | } | ||
+ | space(ans.second);enter(ans.first); | ||
return 0; | return 0; | ||
} | } |