环套树

  • P1453
  • 先把环处理出来,然后对环上每个节点树形DP,然后环形DP就好
//#define LOCAL
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define mem(a, b) memset(a,b,sizeof(a))
#define sz(a) (int)a.size()
#define INF 0x3f3f3f3f
#define DNF 0x7f
#define DBG printf("this is a input\n")
#define fi first
#define se second
#define mk(a, b) make_pair(a,b)
#define pb push_back
#define LF putchar('\n')
#define SP putchar(' ')
#define p_queue priority_queue
#define CLOSE ios::sync_with_stdio(0); cin.tie(0)

template<typename T>
void read(T &x) {x = 0;char ch = getchar();ll f = 1;while(!isdigit(ch)){if(ch == '-')f *= -1;ch = getchar();}while(isdigit(ch)){x = x * 10 + ch - 48; ch = getchar();}x *= f;}
template<typename T, typename... Args>
void read(T &first, Args& ... args) {read(first);read(args...);}
template<typename T>
void write(T arg) {T x = arg;if(x < 0) {putchar('-'); x =- x;}if(x > 9) {write(x / 10);}putchar(x % 10 + '0');}
template<typename T, typename ... Ts>
void write(T arg, Ts ... args) {write(arg);if(sizeof...(args) != 0) {putchar(' ');write(args ...);}}
using namespace std;

ll gcd(ll a, ll b) {
    return b == 0 ? a : gcd(b, a % b);
}

ll lcm(ll a, ll b) {
    return a / gcd(a, b) * b;
}

const int N = 200005;
int n, p[N], head[N], cnt;
double k;
struct node
{
    int next, t;
}edge[N << 1];
void add (int f, int t)
{
    edge[cnt].t = t;
    edge[cnt].next = head[f];
    head[f] = cnt ++;
}
int dfn[N], path[N], times;
int vertex[N], circle[N], tot, dp[N][2], f[N][2];
void dfs1 (int u ,int id)
{
    dfn[u] = ++ times;
    for(int i = head[u] ; i != -1; i = edge[i].next)
    {
        int v = edge[i].t;
        if(i != id)
        {
            if (!dfn[v])
                path[v] = u, dfs1(v,i^1);//前向星选手都懂
            else if (dfn[v] > dfn[u])//tarjan原理
            {
                for(int j = v ; j != u ; j = path[j])
                    vertex[++ tot] = j , circle[j] = 1;
                vertex[++ tot] = u , circle[u] = 1;
            }
        }
    }
}
void dfs2 (int u, int fa)
{
    dp[u][1] = p[u];
    for (int i = head[u] ; i != -1 ; i = edge[i].next)
    {
        int v = edge[i].t;
        if (v != fa && !circle[v])
        {
            dfs2 (v, u);
            dp[u][1] += dp[v][0];
            dp[u][0] += max(dp[v][1], dp[v][0]);
        }
    }
}
int main()
{
    read (n);
    mem(head,-1);
    for (int i = 0 ; i < n ; i ++)
        read (p[i]);
    for (int i = 1 ; i <= n ; i ++)
    {
        int u , v;
        read (u, v);
        add (u , v);
        add (v , u);
    }
    scanf ("%lf",&k);
    dfs1 (1, -1);
    for (int i = 1 ; i <= tot ; i ++)
        dfs2(vertex[i], -1);
    for (int i = 1 ; i <= tot ; i ++)
        vertex[i + tot] = vertex[i];
    int ans = 0;
    mem(f,-0x3f);
    f[1][0] = dp[vertex[1]][0];
    for (int i = 2 ; i <= tot ; i ++)
    {
        f[i][1] = f[i-1][0]+dp[vertex[i]][1];
        f[i][0] = max(f[i-1][0],f[i-1][1])+dp[vertex[i]][0];
    }
    ans = max (f[tot][0], f[tot][1]);
    mem(f,-0x3f);
    f[1][1] = dp[vertex[1]][1];
    for (int i = 2 ; i <= tot ; i ++)
    {
        f[i][1] = f[i-1][0] + dp[vertex[i]][1];
        f[i][0] = max(f[i-1][0],f[i-1][1]) + dp[vertex[i]][0];
    }
    ans = max(ans, f[tot][0]);
    printf("%.1lf\n",ans*k);
}
/*
4
1 2 2 3
0 1 3 1
1 2 3 2
 */
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
r语言可以使用ggtree和ggtreeExtra包来绘制环形进化。首先,你需要准备好的数据和注释数据。然后,可以按照以下步骤进行绘图: 1. 绘制的主体,可以选择特殊的布局(如equal_angle),并为枝添加一些分类颜色信息。 2. 添加外圈注释1,即标签,将标签放置在每个类别分支附近,并为背景颜色使用分类信息。 3. 添加外圈注释2,即点和文字,手动选择一些节点,在枝顶端添加灰色的点和黑色的文字。 下面是一个示例代码,演示如何使用ggtree和ggtreeExtra包绘制环形进化: ```R # 准备数据 tree_df <- read.tree("tree_file.txt") # 替换为你的文件 phy_nodes <- c("A", "B", "C") # 替换为你选择的分支节点 label_node <- c("D", "E", "F") # 替换为你选择的节点 # 绘制的主体 p <- ggtree(tree_df, aes(color = Form), layout = "equal_angle") + geom_treescale(-5, 7, fontsize = 3, linesize = 0.5, width = 1) + scale_color_manual(values = c("black", colors)) + coord_flip() + theme(legend.position = "none") # 添加外圈注释1 p1 <- p + geom_label_repel(data = subset(tree_df, node %in% phy_nodes), mapping = aes(x = x, y = y, label = Form, fill = Form), color = "black", alpha = 0.7) + scale_fill_manual(values = colors) # 添加外圈注释2 p2 <- p1 + geom_point(data = subset(tree_df, node %in% label_node), mapping = aes(x = x * 1.03, y = y * 1.03), color = "grey50") + geom_text_repel(data = subset(tree_df, node %in% label_node), mapping = aes(x = x * 1.05, y = y * 1.05, label = label), color = "black") # 显示绘制结果 print(p2) ``` 希望这个例子能帮助到你!如果你还有其他问题,请随时提问。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值