HDU 4812 D Tree

HDU 4812 

思路:

点分治

先预处理好1e6 + 3以内到逆元

然后用map 映射以分治点为起点的链的值a 成他的下标 u 

然后暴力跑出以分治点儿子为起点的链的值b,然后在map里查找inv[b]*k

代码:

#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pi acos(-1.0)
#define LL long long
//#define mp make_pair
#define pb push_back
#define ls rt<<1, l, m
#define rs rt<<1|1, m+1, r
#define ULL unsigned LL
#define pll pair<LL, LL>
#define pii pair<int, int>
#define mem(a, b) memset(a, b, sizeof(a))
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define fopen freopen("in.txt", "r", stdin);freopen("out.txt", "w", stout);
//headconst int MOD = 1e6 + 3;
const int INF = 0x7f7f7f7f;
const int N = 1e5 + 5;
int inv[MOD + 5], mp[MOD + 5], head[N], mxsz[N], sz[N], v[N], cnt = 0, rt = 0, n, k, ans1, ans2;
int deep[N], dis[N], id[N], top = 0;
bool vis[N];
struct edge {int to, nxt;
}edge[N*2];
void add_edge(int u, int v) {edge[cnt].to = v;edge[cnt].nxt = head[u];head[u] = cnt++;
}
void init() {inv[1] = 1;for (int i = 2; i < MOD; i++) inv[i] = (MOD - MOD/i) * 1LL * inv[MOD%i] % MOD;
}
void update(int x, int y) {int t = (1LL * inv[x] * k) % MOD;int now = mp[t];if(!now) return ;if(now > y) swap(now, y);if(now < ans1 || now == ans1 && y < ans2) ans1 = now, ans2 = y;
}
void get_rt(int o, int u) {sz[u] = 1, mxsz[u] = 0;for (int i = head[u]; ~i; i = edge[i].nxt) {if(edge[i].to != o && !vis[edge[i].to]) {get_rt(u, edge[i].to);sz[u] += sz[edge[i].to];mxsz[u] = max(mxsz[u], sz[edge[i].to]);}}mxsz[u] = max(mxsz[u], n - sz[u]);if(mxsz[u] < mxsz[rt]) rt = u;
}
void get_d(int o, int u) {deep[++top] = dis[u];id[top] = u;for (int i = head[u]; ~i; i = edge[i].nxt) {if(!vis[edge[i].to] && edge[i].to != o) {dis[edge[i].to] = (1LL * dis[u] * v[edge[i].to])%MOD;get_d(u, edge[i].to);}}
}
void solve(int u) {vis[u] = true;mp[v[u]] = u;for (int i = head[u]; ~i; i = edge[i].nxt) {if(!vis[edge[i].to]) {top = 0, dis[edge[i].to] = v[edge[i].to];get_d(u, edge[i].to);for (int j = 1; j <= top; j++) update(deep[j], id[j]);top = 0, dis[edge[i].to] = (1LL * v[u] * v[edge[i].to])%MOD;get_d(u, edge[i].to);for (int j = 1; j <= top; j++) {int t = deep[j];if(!mp[t] || id[j] < mp[t]) mp[t] = id[j];}}}mp[v[u]] = 0;for (int i = head[u]; ~i; i = edge[i].nxt) {if(!vis[edge[i].to]) {top = 0, dis[edge[i].to] = (1LL * v[u] * v[edge[i].to])%MOD;get_d(u, edge[i].to);for (int j = 1; j <= top; j++) mp[deep[j]] = 0;}}for (int i = head[u]; ~i; i = edge[i].nxt) {if(!vis[edge[i].to]) {mxsz[0] = n =  sz[edge[i].to];get_rt(rt = 0, edge[i].to);solve(rt);}}
}
int main() {init();int u, V;while(~scanf("%d%d", &n, &k)) {mem(head, -1);mem(vis, false);mem(mp, 0);cnt = 0;ans1 = ans2 = INF;for (int i = 1; i <= n; i++) scanf("%d", &v[i]);for (int i = 1; i < n; i++) scanf("%d%d", &u, &V), add_edge(u, V), add_edge(V, u);mxsz[0] = n;get_rt(rt = 0, 1);solve(rt);if(ans1 == INF) printf("No solution\n");else printf("%d %d\n", ans1, ans2);}return 0;
}

 

转载于:https://www.cnblogs.com/widsom/p/9106507.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/538564.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Angular CLI 安装

安装Angular 官网的教程&#xff0c;因为国内网络环境原因&#xff0c;访问不了服务器&#xff0c;导致安装失败。 1、先安装NodeJs 安装教程&#xff1a;http://blog.csdn.net/zengmingen/article/details/72650484 2、通过NodeJs中的模块npm 命令行安装 CLI 2.1、设置npm的…

go 写文件_「go」 项目多个文件编程

golang 学习的时候很多sample 讲的都是一个文件的go 文件怎么写&#xff0c;但是现实中不可能所有的实现都写到一个文件里面&#xff0c;按照功能的不同&#xff0c;要么拆分成不同的文件&#xff0c;要么拆分成不同的文件。下面有些个人的经验分享下&#xff0c;如果有问题请指…

CycleGAN 各种变变变

转载自 简单介绍了一下GAN和DCGAN的原理。以及如何使用Tensorflow做一个简单的生成图片的demo。 Ian Goodfellow对GAN一系列工作总结的ppt&#xff0c;确实精彩&#xff0c;推荐&#xff1a;独家 | GAN之父NIPS 2016演讲现场直击&#xff1a;全方位解读生成对抗网络的原理及未来…

pycharm与webstorm 2017 激活破解

原有的方式已经失效&#xff0c;见下面博文&#xff1a; https://blog.csdn.net/justszh/article/details/81484802

mysql blob 比较_与MSSQL对比学习MYSQL的心得(四)--BLOB数据类型

MYSQL里的BLOB数据类型BLOB是一个二进制大对象&#xff0c;用来存储可变数量的数据。BLOB类型分为4种&#xff1a;TinyBlob、Blob、MediumBlob、LongBlob&#xff0c;这几个类型之间的唯一区别是在存储文件的最大大小上不同。MySQL的四种BLOB类型 类型 大小(单位&#xff1…

Webstorm常用快捷键

webstrom 使用 eclipse快键键 File--settings keymap 选择 eclipse 原文链接&#xff1a;http://www.cnblogs.com/yeminglong/p/5995421.html ------------------以下是webstrom默认的----------------------------------- Ctrl/ 或 CtrlShift/ 注释&#xff08;// 或者/…

VirtualBox 上安装Debian 后分辨率设置

VirtualBox 上安装Debian 后分辨率设置 首先要配置source.list打开终端&#xff0c; su 切换成root用户&#xff0c; cd /etc/apt 然后编辑source.list rootdebian:/etc/apt# vi source.list 注释deb cdrom:行&#xff0c;加以下源 deb http://deb.debian.org/debian stretc…

疯狂的程序员_程序员的乐趣是什么?

作者&#xff1a;Java3y我是一个程序员&#xff0c;外行人都以为我是修电脑的&#xff0c;我笑了笑&#xff0c;随意ctrl cctrl v了一把&#xff0c;想象着你们因为我的文章而开心不止&#xff0c;我感到充实而欣慰。想象着你们给我拼命点赞的样子&#xff0c;是多么的滑稽&…

template多行编写的方式

模板是包在 ECMAScript 2015 反引号 () 中的一个多行字符串。 反引号 () — 注意&#xff0c;不是单引号 () — 允许把一个字符串写在多行上&#xff0c; 使 HTML 模板更容易阅读。 反引号&#xff1a;键盘数字键1 旁边的&#xff0c;ESC键下面的键 如果单引号 Component({sel…

sqllite事务和MySQL事务_Android学习---SQLite数据库的增删改查和事务(transaction)调用...

上一篇文章中介绍了手工拼写sql语句进行数据库的CRUD操作,本文将介绍调用sqlite内置的方法实现CRUD操作,其实质也是通过拼写sql语句.首先,创建一个新的android项目:其次,查看代码实现增删查改:1.创建DB工具类MyDBHelper.java(创建数据库的操作)packagecom.amos.android_db;impo…

sqlserver2000给账户授予所有的权限_你的位置信息权限设置对了么?

位置信息权限是众多应用权限中的一种&#xff0c;是应用获取手机地理位置信息的必要凭证。在你首次安装应用并打开时&#xff0c;通常会出现一连串的权限弹框&#xff0c;如果该应用在其运行过程中会用到你的地理位置信息&#xff0c;那么这些弹框中就会包含一个与位置信息有关…

Python之路,Day1 - Python基础1

本节内容 Python介绍发展史Python 2 or 3?安装Hello World程序变量用户输入模块初识.pyc是个什么鬼&#xff1f;数据类型初识数据运算表达式if ...else语句表达式for 循环break and continue 表达式while 循环作业需求 一、 Python介绍 python的创始人为吉多范罗苏姆&#xf…

mysql 范式化_MySQL-范式和反范式

1.第一范式(1NF)(列不能再拆分)原子性&#xff0c;字段不可分(列的信息)&#xff0c;只要是关系型数据库&#xff0c;就自动满足1NF&#xff1b;2.第二范式(2NF)(主键唯一&#xff0c;且被依赖)在第一范式基础上建立的&#xff0c;即满足第二范式的必须先满足第一范式。要求DB表…

端口被占用解决办法

1. 端口被占用解决办法 netstat -ano | findstr 8080(端口号) taskkill -pid (进程pid) –f转载于:https://www.cnblogs.com/xaoco/p/9114773.html

java 判断是否是list_JAVA从头开始一基础梳理(4-3)

大家好&#xff0c;今天我们介绍一下java中常用的集合类型。首先&#xff0c;我们先看一下java中集合类型的结构。以上是集合的继承关系图&#xff0c;通常我们使用的比较多的是 Set , List , Map以及其衍生的子类和接口实现类。首先给大家介绍一下List&#xff0c;List本身是一…

Python2.x还是3.x?

2.x 和 3.x对于程序员的编码来说&#xff0c;没有发生太大的变化&#xff0c;当然也是有变化的&#xff0c;主要是Python内部发生了巨变。 要用3.x的原因是&#xff1a; 1、3.x和2.x版本不兼容。 2、Python库新增的内容不支持2.x了。 3、2.x版本官方支持到2020年结束。 晚改…

前端网页广告无线翻滚_从小白到web前端工程师进阶之路 从0到1到更深

互联网的发展&#xff0c;让web前端技术发生了翻天覆地的变化&#xff0c;前端开发工程师可以让网页内容变得更加生动&#xff0c;为用户带来更好的体验。那么&#xff0c;武汉web前端培训哪个好&#xff1f;web前端好学吗&#xff1f;作为一个合格的Web前端工程师&#xff0c;…

PowerDesigner导出表为Excel(转)

打开脚本运行器CtrlShiftX 导出&#xff1a; ****************************************************************************** Option ExplicitDim rowsNumrowsNum 0 -----------------------------------------------------------------------------Main function -------…

判读一个对象不为空_ArrayList实现分析(一)——对象创建

ArrayList是java中最常用的集合类之一&#xff0c;它的内部实现是基于数组&#xff0c;因此ArryList可以根据索引实现随机访问。ArryList继承了AbstractList类&#xff0c;并且实现了List, RandomAccess, Cloneable接口。下面详细分析一下ArrayList的实现&#xff0c;下面的分析…

AngularJS与Angular的区别

指同一事物&#xff0c;版本的区别&#xff0c;叫法不同 Angular2.0之前的版本&#xff08;1.x&#xff09;叫做AngularJS 1.x的使用是引入AngularJS的js文件到网页。 2.0之后&#xff0c;就是完全不同了。 Angular2.x与Angular1.x 的区别类似 Java 和 JavaScript 或者说是…