这道题一看就是可持久化并查集 然后我就愉快的yy了一波 还是错掉了qwqwqwqwq
方法是对的 就是我每次在树上查询$fa$的时候我还压缩了路径 导致这玩意空间炸掉了
所以要保证时间复杂度 就启发式合并 也就是$size$小的往$size$大的搞
这样子就保证每次合并的时候连通块元素个数每次至少乘以$2$ 也就是保证了层数是$log$级的
代码
#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + 5;
int n,f[90 * N],ls[90 * N],rs[90 * N],size[90 * N],root[N];
int cnt,m,fa[N];
int build(int o, int l ,int r) {
int nd = ++ cnt;
if(l == r) {
f[nd] = l; size[nd] = 1;
return nd;
}
int mid = (l + r) >> 1;
ls[nd] = build(2 * o, l, mid);
rs[nd] = build(2 * o + 1, mid + 1, r);
return nd;
}
void Init( ) {
scanf("%d%d",& n,& m);
root[0] = build(1, 1, n);
}
int modify_fa(int pre, int o, int l, int r, int pos, int del) {
int nd = ++ cnt;
f[nd] = f[pre],ls[nd] = ls[pre],rs[nd] = rs[pre],size[nd] = size[pre];
if(l == r) {
f[nd] = del; return nd;
}
int mid = (l + r) >> 1;
if(pos <= mid) {
ls[nd] = modify_fa(ls[pre], 2 * o, l, mid, pos, del);
}
else rs[nd] = modify_fa(rs[pre], 2 * o + 1, mid + 1, r, pos, del);
return nd;
}
int modify_size(int pre, int o, int l, int r, int pos, int del) {
int nd = ++ cnt;
f[nd] = f[pre],ls[nd] = ls[pre],rs[nd] = rs[pre],size[nd] = size[pre] + del;
if(l == r) {
return nd;
}
int mid = (l + r) >> 1;
if(pos <= mid) {
ls[nd] = modify_size(ls[pre], 2 * o, l, mid, pos, del);
}
else rs[nd] = modify_size(rs[pre], 2 * o + 1, mid + 1, r, pos, del);
return nd;
}
int query_fa(int nd, int o, int l, int r, int pos) {
if(l == r) {
return f[nd];
}
int mid = (l + r) >> 1;
if(pos <= mid) return query_fa(ls[nd], 2 * o, l, mid, pos);
else return query_fa(rs[nd], 2 * o + 1, mid + 1, r, pos);
}
int query_size(int nd, int o, int l, int r, int pos) {
if(l == r) {
return size[nd];
}
int mid = (l + r) >> 1;
if(pos <= mid) return query_size(ls[nd], 2 * o, l, mid, pos);
else return query_size(rs[nd], 2 * o + 1, mid + 1, r, pos);
}
int Find_fa(int x, int M) {
int ff = query_fa(root[x], 1, 1, n, M);
if(ff == M) return ff;
int F = Find_fa(x, ff);
return F;
}
void Solve( ) {
int las = 0;
for(int i = 1;i <= m;i ++)
{
int opt,x,y;
scanf("%d",& opt);
if(opt == 1) {
scanf("%d%d",& x,& y); x = x + las, y = y + las;
int fa1 = Find_fa(i - 1, x), fa2 = Find_fa(i - 1, y);
if(fa1 == fa2) {
root[i] = root[i - 1]; continue;
}
int s1 = query_size(root[i - 1], 1, 1, n, fa1);
int s2 = query_size(root[i - 1], 1, 1, n, fa2);
if(s1 < s2) { swap(s1,s2); swap(fa1, fa2); }
root[i] = modify_fa(root[i - 1], 1, 1, n, fa2, fa1);
root[i] = modify_size(root[i], 1, 1, n, fa1, s2);
}
else {
scanf("%d%d",& x,& y); x = x + las, y = y + las;
root[i] = root[i - 1];
int fat = Find_fa(x, y);
int s = query_size(root[x], 1, 1, n, fat);
printf("%d\n",s); las = s;
}
}
}
int main( ) {
freopen("build.in","r",stdin);
freopen("build.out","w",stdout);
Init( );
Solve( );
}
这道题班上有人用贪心A掉了
而蒟蒻我只会垃圾dp qwqwq 首先我们可以发现这个玩意正着来搞是不行的
因为每个人一旦选择这座城市自己要 那么他的下一个人的选择是确定的
并且后面的选择跟前面的毫无联系 就无法体现最优选择的思想 所以我们考虑倒着搞 维护一个$sum$前缀和
$dp[i][0/1][0/1]$表示到了第$i$个城市 当前是谁的回合 这个人选不选择这个城市 这个人获得的最大收益
$dp[i][now][1] = sum[n] - sum[i] + a[i] - max(dp[i +1][now xor 1][0],dp[i + 1][now xor 1][1])$
表示如果我选择这座城市 相当于把主动权给别人 那么下个人一定会选择它的最优方案 总的收益是一定的 所以就用总收益减去别人的收益就可以了
$dp[i][now][0] = max(dp[i + 1][now][1],dp[i + 1][now][0])$
表示如果我这个位置不选 我的最大收益就是后面的选择所带来的最大收益
代码
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e5 + 5;
int n;
ll dp[N][2][2],sum[N],a[N];
void Init( ) {
scanf("%d",& n);
for(int i = 1;i <= n;i ++) {
scanf("%lld",& a[i]);
sum[i] = sum[i - 1] + a[i];
}
}
void Solve( ) {
dp[n][1][1] = a[n], dp[n][0][1] = a[n];
for(int i = n - 1;i >= 1;i --) {
for(int now = 0;now <= 1;now ++) {
dp[i][now][1] = a[i] + sum[n] - sum[i] - max(dp[i + 1][now ^ 1][0], dp[i + 1][now ^ 1][1]);
dp[i][now][0] = max(dp[i + 1][now][1], dp[i + 1][now][0]);
}
}
printf("%lld\n",max(dp[1][1][1], dp[1][1][0]));
}
int main( ) {
freopen("distribute.in","r",stdin);
freopen("distribute.out","w",stdout);
Init( );
Solve( );
}
然后我考试的时候做不来这道题
emmmmmm考虑要存在合法的圆边那么他肯定存在于一个环里面 在无向图里面 他肯定是一个点双连通分量
所以把每个点双处理出来 如果边数等于点数 他就是一个合法简单环 否则会有不止一个环 就不合法
考虑为什么是点双连通 不是边双呢 因为只有点双才是多个简单环叠加
左边这个虽然是边双 但是他是合法的 若是按照之前的方法 他会被判断成不合法 因为他不是简单环的叠加 而点双就一定是叠加 判断就一定合法
代码
#include <bits/stdc++.h>
using namespace std;
const int N = 1e5 + 5;
int n,head[N],nex[2 * N],tov[2 * N],id[2 * N],c[N];
int dfn[N],low[N],idc,tot = 1,m,stk[2 * N],top,col,que[2 * N];
bool vis[2 * N],isans[2 * N];
void add(int u, int v, int ID) {
tot ++;
nex[tot] = head[u];
tov[tot] = v; id[tot] = ID;
head[u] = tot;
}
void Init( ) {
scanf("%d%d",& n,& m);
for(int i = 1;i <= m;i ++) {
int u,v;
scanf("%d%d",& u,& v);
add(u, v, i); add(v, u, i);
}
}
void tarjan(int u,int from) {
low[u] = dfn[u] = ++idc;
for(int i = head[u];i;i = nex[i]) {
int v = tov[i];
if((i ^ 1) == from) continue;
if(! vis[i]) vis[i] = vis[i ^ 1] = true, stk[++ top] = i;
if(! dfn[v]) {
tarjan(v, i);
low[u] = min(low[u], low[v]);
if(low[v] >= dfn[u]) {
col ++;
int cntn = 0,cnte = 0;
while(1) {
int e = stk[top --];
if(c[tov[e]] != col) c[tov[e]] = col,cntn ++;
if(c[tov[e ^ 1]] != col) c[tov[e ^ 1]] = col, cntn ++;
que[++ cnte] = e;
if(e == i) break;
}
if(cnte == cntn) {
for(int i = 1;i <= cnte;i ++)
isans[que[i]] = isans[que[i] ^ 1] = true;
}
}
}
else if(low[u] > dfn[v]) low[u] = min(low[u], dfn[v]);
}
}
void Solve( ) {
for(int i = 1;i <= n;i ++) {
if(! dfn[i]) tarjan(i, 0);
}
int ans = 0;
for(int i = 2;i <= tot;i += 2) if(isans[i]) ans ++;
printf("%d\n",ans);
for(int i = 2;i <= tot;i += 2)
if(isans[i]) printf("%d ",i / 2);
}
int main( ) {
freopen("find.in","r",stdin);
freopen("find.out","w",stdout);
Init( );
Solve( );
}
原文地址：https://www.cnblogs.com/Rubenisveryhandsome/p/9690596.html
时间： 2024-08-29 17:53:15