传送门.

题解:


最主要的问题是如何判断一个数是否合法,这就需要发现性质了。

这个状态划分还是不太容易想到,

每次加的数\(∈[0,k)\),也就是个位一直在变变变,更高的位每次都是加一,这启发我们状态的划分。

这个时候可以利用数位dp的逐位确定思想,在尝试后,发现可以从高位到低位,因为当高位确定后,高位就不会变了,那么高位的最大值也就确定了。

设\(f[i][p][a]\),\(i\)含义如上,\(i+1\)位后的最大值是p,\(2-i\)位是0,当前个位是\(a\),使第\(i\)位加1后个位变成什么?

\(i=2\)时直接暴力处理,\(f[i]\)可以\(O(k)\)由\(f[i-1]\)推出来,复杂度\(O(n*k^3)\)。

有了f方便处理出\(g[i][p][x][a]\),\(i、p、a\)含义如上,x表示第i位要+x,

这里\(x=0\),g的值就是a,然后g自己推自己,复杂度\(O(n*k^3)\)。

接下来的部分就很傻逼了,带根联通块,直接在dfs序上dp,做到个位的时候,再跳跳看能不能跳到那个位去就好了。

Code:


#include<bits/stdc++.h>
#define fo(i, x, y) for(int i = x, B = y; i <= B; i ++)
#define ff(i, x, y) for(int i = x, B = y; i < B; i ++)
#define fd(i, x, y) for(int i = x, B = y; i >= B; i --)
#define ll long long
#define pp printf
#define hh pp("\n")
using namespace std; const int N = 505; int n, k;
int d[N], x, y;
vector<int> e[N];
#define pb push_back
#define si size() int p[N], q[N], np[N], p0, fa[N]; void dg(int x) {
p[x] = ++ p0; np[p0] = x;
ff(j, 0, e[x].si) {
int y = e[x][j];
if(fa[x] != y) fa[y] = x, dg(y);
}
q[x] = p0;
} const int mo = 1e9 + 7; int f[N][10][10], g[N][10][10][10];
ll h[N][N][10][10]; void add(ll &x, ll y) {
(x += y) >= mo ? x -= mo : 0;
} int T; int main() {
// freopen("buried.in", "r", stdin);
// freopen("buried.out", "w", stdout);
k = 10;
for(scanf("%d", &T); T; T --) {
fo(i, 1, n) e[i].clear();
memset(h, 0, sizeof h);
fo(i, 1, n) fa[i] = 0; p0 = 0;
scanf("%d", &n);
fo(i, 1, n - 1) {
scanf("%d %d", &x, &y);
e[x].pb(y); e[y].pb(x);
}
fo(i, 1, n) scanf("%d", &d[i]);
fo(i, 1, n) sort(e[i].begin(), e[i].end());
dg(1);
//Initialization of F
ff(p, 0, k) {
ff(a, 0, k) if(p || a) {
int x = a;
while(x < k) {
x += max(x, p);
}
f[1][p][a] = x % k;
}
}
//dp F
fo(i, 2, n - 1) {
ff(p, 0, k) {
ff(a, 0, k) {
int x = a;
fo(c, 1, k) x = f[i - 1][max(p, c - 1)][x];
f[i][p][a] = x;
}
}
}
//dp g
fo(i, 2, n) {
ff(p, 0, k) {
ff(a, 0, k) {
g[i][p][0][a] = a;
ff(x, 1, k) g[i][p][x][a] = f[i - 1][max(x - 1, p)][g[i][p][x - 1][a]];
}
}
}
//dp h
ll ans = 0;
fo(j, 1, n) h[1][j][0][1] = 1;
fo(i, 1, n) {
int num = d[np[i]];
int ni = q[np[i]] + 1;
fo(j, 2, n) {
ff(p, 0, k) ff(a, 0, k) if(h[i][j][p][a] && (p || a)) {
add(h[i + 1][j - 1][max(p, num)][g[j][p][num][a]], h[i][j][p][a]);
add(h[ni][j][p][a], h[i][j][p][a]);
}
}
ff(p, 0, k) ff(a, 0, k) if(h[i][1][p][a] && (p || a)) {
int x = a;
while(x < num) x += max(x, p);
if(x == num) add(ans, h[i][1][p][a]);
add(h[ni][1][p][a], h[i][1][p][a]);
}
}
pp("%lld\n", ans);
}
}

CodeChef Max-digit Tree(动态规划)的更多相关文章

  1. 【Codeforces715C&716E】Digit Tree 数学 + 点分治

    C. Digit Tree time limit per test:3 seconds memory limit per test:256 megabytes input:standard input ...

  2. Codeforces 716 E Digit Tree

    E. Digit Tree time limit per test 3 seconds memory limit per test 256 megabytes input standard input ...

  3. 【题解】Digit Tree

    [题解]Digit Tree CodeForces - 716E 呵呵以为是数据结构题然后是淀粉质还行... 题目就是给你一颗有边权的树,问你有多少路径,把路径上的数字顺次写出来,是\(m\)的倍数. ...

  4. Codechef Union on Tree

    Codechef Union on Tree https://www.codechef.com/problems/BTREE 简要题意: 给你一棵树,\(Q\)次询问,每次给出一个点集和每个点的\(r ...

  5. 【Codeforces 715C】Digit Tree(点分治)

    Description 程序员 ZS 有一棵树,它可以表示为 \(n\) 个顶点的无向连通图,顶点编号从 \(0\) 到 \(n-1\),它们之间有 \(n-1\) 条边.每条边上都有一个非零的数字. ...

  6. Codechef Observing the Tree

    Home » Practice(Hard) » Observing the Tree   https://www.codechef.com/problems/QUERY Observing the T ...

  7. HDU 1003 Max Sum【动态规划求最大子序列和详解 】

    Max Sum Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Sub ...

  8. CF715C:Digit Tree

    传送门 一句话怎么说来着 算法+高级数据结构=OI 现在我感觉到的是 我会的算法+我会的高级数据结构=WA 这道题提交了三四十次,从刚看题到完全写好花了好几天..,主要死于看错费马小定理的适用条件. ...

  9. CF 716E. Digit Tree [点分治]

    题意:一棵树,边上有一个个位数字,走一条路径会得到一个数字,求有多少路径得到的数字可以整除\(P\) 路径统计一般就是点分治了 \[ a*10^{deep} + b \ \equiv \pmod P\ ...

  10. hdu 1003 Max Sum (动态规划)

    转载于acm之家http://www.acmerblog.com/hdu-1003-Max-Sum-1258.html Max Sum Time Limit: 2000/1000 MS (Java/O ...

随机推荐

  1. react 类样式的一些使用方法

    在 css类不想使用穿透的状态,可以再webpack配置  modules:true,它位于css-loader下,此外调用css的类时可能会自动生成一个hash值,这时候如果想显示本来的名字,可以打 ...

  2. 【C++进阶:移位运算符的用法】

    数据在计算机中以补码存储 移位运算符:<<  左移运算符 >>  右移运算符 一丶 << 左移运算符 移位规则:左边抛弃,右边补零 int num=10; num& ...

  3. ConcurrentLinkedQueue 源码分析

    ConcurrentLinkedQueue ConcurrentLinkedQueue 能解决什么问题?什么时候使用 ConcurrentLinkedQueue? 1)ConcurrentLinked ...

  4. 阶段1 语言基础+高级_1-3-Java语言高级_04-集合_05 List集合_3_LinkedList集合

    链表实现,查询慢,增删快 addFirst 集合的开头添加 push等效于addFirst addLast getFirst和getLast clear 清空后 又输出 isEmpty 不为空才去获取 ...

  5. django-xadmin常用内容记录

    自定义菜单名称: 1 修改app下的 apps.py文件 添加 class OperationConfig(AppConfig): name = 'operation' verbose_name = ...

  6. Linux-部署ftp

    通过外部window ftp 客户端 访问linux 有两种方法 方法一:Linux系统未安装vsftp 服务 这个是本人使用的ftp客户端的版本号 启动ftp客户端,填写ip ,账号,密码 问题:当 ...

  7. vue集成汉字转拼音或提取首字母

    需求:             有时我们为了节省用户的维护量,需要根据中文生成出相应的拼音和缩写 解决:            此方法是利用汉字和Unicode编码对应找到相应字母 一.编写汉字和编码 ...

  8. 创建客户端项目并读取服务化的配置中心(Consul + Spring Cloud Config)

    创建客户端项目并读取服务化的配置中心 将配置中心注册到服务中心(Consul) POM文件添加依赖: <dependency> <groupId>org.springframe ...

  9. [2019杭电多校第七场][hdu6646]A + B = C(hash)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6646 题意为求a*10x+b*10y=c*10z满足公式的任意一组解x,y,z. 因为c有可能会由a+ ...

  10. Qt5 对xml文件常用的操作(读写,增删改查)

    转自:https://blog.csdn.net/hpu11/article/details/80227093 项目配置 pro文件里面添加QT+=xml include <QtXml>, ...