# |
Submission time |
Handle |
Problem |
Language |
Result |
Execution time |
Memory |
607760 |
2022-07-27T03:13:58 Z |
tmt514 |
Wells (CEOI21_wells) |
C++14 |
|
24 ms |
47304 KB |
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1500001;
const long long MOD = 1e9+7;
vector<int> a[MAXN];
int e[MAXN][2];
int dp[MAXN][2];
int isfree[MAXN];
int depth[MAXN];
int N, K;
// compute longest path on edge i direction d
int go(int i, int d) {
if(dp[i][d] != -1) return dp[i][d];
// fprintf(stderr, "i=%d, d=%d\n", i, d);
int x = e[i][d], y = e[i][d^1];
dp[i][d] = 2;
for(int j : a[y]) {
if (j != i) {
int d2 = (y == e[j][1]);
dp[i][d] = max(dp[i][d], go(j, d2) + 1);
}
}
return dp[i][d];
}
// test
bool test(int x, int p, int level) {
if (isfree[x]) return true;
bool success = true;
if (level == 0) {
for(int i : a[x]) {
int d = (x == e[i][1]);
int y = e[i][d^1];
if (y == p) continue;
success &= test(y, x, K-1);
if (!success) return false;
}
} else {
int deepmost = 0, deeptwo = 0;
for(int i : a[x]) {
int d = (x == e[i][1]);
int y = e[i][d^1];
if (y == p) continue;
success &= test(y, x, level - 1);
if (!success) return false;
if (dp[i][d] > deepmost) {
deeptwo = deepmost;
deepmost = dp[i][d] - 1;
} else if (dp[i][d] > deeptwo) {
deeptwo = dp[i][d] - 1;
}
}
if (deepmost >= level && deeptwo >= level && level * 2 + 1 <= K) {
return false;
}
if (min(level-1, deepmost) + min(level-1, deeptwo) + 1 >= K) {
return false;
}
}
return success;
}
int main() {
scanf("%d%d", &N, &K);
for(int i=1;i<N;i++) {
int x, y;
scanf("%d%d", &x, &y);
a[x].push_back(i);
a[y].push_back(i);
e[i][0] = x;
e[i][1] = y;
}
memset(dp, -1, sizeof(dp));
for(int i=1;i<N;i++) {
for(int d=0;d<2;d++) go(i, d);
}
// check if each vertex is a "free vertex"
int freecount=0;
for(int x=1;x<=N;x++) {
vector<int> lens;
for(int i : a[x]) {
int d = (x == e[i][1]);
lens.push_back(dp[i][d]-1);
}
sort(lens.begin(), lens.end(), greater<>());
int big = (lens.size()>=1? lens[0]:0) + (lens.size()>=2? lens[1]:0);
if(big+1<K) {
isfree[x]=1;
++freecount;
}
}
// trivial case: no vertex belongs to a path of length K.
if (freecount == N) {
long long ans = 1;
for(int i=0;i<N;i++) ans=ans*2%MOD;
printf("YES\n%lld\n", ans);
return 0;
}
// find one non-free vertex, and do DFS.
// find a path of length k.
auto get_one_path_start = [&]() {
for(int x=1;x<=N;x++) {
for(int i: a[x]) {
int d = (x==e[i][1]);
if(dp[i][d] >= K)
return x;
}
}
return -1;
};
int x = get_one_path_start();
vector<int> path;
while(path.size() < K-1) {
for(int i:a[x]) {
int d = (x==e[i][1]);
if((path.size()==0 || e[i][d^1]!=path.back()) && dp[i][d] >= K - (int)path.size()) {
path.push_back(x);
x = e[i][d^1];
break;
}
}
}
path.push_back(x);
// fprintf(stderr, "path: ");
// for(int i:path) fprintf(stderr, "%d, ", i);
// fprintf(stderr, "\n");
int possible = 0;
for(int y : path) {
if (test(y, -1, 0)) ++possible;
}
if (!possible) {
printf("NO\n0\n");
} else {
for(int i=0;i<freecount;i++) {
possible = possible*2%MOD;
}
printf("YES\n%d\n", possible);
}
return 0;
}
Compilation message
wells.cpp: In function 'int go(int, int)':
wells.cpp:18:7: warning: unused variable 'x' [-Wunused-variable]
18 | int x = e[i][d], y = e[i][d^1];
| ^
wells.cpp: In function 'int main()':
wells.cpp:123:21: warning: comparison of integer expressions of different signedness: 'std::vector<int>::size_type' {aka 'long unsigned int'} and 'int' [-Wsign-compare]
123 | while(path.size() < K-1) {
| ~~~~~~~~~~~~^~~~~
wells.cpp:71:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
71 | scanf("%d%d", &N, &K);
| ~~~~~^~~~~~~~~~~~~~~~
wells.cpp:74:10: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
74 | scanf("%d%d", &x, &y);
| ~~~~~^~~~~~~~~~~~~~~~
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
21 ms |
47188 KB |
Output is correct |
2 |
Partially correct |
21 ms |
47304 KB |
Output is partially correct |
3 |
Partially correct |
21 ms |
47200 KB |
Output is partially correct |
4 |
Correct |
23 ms |
47188 KB |
Output is correct |
5 |
Partially correct |
24 ms |
47284 KB |
Output is partially correct |
6 |
Partially correct |
21 ms |
47244 KB |
Output is partially correct |
7 |
Correct |
21 ms |
47284 KB |
Output is correct |
8 |
Correct |
21 ms |
47188 KB |
Output is correct |
9 |
Correct |
24 ms |
47164 KB |
Output is correct |
10 |
Correct |
22 ms |
47188 KB |
Output is correct |
11 |
Incorrect |
22 ms |
47188 KB |
Output isn't correct |
12 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
21 ms |
47188 KB |
Output is correct |
2 |
Partially correct |
21 ms |
47304 KB |
Output is partially correct |
3 |
Partially correct |
21 ms |
47200 KB |
Output is partially correct |
4 |
Correct |
23 ms |
47188 KB |
Output is correct |
5 |
Partially correct |
24 ms |
47284 KB |
Output is partially correct |
6 |
Partially correct |
21 ms |
47244 KB |
Output is partially correct |
7 |
Correct |
21 ms |
47284 KB |
Output is correct |
8 |
Correct |
21 ms |
47188 KB |
Output is correct |
9 |
Correct |
24 ms |
47164 KB |
Output is correct |
10 |
Correct |
22 ms |
47188 KB |
Output is correct |
11 |
Incorrect |
22 ms |
47188 KB |
Output isn't correct |
12 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
21 ms |
47188 KB |
Output is correct |
2 |
Partially correct |
21 ms |
47304 KB |
Output is partially correct |
3 |
Partially correct |
21 ms |
47200 KB |
Output is partially correct |
4 |
Correct |
23 ms |
47188 KB |
Output is correct |
5 |
Partially correct |
24 ms |
47284 KB |
Output is partially correct |
6 |
Partially correct |
21 ms |
47244 KB |
Output is partially correct |
7 |
Correct |
21 ms |
47284 KB |
Output is correct |
8 |
Correct |
21 ms |
47188 KB |
Output is correct |
9 |
Correct |
24 ms |
47164 KB |
Output is correct |
10 |
Correct |
22 ms |
47188 KB |
Output is correct |
11 |
Incorrect |
22 ms |
47188 KB |
Output isn't correct |
12 |
Halted |
0 ms |
0 KB |
- |
# |
Verdict |
Execution time |
Memory |
Grader output |
1 |
Correct |
21 ms |
47188 KB |
Output is correct |
2 |
Partially correct |
21 ms |
47304 KB |
Output is partially correct |
3 |
Partially correct |
21 ms |
47200 KB |
Output is partially correct |
4 |
Correct |
23 ms |
47188 KB |
Output is correct |
5 |
Partially correct |
24 ms |
47284 KB |
Output is partially correct |
6 |
Partially correct |
21 ms |
47244 KB |
Output is partially correct |
7 |
Correct |
21 ms |
47284 KB |
Output is correct |
8 |
Correct |
21 ms |
47188 KB |
Output is correct |
9 |
Correct |
24 ms |
47164 KB |
Output is correct |
10 |
Correct |
22 ms |
47188 KB |
Output is correct |
11 |
Incorrect |
22 ms |
47188 KB |
Output isn't correct |
12 |
Halted |
0 ms |
0 KB |
- |