HYSBZ3531-旅行

3531: [Sdoi2014]旅行

Time Limit: 40 Sec   Memory Limit: 512 MB
Submit: 2451   Solved: 1067
[ Submit][ Status][ Discuss]

Description

 S国有N个城市,编号从1到N。城市间用N-1条双向道路连接,满足
从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教,如飞天面条神教、隐形独角兽教、绝地教都是常见的信仰。为了方便,我们用不同的正整数代表各种宗教,  S国的居民常常旅行。旅行时他们总会走最短路,并且为了避免麻烦,只在信仰和他们相同的城市留宿。当然旅程的终点也是信仰与他相同的城市。S国政府为每个城市标定了不同的旅行评级,旅行者们常会记下途中(包括起点和终点)留宿过的城市的评级总和或最大值。
    在S国的历史上常会发生以下几种事件:
”CC x c”:城市x的居民全体改信了c教;
”CW x w”:城市x的评级调整为w;
”QS x y”:一位旅行者从城市x出发,到城市y,并记下了途中留宿过的城市的评级总和;
”QM x y”:一位旅行者从城市x出发,到城市y,并记下了途中留宿过的城市的评级最大值。
    由于年代久远,旅行者记下的数字已经遗失了,但记录开始之前每座城市的信仰与评级,还有事件记录本身是完好的。请根据这些信息,还原旅行者记下的数字。    为了方便,我们认为事件之间的间隔足够长,以致在任意一次旅行中,所有城市的评级和信仰保持不变。

Input

    输入的第一行包含整数N,Q依次表示城市数和事件数。
    接下来N行,第i+l行两个整数Wi,Ci依次表示记录开始之前,城市i的
评级和信仰。
    接下来N-1行每行两个整数x,y表示一条双向道路。
    接下来Q行,每行一个操作,格式如上所述。

Output

    对每个QS和QM事件,输出一行,表示旅行者记下的数字。

Sample Input

5 6
3 1
2 3
1 2
3 3
5 1
1 2
1 3
3 4
3 5
QS 1 5
CC 3 1
QS 1 5
CW 3 3
QS 1 5
QM 2 4

Sample Output

8
9
11
3

HINT

N,Q < =10^5    , C < =10^5


 数据保证对所有QS和QM事件,起点和终点城市的信仰相同;在任意时

刻,城市的评级总是不大于10^4的正整数,且宗教值不大于C。

Source



解题思路:树链剖分,对每种宗教都建一棵动态树


#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <queue>
#include <vector>
#include <bitset>
#include <functional>

using namespace std;

#define LL long long
const int INF = 0x3f3f3f3f;

int n, q, x, y, tot;
int a[100009], b[100009], g[100009], top[100009];
int s[100009], nt[200009], e[200009], cnt;
int ct[100009], mx[100009], fa[100009], dep[100009];
int ss[100009], L[100009 * 20], R[100009 * 20];
LL sum[100009 * 20], ma[100009 * 20];
char ch[5];

void dfs(int k, int f)
{
	dep[k] = dep[f] + 1;
	fa[k] = f, ct[k] = 1, mx[k] = 0;
	for (int i = s[k]; ~i; i = nt[i])
	{
		if (e[i] == f) continue;
		dfs(e[i], k);
		ct[k] += ct[e[i]];
		if (ct[e[i]] > ct[mx[k]]) mx[k] = e[i];
	}
}

void update(int &k, int l, int r, int p, int val)
{
	if (!k) { k = ++tot, L[tot] = R[tot] = 0; }
	if (l == r) { sum[k] = ma[k] = val; return; }
	int mid = (l + r) >> 1;
	if (mid >= p) update(L[k], l, mid, p, val);
	else update(R[k], mid + 1, r, p, val);
	ma[k] = max(ma[L[k]], ma[R[k]]);
	sum[k] = sum[L[k]] + sum[R[k]];
}

void Dfs(int k, int t)
{
	top[k] = !t ? k : top[fa[k]];
	g[k] = ++cnt;
	update(ss[b[k]], 1, n, g[k], a[k]);
	if (mx[k]) Dfs(mx[k], 1);
	for (int i = s[k]; ~i; i = nt[i])
	{
		if (e[i] == fa[k] || e[i] == mx[k]) continue;
		Dfs(e[i], 0);
	}
}

LL query(int k, int l, int r, int ll, int rr, int flag)
{
	if (l >= ll&&r <= rr) return flag ? sum[k] : ma[k];
	int mid = (l + r) >> 1;
	LL ans = 0;
	if (ll <= mid)
	{
		if (flag) ans += query(L[k], l, mid, ll, rr, flag);
		else ans = max(ans, query(L[k], l, mid, ll, rr, flag));
	}
	if (rr > mid)
	{
		if (flag) ans += query(R[k], mid + 1, r, ll, rr, flag);
		else ans = max(ans, query(R[k], mid + 1, r, ll, rr, flag));
	}
	return ans;
}

LL solve(int x, int y, int flag)
{
	LL ans = 0;
	int k = ss[b[x]];
	while (top[x] != top[y])
	{
		if (dep[top[x]] < dep[top[y]]) swap(x, y);
		if (flag) ans += query(k, 1, n, g[top[x]], g[x], flag);
		else ans = max(ans, query(k, 1, n, g[top[x]], g[x], flag));
		x = fa[top[x]];
	}
	if (dep[x] > dep[y]) swap(x, y);
	if (flag) ans += query(k, 1, n, g[x], g[y], flag);
	else ans = max(ans, query(k, 1, n, g[x], g[y], flag));
	return ans;
}

int main()
{
	//freopen("input.txt", "r", stdin);
	//freopen("output.txt", "w", stdout);
	while (~scanf("%d%d", &n, &q))
	{
		for (int i = 1; i <= n; i++) scanf("%d%d", &a[i], &b[i]);
		memset(s, -1, sizeof s);
		memset(ss, 0, sizeof ss);
		dep[0] = ct[0] = cnt = 0;
		int u, v;
		for (int i = 1; i < n; i++)
		{
			scanf("%d%d", &u, &v);
			nt[cnt] = s[u], s[u] = cnt, e[cnt++] = v;
			nt[cnt] = s[v], s[v] = cnt, e[cnt++] = u;
		}
		dfs(1, 0);
		Dfs(1, cnt = tot = 0);
		while (q--)
		{
			scanf("%s%d%d", ch, &x, &y);
			if (ch[1] == 'C')
			{
				update(ss[b[x]], 1, n, g[x], 0);
				update(ss[b[x] = y], 1, n, g[x], a[x]);
			}
			else if (ch[1] == 'W') update(ss[b[x]], 1, n, g[x], a[x] = y);
			else printf("%lld\n", solve(x, y, ch[1] == 'S'));
		}
	}
	return 0;
}

内容概要:本文详细介绍了扫描单分子定位显微镜(scanSMLM)技术及其在三维超分辨体积成像中的应用。scanSMLM通过电调透镜(ETL)实现快速轴向扫描,结合4f检测系统将不同焦平面的荧光信号聚焦到固定成像面,从而实现快速、大视场的三维超分辨成像。文章不仅涵盖了系统硬件的设计与实现,还提供了详细的软件代码实现,包括ETL控制、3D样本模拟、体积扫描、单分子定位、3D重建和分子聚类分析等功能。此外,文章还比较了循环扫描与常规扫描模式,展示了前者在光漂白效应上的优势,并通过荧光珠校准、肌动蛋白丝、线粒体网络和流感A病毒血凝素(HA)蛋白聚类的三维成像实验,验证了系统的性能和应用潜力。最后,文章深入探讨了HA蛋白聚类与病毒感染的关系,模拟了24小时内HA聚类的动态变化,提供了从分子到细胞尺度的多尺度分析能力。 适合人群:具备生物学、物理学或工程学背景,对超分辨显微成像技术感兴趣的科研人员,尤其是从事细胞生物学、病毒学或光学成像研究的科学家和技术人员。 使用场景及目标:①理解和掌握scanSMLM技术的工作原理及其在三维超分辨成像中的应用;②学习如何通过Python代码实现完整的scanSMLM系统,包括硬件控制、图像采集、3D重建和数据分析;③应用于单分子水平研究细胞内结构和动态过程,如病毒入侵机制、蛋白质聚类等。 其他说明:本文提供的代码不仅实现了scanSMLM系统的完整工作流程,还涵盖了多种超分辨成像技术的模拟和比较,如STED、GSDIM等。此外,文章还强调了系统在硬件改动小、成像速度快等方面的优势,为研究人员提供了从理论到实践的全面指导。
内容概要:本文详细介绍了基于Seggiani提出的渣层计算模型,针对Prenflo气流床气化炉中炉渣的积累和流动进行了模拟。模型不仅集成了三维代码以提供气化炉内部的温度和浓度分布,还探讨了操作条件变化对炉渣行为的影响。文章通过Python代码实现了模型的核心功能,包括炉渣粘度模型、流动速率计算、厚度更新、与三维模型的集成以及可视化展示。此外,还扩展了模型以考虑炉渣组成对特性的影响,并引入了Bingham流体模型,更精确地描述了含未溶解颗粒的熔渣流动。最后,通过实例展示了氧气-蒸汽流量增加2%时的动态响应,分析了温度、流动特性和渣层分布的变化。 适合人群:从事煤气化技术研究的专业人士、化工过程模拟工程师、以及对工业气化炉操作优化感兴趣的科研人员。 使用场景及目标:①评估不同操作条件下气化炉内炉渣的行为变化;②预测并优化气化炉的操作参数(如温度、氧煤比等),以防止炉渣堵塞;③为工业气化炉的设计和操作提供理论支持和技术指导。 其他说明:该模型的实现基于理论公式和经验数据,为确保模型准确性,实际应用中需要根据具体气化炉的数据进行参数校准。模型还考虑了多个物理场的耦合,包括质量、动量和能量守恒方程,能够模拟不同操作条件下的渣层演变。此外,提供了稳态求解器和动态模拟工具,可用于扰动测试和工业应用案例分析。
好的,你可以参考以下示例代码实现Java销售额查询和营业额统计: ```java import java.util.ArrayList; import java.util.HashMap; import java.util.List; import java.util.Map; public class SalesStatistics { private Map<String, Double> salesData = new HashMap<>(); public SalesStatistics() { // 初始化销售数据 salesData.put("2021-01-01", 1000.0); salesData.put("2021-01-02", 1500.0); salesData.put("2021-01-03", 2000.0); salesData.put("2021-01-04", 1200.0); salesData.put("2021-01-05", 1800.0); } // 查询指定日期的销售额 public double getSalesByDate(String date) { if (salesData.containsKey(date)) { return salesData.get(date); } else { return 0.0; } } // 统计指定日期范围内的营业额 public double getTotalSales(String startDate, String endDate) { double totalSales = 0.0; for (String date : salesData.keySet()) { if (date.compareTo(startDate) >= 0 && date.compareTo(endDate) <= 0) { totalSales += salesData.get(date); } } return totalSales; } // 获取所有销售日期 public List<String> getAllSalesDates() { return new ArrayList<>(salesData.keySet()); } } ``` 使用示例: ```java public static void main(String[] args) { SalesStatistics salesStatistics = new SalesStatistics(); // 查询指定日期的销售额 double sales = salesStatistics.getSalesByDate("2021-01-02"); System.out.println("2021-01-02的销售额为:" + sales); // 统计指定日期范围内的营业额 double totalSales = salesStatistics.getTotalSales("2021-01-02", "2021-01-04"); System.out.println("2021-01-02到2021-01-04的营业额为:" + totalSales); // 获取所有销售日期 List<String> salesDates = salesStatistics.getAllSalesDates(); System.out.println("所有销售日期为:" + salesDates); } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值