HDU - 4812 D Tree 点分治
http://acm.hdu.edu.cn/showproblem.php?pid=4812
题意:有一棵树,每个点有一个权值要求找最小的一对点,路径上的乘积mod1e6+3为k
题解:点分治,挨个把子树更新,每次把子树和现有的map里找满足条件的点对,然后更新子树到map里,map维护的是每个到根的乘积值最小的点.
//#pragma comment(linker, "/stack:200000000")
//#pragma GCC optimize("Ofast,no-stack-protector")
//#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")
//#pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
#define fi first
#define se second
#define mp make_pair
#define pb push_back
#define pi acos(-1.0)
#define ll long long
#define vi vector<int>
#define mod 1000003
#define C 0.5772156649
#define ls l,m,rt<<1
#define rs m+1,r,rt<<1|1
#define pil pair<int,ll>
#define pli pair<ll,int>
#define pii pair<int,int>
#define cd complex<double>
#define ull unsigned long long
#define base 1000000000000000000
#define fio ios::sync_with_stdio(false);cin.tie(0)
using namespace std;
const double eps=1e-6;
const int N=100000+10,maxn=1000000+10,inf=0x3f3f3f3f,INF=0x3f3f3f3f3f3f3f3f;
struct edge{
int to,Next;
}e[N*2];
int cnt,head[N];
void add(int u,int v)
{
e[cnt].to=v;
e[cnt].Next=head[u];
head[u]=cnt++;
}
int n,k,inv[maxn];
bool vis[N];
int sz[N],zx[N],a[N];
map<int,int>ma;
int x,y;
vector<pair<ll,int> >te;
void init()
{
inv[1]=1;
for(int i=2;i<maxn;i++)
inv[i]=(mod-mod/i)*1ll*inv[mod%i]%mod;
}
void dfssz(int u,int f)
{
sz[u]=1;
for(int i=head[u];~i;i=e[i].Next)
{
int x=e[i].to;
if(x!=f && !vis[x])
{
dfssz(x,u);
sz[u]+=sz[x];
}
}
}
void dfszx(int u,int f,int root,int &ans)
{
zx[u]=sz[root]-sz[u];
for(int i=head[u];~i;i=e[i].Next)
{
int x=e[i].to;
if(x!=f && !vis[x])
{
zx[u]=max(zx[u],sz[x]);
dfszx(x,u,root,ans);
}
}
if(zx[u]<zx[ans])ans=u;
}
int findzx(int root)
{
dfssz(root,-1);
int ans=root;
zx[ans]=inf;
dfszx(root,-1,root,ans);
return ans;
}
void dfs1(int u,int f,int root,ll ans)
{
if(a[root]*ans%mod==k)
{
int xx=root,yy=u;
if(xx>yy)swap(xx,yy);
if(xx<x || (xx==x&&yy<y))x=xx,y=yy;
}
ll te=1ll*k*inv[ans*a[root]%mod]%mod;
if(ma.find(te)!=ma.end())
{
int xx=ma[te],yy=u;
if(xx>yy)swap(xx,yy);
if(xx<x || (xx==x&&yy<y))x=xx,y=yy;
}
for(int i=head[u];~i;i=e[i].Next)
{
int x=e[i].to;
if(x!=f && !vis[x])dfs1(x,u,root,ans*a[x]%mod);
}
}
void dfs2(int u,int f,ll ans)
{
if(ma.find(ans)==ma.end())ma[ans]=u;
else if(u<ma[ans])ma[ans]=u;
for(int i=head[u];~i;i=e[i].Next)
{
int x=e[i].to;
if(x!=f && !vis[x])dfs2(x,u,ans*a[x]%mod);
}
}
void solve(int root)
{
int zx=findzx(root);
ma.clear();
vis[zx]=1;
for(int i=head[zx];~i;i=e[i].Next)
{
int x=e[i].to;
if(!vis[x])
{
te.clear();
dfs1(x,zx,zx,a[x]);
dfs2(x,zx,a[x]);
}
}
for(int i=head[zx];~i;i=e[i].Next)
if(!vis[e[i].to])
solve(e[i].to);
}
int main()
{
init();
while(~scanf("%d%d",&n,&k))
{
cnt=0;
memset(head,-1,sizeof head);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
vis[i]=0;
}
for(int i=1;i<n;i++)
{
int a,b;
scanf("%d%d",&a,&b);
add(a,b);add(b,a);
}
x=n+1,y=n+1;
solve(1);
if(x==n+1)puts("No solution");
else printf("%d %d\n",x,y);
}
return 0;
}
/********************
5 5
2 5 2 3 3
1 2
1 3
2 4
2 5
********************/
HDU - 4812 D Tree 点分治的更多相关文章
- HDU 4812 D Tree 树分治+逆元处理
D Tree Problem Description There is a skyscraping tree standing on the playground of Nanjing Unive ...
- HDU 4812 D Tree 树分治
题意: 给出一棵树,每个节点上有个权值.要找到一对字典序最小的点对\((u, v)(u < v)\),使得路径\(u \to v\)上所有节点权值的乘积模\(10^6 + 3\)的值为\(k\) ...
- HDU 4812 D Tree
HDU 4812 思路: 点分治 先预处理好1e6 + 3以内到逆元 然后用map 映射以分治点为起点的链的值a 成他的下标 u 然后暴力跑出以分治点儿子为起点的链的值b,然后在map里查找inv[b ...
- hdu 4812 D Tree(树的点分治)
D Tree Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 102400/102400 K (Java/Others) Total ...
- HDU 4812 D Tree 树分区+逆+hash新位置
意甲冠军: 特定n点树 K 以下n号码是正确的点 以下n-1行给出了树的侧. 问: 所以,如果有在正确的道路点图的路径 % mod = K 如果输出路径的两端存在. 多条路径则输出字典序最小的一条. ...
- HDU 4871 Shortest-path tree 最短路 + 树分治
题意: 输入一个带权的无向连通图 定义以顶点\(u\)为根的最短路生成树为: 树上任何点\(v\)到\(u\)的距离都是原图最短的,如果有多条最短路,取字典序最小的那条. 然后询问生成树上恰好包含\( ...
- hdu 4812 DTree (点分治)
D Tree Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 102400/102400 K (Java/Others)Total S ...
- HDU 5513 Efficient Tree
HDU 5513 Efficient Tree 题意 给一个\(N \times M(N \le 800, M \le 7)\)矩形. 已知每个点\((i-1, j)\)和\((i,j-1)\)连边的 ...
- hdu 3842 Machine Works(cdq分治维护凸壳)
题目链接:hdu 3842 Machine Works 详细题解: HDU 3842 Machine Works cdq分治 斜率优化 细节比较多,好好体会一下. 在维护斜率的时候要考虑x1与x2是否 ...
随机推荐
- java个人所得税计算器
class Caculate{ private String name; private double money; private double actual; /** * @param usern ...
- hibernate基于注解实现映射关系的配置
关联关系的配置步骤 ①要理清楚管理关系 ②确定管理依赖关系的哪一方 1一对一例如:person 和IdCard ①确定依赖关系:一对一 ②依赖关系由person类管理代码如下: person: @En ...
- Code Forces 652A Gabriel and Caterpillar
A. Gabriel and Caterpillar time limit per test 1 second memory limit per test 256 megabytes input st ...
- c# 自定义控件之小小进度条
先看效果图: 非常简洁的一个进度条. 完整项目源码下载:http://files.cnblogs.com/files/tuzhiyuan/%E8%BF%9B%E5%BA%A6%E6%9D%A1%E6% ...
- 知乎live 我的读书经验 总结
https://www.zhihu.com/lives/757587093366009856/messages 碎片化阅读没有意义, 捡硬币捡成富翁 kindle不能全文检索 短篇文章的阅读是否有 ...
- 7.1 - CRM系统
一.简介 crm 客户关系管理软件 ( Customer Relationship Management ) ( 详细内容 ) stark组件(仿admin组件)( 详细内容 ) rbac组件(基于角 ...
- mysql实现开窗函数、Mysql实现分析函数
关键字:mysql实现开窗函数.Mysql实现分析函数.利用变量实现窗口函数 注意,变量是从左到右顺序执行的 --测试数据 CREATE TABLE `tem` ( `id` ) NOT NULL A ...
- Unity3d GUI适应分辨率
float sourceWidth = 1024f; float sourceHeight = 768f; float m_fScaleWidth; float m_fScaleHeight; voi ...
- 上传文件小的oK,大一点的传不了,显示 (failed) net::ERR_CONNECTION_RESET
我很确定已经修改了php.ini中的文件上传限制,文件权限可写. 修改php.ini file_uploads = on ;是否允许通过HTTP上传文件的开关.默认为ON即是开 upload_tmp_ ...
- 微博开源框架Motan初体验
前两天,我在开源中国的微信公众号看到新浪微博的轻量Rpc框架--Motan开源了.上网查了下,才得知这个Motan来头不小,支撑着新浪微博的千亿调用,曾经在2014年的春晚中有着千亿次的调用,对抗了春 ...