A. Drazil and Factorial
题意:
F(x)定义为x各个位阶乘的乘积,给定一个数x要求求出F(x)=F(y)的最大的数y, 其中y的各个位不能包含0和1.
思路:
对于x的各个位, 先把阶乘分解质因子, 然后有多少个质因子7就弄成多少个7个阶乘, 5,3和2也类似。
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <stack>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <queue>
#include <ctime>
#include <cstdlib>
using namespace std;
#define mxn 100020
#define mxe 200020
#define N 10020
#define M 2000020
#define mod 2147483647
#define LL long long
#define inf 0x3f3f3f3f
#define vi vector<int>
#define PB push_back
#define MP make_pair
#define pii pair<int, int>
#pragma comment(linker,"/STACK:1024000000,1024000000")
#define ls (i << 1)
#define rs (ls | 1)
#define md ((ll + rr) >> 1)
LL n;
int cnt[10];
void f(int k) {
if(k == 0 || k == 1) return;
for(int i = 2; i * i <= k; ++i) {
if(k % i) continue;
while(k % i == 0) {
cnt[i]++;
k /= i;
}
}
if(k > 1)
cnt[k]++;
}
int ans[1000], len;
char s[222];
int main() {
// freopen("tt.txt", "r", stdin);
while(cin >> n) {
scanf("%s", s);
memset(cnt, 0, sizeof cnt);
for(int i = 0; s[i]; ++i) {
int t = s[i] - '0';
for(int i = 2; i <= t; ++i)
f(i);
}
len = 0;
while(cnt[7] > 0) {
cnt[2] -= 4, cnt[3] -=2, cnt[5]--, cnt[7]--;
ans[len++] = 7;
}
while(cnt[5] > 0) {
cnt[2] -= 3, cnt[3]--, cnt[5]--;
ans[len++] = 5;
}
while(cnt[3] > 0) {
cnt[2] --, cnt[3]--;
ans[len++] = 3;
}
while(cnt[2] > 0) {
--cnt[2];
ans[len++] = 2;
}
for(int i = 0; i < len; ++i)
printf("%d", ans[i]);
puts("");
}
return 0;
}
B.Drazil and Tiles
题意:
给定一个n*m的矩阵, 要求在矩阵的空位上放置1*2的牌子, 要把所有空位放满, 如果存在多种方案或者没有方案就输出Not unique, 否则输出方案。
思路:
把矩阵的空位看成点, 上下左右相邻的空位连一条边, 对于那些度数为1的点是只有一种方案的, 如果有度数为0的空位,就输出not unique, 否则一直找度数为1的点放牌子, 最后如果剩下的空位度数都大于1的话肯定存在不同方案, 因为剩下的空位肯定会构成环, 所以输出not unique就行。
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <stack>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <queue>
#include <ctime>
#include <cstdlib>
using namespace std;
#define mxn 2020
#define mxe 200020
#define N 10020
#define M 2000020
#define mod 2147483647
#define LL long long
#define inf 0x3f3f3f3f
#define vi vector<int>
#define PB push_back
#define MP make_pair
#define pii pair<int, int>
#pragma comment(linker,"/STACK:1024000000,1024000000")
#define ls (i << 1)
#define rs (ls | 1)
#define md ((ll + rr) >> 1)
int n, m;
char s[mxn][mxn];
int dir[4][2] = {0, 1, 1, 0, 0, -1, -1, 0};
bool in(int x, int y) {
return x >= 0 && y >= 0 && x < n && y < m && s[x][y] != '*';
}
bool vis[mxn][mxn];
int ind[mxn][mxn];
void col(int x, int y, int dx, int dy) {
if(x > dx) swap(x, dx);
if(y > dy) swap(y, dy);
if(x == dx)
s[x][y] = '<', s[dx][dy] = '>';
else
s[x][y] = '^', s[dx][dy] = 'v';
}
int main() {
// freopen("tt.txt", "r", stdin);
while(scanf("%d%d", &n, &m) != EOF) {
for(int i = 0; i < n; ++i) scanf("%s", s[i]);
memset(vis, 0, sizeof vis);
memset(ind, 0, sizeof ind);
queue<int> q;
for(int i = 0; i < n; ++i) {
for(int j = 0; j < m; ++j) {
if(s[i][j] == '*') continue;
for(int k = 0; k < 4; ++k) {
int dx = i + dir[k][0];
int dy = j + dir[k][1];
if(!in(dx, dy)) continue;
ind[i][j]++;
}
if(ind[i][j] <= 1)
q.push(i), q.push(j);
}
}
bool flag = 1;
while(!q.empty()) {
int x = q.front(); q.pop();
int y = q.front(); q.pop();
// printf("%d %d\n", x, y);
if(vis[x][y]) continue;
if(ind[x][y] != 1) {
flag = 0;
break;
}
for(int i = 0; i < 4; ++i) {
int dx = x + dir[i][0];
int dy = y + dir[i][1];
if(!in(dx, dy) || vis[dx][dy]) continue;
vis[x][y] = vis[dx][dy] = 1;
col(x, y, dx, dy);
for(int j = 0; j < 4; ++j) {
int xx = dx + dir[j][0];
int yy = dy + dir[j][1];
if(!in(xx, yy) || vis[xx][yy]) continue;
ind[xx][yy]--;
if(ind[xx][yy] <= 1)
q.push(xx), q.push(yy);
}
}
}
for(int i = 0; i < n; ++i) {
for(int j = 0; j < m; ++j)
if(in(i, j) && vis[i][j] == 0) {
flag = 0;
break;
}
}
if(!flag)
puts("Not unique");
else
for(int i = 0; i < n; ++i)
printf("%s\n", s[i]);
}
return 0;
}
C.Drazil and Park
题意:
有n棵树都成一个环, 第i个树和第i+1个树之间的距离是d[i],每个树有一定的高度h[i], m次询问,每次询问一段区间内任意两个端点u,v, 2 * h[u] + 2 * h[v] + d[u] + d[u+1] + ... +d[v]的最大值。
思路:
先把环转化成线性的问题,就是1,2,3...n, 1, 2, 3, ....n。
每次询问如果a <= b, 就转化成询问区间[b+1, a + n - 1],否则询问区间[b+1,a-1].
对于每次询问[l,r]:
2 * h[u] + 2 * h[v] + dist(u, v) = 2 * h[v] + d[1] + d[2] + ... + d[v-1] + 2 * h[u] - (d[1] + d[2] + ... + d[u-1]).
设A[v] = 2 * h[v] + d[1] + d[2] + ... + d[v-1], B[u] = 2 * h[u] - (d[1] + d[2] + ... + d[u-1]).
由于u < v 时, A[u] + B[v] < A[v] + B[u], 所以就转化成求[l,r]区间内A + B的最大值,用线段树可以logn求出。
但是由于u不能等于v, 所以线段树需要保存3个变量, 就是区间的A的最大值, 区间内B的最大值还有A+B并且u!=v的最大值。
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <stack>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <queue>
#include <ctime>
#include <cstdlib>
using namespace std;
#define mxn 200020
#define mxe 200020
#define M 2000020
#define mod 2147483647
#define LL long long
#define inf 0x3f3f3f3f
#define vi vector<int>
#define PB push_back
#define MP make_pair
#define pii pair<int, int>
#pragma comment(linker,"/STACK:1024000000,1024000000")
#define ls (i << 1)
#define rs (ls | 1)
#define md ((ll + rr) >> 1)
int n, m;
LL d[mxn], h[mxn];
LL s[mxn<<2], a[mxn<<2], b[mxn<<2];
void pushUp(int ll, int rr, int i) {
s[i] = max(s[ls], s[rs]);
s[i] = max(s[i], a[ls] + b[rs]);
s[i] = max(s[i], a[rs] + b[ls]);
a[i] = max(a[ls], a[rs]);
b[i] = max(b[ls], b[rs]);
}
void build(int ll, int rr, int i) {
if(ll == rr) {
s[i] = 0;
a[i] = 2 * h[ll] + d[ll-1];
b[i] = 2 * h[ll] - d[ll-1];
return;
}
build(ll, md, ls), build(md + 1, rr, rs);
pushUp(ll, rr, i);
}
LL query(LL &x, LL &y, int l, int r, int ll, int rr, int i) {
if(ll == l && rr == r) {
x = a[i], y = b[i];
return s[i];
}
if(r <= md) {
return query(x, y, l, r, ll, md, ls);
}
if(l > md) return query(x, y, l, r, md + 1, rr, rs);
LL dx, dy, ix, iy;
LL tmp1 = query(dx, dy, l, md, ll, md, ls);
LL tmp2 = query(ix, iy, md + 1, r, md + 1, rr, rs);
LL ret = max(tmp1, tmp2);
ret = max(ret, dx + iy);
ret = max(ret, ix + dy);
x = max(dx, ix);
y = max(dy, iy);
return ret;
}
int main() {
while(scanf("%d%d", &n, &m) != EOF) {
for(int i = 1; i <= n; ++i) scanf("%I64d", &d[i]);
for(int i = 1; i <= n; ++i) scanf("%I64d", &h[i]);
for(int i = 1; i <= n; ++i) {
d[i+n] = d[i];
h[i+n] = h[i];
}
for(int i = 1; i <= 2 * n; ++i)
d[i] = d[i-1] + d[i];
int N = 2 * n;
build(1, N, 1);
while(m--) {
int l, r;
scanf("%d%d", &l, &r);
if(l > r) {
swap(l, r);
++l, --r;
}
else {
swap(l, r);
++l, r = r + n - 1;
}
LL x, y;
printf("%I64d\n", query(x, y, l, r, 1, N, 1));
}
}
return 0;
}
D.Drazil and Morning Exercise
题意:
给定一颗带权树,有q个询问, 每次需要在树上求出一个最大的联通块,并且联通块每个节点距离所有叶子的距离最大值的最大值和最小值的差不能超过L。
思路:
先通过两次dfs求出每个节点到叶子节点的距离的最大值far[i]。
然后取最小的far的节点做为根节点, 把树转化成有根树。
可以证明,u节点的子树中任意节点的far值都大于或者等于far[u]。
假设根节点是root, 距离它最远的叶子节点是v, v在u的儿子w的子树中,对于u的除了w的其他子树的节点, 距离它们最远的叶子节点肯定是v, 所以他们的far值是递增的。
而对于w子树, 距离它们最远的叶子节点一定不是v,否则他们的far值会比u还小。
求出far后, 每个节点如果作为far最小的节点, 那个所选的节点一定是在他的子树中, 假设选取的是u子树, 那个要尽量贪心地往下选取节点, 只要far[v] + far[u] <= L,所以对于所有节点u在祖先中都有一定的作用范围, 也就是有一段连续的祖先都能够选它, 所以可以类似lca那样求出第一个不能选它的祖先, 并且在这个祖先标记一下。然后做dfs的时候可以把那个第一个不能够选它的祖先减去1。
#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <stack>
#include <cmath>
#include <set>
#include <map>
#include <vector>
#include <queue>
#include <ctime>
#include <cstdlib>
using namespace std;
#define mxn 100020
#define mxe 200020
#define M 2000020
#define mod 2147483647
#define LL long long
#define inf 0x3f3f3f3f
#define vi vector<int>
#define PB push_back
#define MP make_pair
#define pii pair<int, int>
#pragma comment(linker,"/STACK:1024000000,1024000000")
#define ls (i << 1)
#define rs (ls | 1)
#define md ((ll + rr) >> 1)
int n, m;
int fst[mxn], nxt[mxe], cost[mxe], to[mxe], e;
void init() {
memset(fst, -1, sizeof fst);
e = 0;
}
void add(int u, int v, int c) {
to[e] = v, nxt[e] = fst[u], cost[e] = c, fst[u] = e++;
}
LL mx[mxn][3];
LL far[mxn];
int fa[22][mxn];
bool cmp(LL x, LL y) {
return x > y;
}
void f(int u, LL x) {
mx[u][2] = x;
sort(mx[u], mx[u] + 3, cmp);
}
void dfs1(int u, int p, LL d) {
mx[u][0] = mx[u][1] = 0;
for(int i = fst[u]; ~i; i = nxt[i]) {
int v = to[i];
if(v == p) continue;
dfs1(v, u, d + cost[i]);
f(u, cost[i] + mx[v][0]);
}
// printf("%d %I64d %I64d\n", u, mx[u][0], mx[u][1]);
}
void dfs2(int u, int p, LL d) {
far[u] = max(d, mx[u][0]);
for(int i = fst[u]; ~i; i = nxt[i]) {
int v = to[i], c = cost[i];
if(v == p) continue;
if(mx[u][0] == c + mx[v][0])
dfs2(v, u, max(d, mx[u][1]) + c);
else
dfs2(v, u, max(d, mx[u][0]) + c);
}
}
LL L;
int ans, cnt;
int cc[mxn];
void dfs3(int u, int p) {
fa[0][u] = p;
for(int i = fst[u]; ~i; i = nxt[i]) {
int v = to[i];
if(v != p)
dfs3(v, u);
}
}
int f(int u) {
int tmp = u;
for(int k = 21; k >= 0; --k) {
if(fa[k][tmp] != -1) {
int v = fa[k][tmp];
if(far[u] - far[v] <= L)
tmp = v;
}
}
// printf("in f() %d %d\n", tmp, u);
return fa[0][tmp];
}
void dfs4(int u, int p) {
// printf("pre %d %d\n", u, cnt);
int tmp = cnt++;
for(int i = fst[u]; ~i; i = nxt[i]) {
int v = to[i];
if(v == p) continue;
dfs4(v, u);
}
cnt -= cc[u];
// printf("cc[%d]: %d\n", u, cc[u]);
int t = f(u);
if(t != -1)
cc[t]++;
// printf("nxt %d %d %d\n", u, cnt, tmp);
ans = max(ans, cnt - tmp);
}
int main() {
// freopen("tt.txt", "r", stdin);
while(scanf("%d", &n) != EOF) {
init();
for(int i = 1; i < n; ++i) {
int u, v, c;
scanf("%d%d%d", &u, &v, &c);
add(u, v, c);
add(v, u, c);
}
dfs1(1, -1, 0);
dfs2(1, -1, 0);
int root = 1;
for(int i = 2; i <= n; ++i)
if(far[root] > far[i])
root = i;
dfs3(root, -1);
// for(int i = 1; i <= n; ++i)
// printf("%I64d ", far[i]);
// puts("");
for(int k = 0; k < 21; ++k)
for(int i = 1; i <= n; ++i) {
if(fa[k][i] != -1)
fa[k+1][i] = fa[k][fa[k][i]];
else
fa[k+1][i] = -1;
}
scanf("%d", &m);
while(m--) {
scanf("%I64d", &L);
cnt = ans = 0;
memset(cc, 0, sizeof cc);
dfs4(root, -1);
printf("%d\n", ans);
}
}
return 0;
}