链接:2018 Multi-University Training Contest 3
过题数:5
排名:146/815
A. Ascending Rating
题意
给定一个长度为 n n 的序列 ,对每一个长度为 m m 的区间 ,都有以下操作:
- 最开始有一个 maxratingi m a x r a t i n g i 与 counti c o u n t i 值,初始值分别为 −1 − 1 和 0 0 ;
- 从 到 r r ,每扫描一个数字,如果当前数字大于 ,就将 maxratingi m a x r a t i n g i 更新为当前数字,并将 counti+1 c o u n t i + 1 。
每一个连续的区间 maxratingi m a x r a t i n g i 和 counti c o u n t i 的初始值都为 −1 − 1 和 0 0 ,求每个区间最终的 值和 counti c o u n t i 值。
输入
第一行为一个整数 T (1≤T≤2000) T ( 1 ≤ T ≤ 2000 ) ,接下来有 T T 组数据,每组数据第一行为 个整数 n,m,k,p,q,r,MOD(1≤m,k≤n≤107,5≤p,q,r,MOD≤109) n , m , k , p , q , r , M O D ( 1 ≤ m , k ≤ n ≤ 10 7 , 5 ≤ p , q , r , M O D ≤ 10 9 ) ,第二行为 k k 个整数 ,表示序列的前 k k 项,第 到第 n n 项由以下公式生成:
数据保证 ∑n≤7×107 ∑ n ≤ 7 × 10 7 且 ∑k≤2×106 ∑ k ≤ 2 × 10 6 。
输出
每组数据输出两个数字 A A 和 ,其中:
A=∑i=1n−m+1maxratingi⨁iB=∑i=1n−m+1counti⨁i A = ∑ i = 1 n − m + 1 m a x r a t i n g i ⨁ i B = ∑ i = 1 n − m + 1 c o u n t i ⨁ i
样例
输入 |
---|
1 10 6 10 5 5 5 5 3 2 2 1 5 7 6 8 2 9 |
输出 |
46 11 |
题解
从后往前维护一个单调队列,维护到第 i i 个区间时队列的大小即 ,队列中的最大值就是 maxratingi m a x r a t i n g i ,扫一遍只需要 O(n) O ( n ) 。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
#include <unordered_set>
#include <unordered_map>
using namespace std;
#define LL long long
const int maxn = 10000000 + 100;
int T, n, m, k, p, q, r, MOD, head, tail;
LL A, B;
int num[maxn], que[maxn];
int add(int a, int b) {
int ret = a + b;
if(ret >= MOD) {
return ret - MOD;
}
return ret;
}
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(true);
scanf("%d", &T);
while(T--) {
head = tail = 0;
scanf("%d%d%d%d%d%d%d", &n, &m, &k, &p, &q, &r, &MOD);
for(int i = 1; i <= k; ++i) {
scanf("%d", &num[i]);
}
for(int i = k + 1; i <= n; ++i) {
num[i] = add((LL)p * num[i - 1] % MOD, (LL)q * i % MOD);
num[i] = add(num[i], r);
}
for(int i = 0; i < m; ++i) {
int Index = n - i;
while(head != tail && num[Index] >= num[que[head]]) {
--head;
}
que[++head] = Index;
}
A = num[que[tail + 1]] ^ (n - m + 1);
B = (head - tail) ^ (n - m + 1);
for(int i = n - m; i >= 1; --i) {
while(head != tail && num[i] >= num[que[head]]) {
--head;
}
que[++head] = i;
if(que[tail + 1] >= i + m) {
++tail;
}
A += num[que[tail + 1]] ^ i;
B += (head - tail) ^ i;
}
printf("%I64d %I64d\n", A, B);
}
return 0;
}
C. Dynamic Graph Matching
题意
定义一种图的匹配:在图上选择任意条边,在被选择的边中任意两条边之间没有公共点,则这些边就是图上的一种匹配。
在一个 n n 个节点的图上,最初这张图上的边数为 ,有 m m 次操作,每次操作为以下两种中的一种:
- :在节点 u u 和节点 之间添加一条边;
- − u v − u v :删掉一条在节点 u u 和 之间的边,在删边之前保证节点 u,v u , v 之间至少存在一条边。
在每次操作之后,对于每一个 k (k∈[1,n2]) k ( k ∈ [ 1 , n 2 ] ) ,输出图上边数为 k k 的匹配的数量,两个节点之间的多条边视为不同的边。
输入
第一行为一个整数 ,接下去有 T T 组数据,每组数据第一行为两个整数 ,接下去 m m 行每行由一个字符 和两个数字 u,v (1≤u,v≤n) u , v ( 1 ≤ u , v ≤ n ) 组成,表示一次操作。
输出
对于每一次操作,输出 n2 n 2 个整数,分别表示当 k∈[1,n2] k ∈ [ 1 , n 2 ] 时,图上边数为 k k 的匹配数量对 取模的结果。
样例
输入 |
---|
1 4 8 + 1 2 + 3 4 + 1 3 + 2 4 - 1 2 - 3 4 + 1 2 + 3 4 |
输出 |
1 0 2 1 3 1 4 2 3 1 2 1 3 1 4 2 |
题解
我们用 2n 2 n 个整数对应图上点的选择方案,整数的第 i i 个二进制位为 表示不选择第 i i 个点,否则表示选择第 个点,则一个数字表示一种选择点的状态。用 dp[i] d p [ i ] 记录某种匹配的边覆盖的点的状态为 i i 的情况下,匹配的数量,则 中 1 1 的个数必须为偶数,否则匹配数量就为 。
每当节点 u u 和 之间的边的数量 +1 + 1 或 −1 − 1 ,就会影响所有选点状态同时包含 u,v u , v 两点的 dp d p 值,如果 i i 的第 和 v v 位都为 ,则 dp[i]+=d×dp[i−(1<<v)−(1<<u)] d p [ i ] + = d × d p [ i − ( 1 << v ) − ( 1 << u ) ] ,当 ch c h 为+
时 d=1 d = 1 ,否则 d=−1 d = − 1 。
初始状态 dp[0]=1 d p [ 0 ] = 1 ,用 ans[i] a n s [ i ] 记录在匹配边的数量为 i i 时的总方案数,设 为数字 i i 的二进制位中 的数量,则 ans[j]=∑2ni=0dp[i] a n s [ j ] = ∑ i = 0 2 n d p [ i ] ,其中 j=biti2 j = b i t i 2 。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
#include <unordered_set>
#include <unordered_map>
using namespace std;
#define LL long long
const int maxn = 1000 + 100;
const int MOD = 1000000000 + 7;
int T, n, m, u, v, d;
char ch[2];
int dp[maxn], ans[maxn], bit[maxn];
int add(int a, int b) {
int ret = a + b;
if(ret < 0) {
return ret + MOD;
}
if(ret >= MOD) {
return ret - MOD;
}
return ret;
}
void solve(int u, int v, int d) {
int tmp = (1 << u) | (1 << v);
for(int i = 1; i < (1 << n); ++i) {
if((bit[i] & 1) == 1) {
continue;
}
if((i & tmp) == tmp) {
dp[i] = add(dp[i], d * dp[i ^ tmp]);
ans[bit[i] >> 1] = add(ans[bit[i] >> 1], d * dp[i ^ tmp]);
}
}
}
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(true);
for(int i = 0; i < maxn; ++i) {
bit[i] = 0;
for(int j = 0; j < 10; ++j) {
if(((i >> j) & 1) == 1) {
++bit[i];
}
}
}
scanf("%d", &T);
while(T--) {
scanf("%d%d", &n, &m);
memset(dp, 0, sizeof(int) * (1 << n));
memset(ans + 1, 0, sizeof(int) * (n >> 1));
dp[0] = 1;
for(int i = 0; i < m; ++i) {
scanf("%s%d%d", ch, &u, &v);
d = ch[0] == '+'? 1: -1;
solve(u - 1, v - 1, d);
for(int j = 1; j <= (n >> 1); ++j) {
if(j != 1) {
printf(" ");
}
printf("%d", ans[j]);
}
printf("\n");
}
}
return 0;
}
D. Euler Function
题意
输出第 k k 小的 为合数的 n n 的值,其中 为欧拉函数( ϕ(x) ϕ ( x ) 表示 i∈[1,x] i ∈ [ 1 , x ] 中 gcd(i,x)=1 gcd ( i , x ) = 1 的数字个数)。
输入
第一行为一个整数 T (1≤T≤105) T ( 1 ≤ T ≤ 10 5 ) ,接下去 T T 行每行一个整数 。
输出
对于每一个 k k ,输出第 小的合法的数字。
样例
输入 |
---|
2 1 2 |
输出 |
5 7 |
题解
通过打表找规律可以发现,除了 k=1 k = 1 时答案为 5 5 外,其他答案都为 。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cfloat>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <algorithm>
using namespace std;
#define LL long long
int main() {
#ifdef Dmaxiya
freopen("test.txt", "r", stdin);
#endif // Dmaxiya
ios::sync_with_stdio(false);
int T, k;
scanf("%d", &T);
while(T--) {
scanf("%d", &k);
if(k == 1) {
printf("5\n");
} else {
printf("%d\n", k + 5);
}
}
return 0;
}
F. Grab The Tree
题意
小 Q Q 和小 玩一个游戏,在一棵 n n 个节点的树上,每个节点都有一个权值 ,由小 Q Q 先开始选点,小 选择的点中,任意两个节点不能同时为某条边的两个端点,之后剩下的点全是小 T T 的,小 的最终分数就是他所选择的所有点的权值异或值,小 T T 的分数也是他拥有的点的权值异或值,问在小 采取最优策略的情况下,谁会获胜,还是平局。
输入
第一行包含一个整数 T (1≤T≤20) T ( 1 ≤ T ≤ 20 ) ,接下去有 T T 组数据,每组数据第一行为一个整数 ,第二行为 n n 个整数 ,接下去 n−1 n − 1 行每行两个整数 u,v (1≤u,v≤n) u , v ( 1 ≤ u , v ≤ n ) ,表示节点 u,v u , v 之间有一条边。
输出
如果小 Q Q 会获胜,则输出 ,如果小 T T 会获胜,则输出 ,否则输出 D D 。
样例
输入 |
---|
1 3 2 2 2 1 2 1 3 |
输出 |
Q |
题解
如果树上所有点的异或值为 ,无论小 Q Q 如何选点,小 的分数都会和小 Q Q 的分数相等,如果树上所有点的异或值非 ,则小 Q Q 只需要选择最大值,就可以让树上所有其他点异或值的最高位小于小 选择的点的权值,小 Q Q 必胜。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
#include <unordered_set>
#include <unordered_map>
using namespace std;
#define LL long long
int T, n, num, Xor, u, v;
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(true);
scanf("%d", &T);
while(T--) {
Xor = 0;
scanf("%d", &n);
for(int i = 0; i < n; ++i) {
scanf("%d", &num);
Xor ^= num;
}
for(int i = 1; i < n; ++i) {
scanf("%d%d", &u, &v);
}
if(Xor == 0) {
printf("D\n");
} else {
printf("Q\n");
}
}
return 0;
}
G. Interstellar Travel
题意
在一个二维平面上总共有 个点,坐标分别为 (xi,yi) ( x i , y i ) ,小 Q Q 要从坐标为 的点到达坐标 (xn,yn) ( x n , y n ) 的点,他每次只能从横坐标小的点到达横坐标大的点,即 xstart<xto x s t a r t < x t o ,如果他要从坐标为 (xi,yi) ( x i , y i ) 的点到达 (xj,yj) ( x j , y j ) 的点,则需要的代价为 xi×yj−xj×yi x i × y j − x j × y i ,问要完成目标的最小代价。
输入
第一行为一个整数 T (1≤T≤10) T ( 1 ≤ T ≤ 10 ) ,接下去有 T T 组数据,每组数据第一行为一个整数 ,接下去 n n 行每行两个整数 ,数据保证 y1=yn=0,0=x1<x2,x3,⋯,xn−1<xn y 1 = y n = 0 , 0 = x 1 < x 2 , x 3 , ⋯ , x n − 1 < x n 。
输出
每组数据输出一个序列 p1,p2,⋯,pm (1≤pi≤n) p 1 , p 2 , ⋯ , p m ( 1 ≤ p i ≤ n ) ,表示小 Q Q 选择的路径为 , pi p i 表示路径上第 i i 个点的编号为 ,点的编号按输入顺序从 1 1 到 确定,其中 p1=1,pm=n p 1 = 1 , p m = n ,如果有多组解,输出字典序最小的一组。
样例
输入 |
---|
1 3 0 0 3 0 4 0 |
输出 |
1 2 3 |
题解
起点为 (0,0) ( 0 , 0 ) ,且经过两个点的代价为两个点对应的向量叉积,即这两点与原点构成的三角形面积,若从向量 i i 到向量 是以顺时针旋转,则这两个向量的叉积就为负数,为了使代价最小,也就是使负数的绝对值——三角形的面积最大,因此这条路径就是从原点到终点的一个凸壳。
首先对所有点去重,然后构造凸壳,构造过程中保留凸壳上的共线点,所有拐点必然是要作为答案的,但是共线点上可能存在比下一个拐点编号更小的点,为使答案的字典序最小,必须选择这些共线点上编号小于下一个拐点的一系列点,且这一系列点的字典序也要最小,在这些共线点中可以用单调栈来维护需要选择的字典序最小的点的编号。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;
#define LL long long
const int maxn = 200000 + 100;
struct Point {
int x, y;
int Index;
Point() {}
Point(int xx, int yy) {
x = xx;
y = yy;
}
};
bool operator<(const Point &a, const Point &b) {
return a.x == b.x? a.y < b.y: a.x < b.x;
}
LL operator^(const Point &a, const Point &b) {
return (LL)a.x * b.y - (LL)a.y * b.x;
}
Point operator-(const Point &a, const Point &b) {
return Point(a.x - b.x, a.y - b.y);
}
bool operator!=(const Point &a, const Point &b) {
return a.x != b.x || a.y != b.y;
}
int T, n, top, cnt, Mon_top;
Point point[maxn], sta[maxn], p[maxn], Mon_sta[maxn];
map<int, int> mp;
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(false);
scanf("%d", &T);
while(T--) {
mp.clear();
cnt = 0;
scanf("%d", &n);
for(int i = 1; i <= n; ++i) {
scanf("%d%d", &point[i].x, &point[i].y);
point[i].Index = i;
mp[point[i].x] = max(mp[point[i].x], point[i].y);
}
sort(point + 1, point + 1 + n);
p[0].x = p[0].y = -1;
for(int i = 1; i <= n; ++i) {
if(point[i].y == mp[point[i].x]) {
if(point[i] != p[cnt]) {
p[++cnt] = point[i];
} else if(point[i].Index < p[cnt].Index) {
p[cnt] = point[i];
}
}
}
top = 0;
for(int i = 1; i <= cnt; ++i) {
while(top > 1 && ((sta[top - 1] - sta[top - 2]) ^ (p[i] - sta[top - 2])) > 0) {
--top;
}
sta[top++] = p[i];
}
printf("1");
for(int i = 0; i < top - 1; ) {
int End = i + 1;
while(End + 1 < top && ((sta[End + 1] - sta[End]) ^ (sta[End] - sta[i])) == 0) {
++End;
}
Mon_top = 0;
for(int j = i + 1; j <= End; ++j) {
while(Mon_top > 0 && sta[j].Index < Mon_sta[Mon_top - 1].Index) {
--Mon_top;
}
Mon_sta[Mon_top++] = sta[j];
}
for(int j = 0; j < Mon_top; ++j) {
printf(" %d", Mon_sta[j].Index);
}
i = End;
}
printf("\n");
}
return 0;
}
H. Monster Hunter
题意
小 Q Q 在一棵 个节点的树上,他最初的位置在节点 1 1 处,除了 以外,其他节点都有一个怪兽,小 Q Q 需要沿着树上路径行走,每到达一个节点 ,就要与这个节点上的怪兽进行战斗,战斗将会损失 ai a i 点 HP H P ,战斗过程中若 HP H P 值小于 0 0 则游戏结束,战斗结束后若小 的 HP H P 值不小于 0 0 ,则他的 值将会增加 bi b i ,每个节点上的怪兽被消灭后,再次走到这个节点时就不需要再进行战斗,问如果要消灭树上所有的怪兽,小 Q Q 最初至少需要有多少点 值。
输入
第一行为一个整数 T (1≤T≤2000) T ( 1 ≤ T ≤ 2000 ) ,接下去有 T T 组数据,每组数据第一行为一个整数 ,接下去 n−1 n − 1 行每行两个整数 ai,bi (0≤ai,bi≤109) a i , b i ( 0 ≤ a i , b i ≤ 10 9 ) ,第 i i 行整数表示节点 上的战斗 HP H P 增减值,接下去 n−1 n − 1 行每行两个整数 ui,vi (1≤ui,vi≤n) u i , v i ( 1 ≤ u i , v i ≤ n ) ,表示节点 ui u i 和 vi v i 之间有一条边,数据保证 ∑n≤106 ∑ n ≤ 10 6 。
输出
输出通关所需要的最小的初始 HP H P 值。
样例
输入 |
---|
1 4 2 6 5 4 6 2 1 2 2 3 3 4 |
输出 |
3 |
题解
首先不考虑树上路径的限制,如果每个怪兽可以以任意次序进行战斗, a<b a < b 的怪兽必然要比 a≥b a ≥ b 的怪兽先进行攻击,这样才能获得最大的 HP H P ,在 a<b a < b 的所有怪兽中,必然要先攻击 a a 值小的,在 大的所有怪兽中,如果先攻击第 i i 只怪兽再攻击第 只怪兽,在攻击后一只怪兽的战斗过程中的 HP H P 值为 HP−ai+bi−aj H P − a i + b i − a j ,如果反过来先攻击第 j j 只怪兽,在后一只怪兽的战斗过程中的 值为 HP−aj+bj−ai H P − a j + b j − a i ,可以发现攻击的次序与 ai a i 和 aj a j 无关,而只与 bi b i 和 bj b j 有关,必然要先攻击 b b 值大的。
通过以上规则可以确定一个最优的攻击怪兽的次序 ,现在考虑树上路径的限制,攻击每只怪兽之前必须要先打败它所有的祖父节点,因此对于每个节点 pi p i ,判断该节点的父节点是否为 1 1 ,若为 ,则可以直接攻击,否则将 pi p i 的 ai,bi a i , b i 值合并到它的父节点的 af,bf a f , b f 上表示在攻击完它的父节点后必须立即攻击第 i i 个节点,父节点的 值应更新为
a′f=−min(−af,−af+bf−ai)b′f=−af+bf−ai+bi+a′f a f ′ = − min ( − a f , − a f + b f − a i ) b f ′ = − a f + b f − a i + b i + a f ′然后将所有 ai a i 的子节点的父节点改为 f f ,更新父节点并快速排序的过程可以用堆来维护,更改父节点的操作可以用并查集来维护。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;
#define LL long long
const int maxn = 100000 + 100;
struct Node {
LL a, b;
int fa, Index;
};
bool operator<(const Node &a, const Node &b) {
bool flag_a = a.a < a.b;
bool flag_b = b.a < b.b;
if(flag_a != flag_b) {
return flag_a < flag_b;
}
if(flag_a) {
return a.a > b.a;
}
return a.b < b.b;
}
bool operator!=(const Node &a, const Node &b) {
return a.a != b.a || a.b != b.b || a.fa != b.fa || a.Index != b.Index;
}
int T, n, u, v;
int fa[maxn];
Node node[maxn];
vector<int> G[maxn];
priority_queue<Node> que;
void Init() {
for(int i = 1; i <= n; ++i) {
fa[i] = i;
G[i].clear();
}
}
int Find(int x) {
return x == fa[x]? x: fa[x] = Find(fa[x]);
}
void unit(int x, int y) {
int xx = Find(x);
int yy = Find(y);
fa[xx] = yy;
}
void dfs(int f, int x) {
node[x].fa = f;
int len = G[x].size();
for(int i = 0; i < len; ++i) {
int pos = G[x][i];
if(pos != f) {
dfs(x, pos);
}
}
}
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(false);
scanf("%d", &T);
while(T--) {
scanf("%d", &n);
Init();
node[1].a = node[1].b = 0;
node[1].Index = 1;
for(int i = 2; i <= n; ++i) {
scanf("%I64d%I64d", &node[i].a, &node[i].b);
node[i].Index = i;
}
for(int i = 1; i < n; ++i) {
scanf("%d%d", &u, &v);
G[u].push_back(v);
G[v].push_back(u);
}
dfs(1, 1);
for(int i = 2; i <= n; ++i) {
que.push(node[i]);
}
while(!que.empty()) {
Node tmp = que.top();
que.pop();
if(node[tmp.Index] != tmp || tmp.Index == 1) {
continue;
} int f = Find(tmp.fa);
unit(tmp.Index, f);
LL a = -min(-node[f].a, -node[f].a + node[f].b - tmp.a);
LL b = -node[f].a + node[f].b - tmp.a + tmp.b + a;
node[f].a = a;
node[f].b = b;
que.push(node[f]);
}
printf("%I64d\n", node[1].a);
}
return 0;
}
I. Random Sequence
题意
有一个长度为 的序列 a1,a2,⋯,an (0≤ai≤m) a 1 , a 2 , ⋯ , a n ( 0 ≤ a i ≤ m ) 和长度为 m m 的序列 ,如果 ai=0 a i = 0 ,则 ai a i 的值可能为 1 1 到 之间的任意值,取得任意值的概率为 1m 1 m ,求公式:
∏i=1n−3vgcd(ai,ai+1,ai+2,ai+3) ∏ i = 1 n − 3 v gcd ( a i , a i + 1 , a i + 2 , a i + 3 )的期望值。
输入
第一行为一个整数 T (1≤T≤10) T ( 1 ≤ T ≤ 10 ) ,接下去有 T T 组数据,每组数据第一行为两个整数 ,第二行为 n n 个整数 ,第三行为 m m 个整数 。
输出
输出所求公式的期望值对 109+7 10 9 + 7 取模后的结果,若期望值无法用整数表示,则先将期望值化简为 AB A B 后,输出 C (0≤C<109+7) C ( 0 ≤ C < 10 9 + 7 ) 的值,其中 C C 满足 。
样例
输入 |
---|
2 6 8 4 8 8 4 6 5 10 20 30 40 50 60 70 80 4 3 0 0 0 0 3 2 4 |
输出 |
8000 3 |
题解
由于任意两个 v v 之间的选择是相互独立的,所以 。
如果定义状态 dp[i][x][y][z] d p [ i ] [ x ] [ y ] [ z ] 表示到第 i i 位 时 ∏ij=4vgcd(aj,aj−1,aj−2,aj−3) ∏ j = 4 i v gcd ( a j , a j − 1 , a j − 2 , a j − 3 ) 的所有合法情况的和,则可以通过枚举第 ai+1 a i + 1 位的所有合法值 xx x x 计算 dp[i+1][xx][x][y] d p [ i + 1 ] [ x x ] [ x ] [ y ] 的值,其递推公式为:
dp[i+1][xx][x][y]=∑x,y,z∈[1,m]dp[i][x][y][z]×vgcd(xx,x,y,z) d p [ i + 1 ] [ x x ] [ x ] [ y ] = ∑ x , y , z ∈ [ 1 , m ] d p [ i ] [ x ] [ y ] [ z ] × v gcd ( x x , x , y , z )注意当 ai a i 的值非 0 0 时 的值只能取 ai a i ,只有 ai a i 的值为 0 0 时 才可以从 1 1 到 枚举所有值,对 xx,y,z x x , y , z 同理。这种递推式的时间复杂度为 O(nm4) O ( n m 4 ) ,必然超时。
我们可以定义状态 dp[i][x][y][z] d p [ i ] [ x ] [ y ] [ z ] 表示到第 i i 位 时所有合法情况下 ∏ij=4vgcd(aj,aj−1,aj−2,aj−3) ∏ j = 4 i v gcd ( a j , a j − 1 , a j − 2 , a j − 3 ) 的和,由公约数关系可以得到, x,y,z x , y , z 必须满足 z z 能整除 , y y 整除 ,于是有下列递推式:
dp[i+1][xx][gcd(xx,x)][gcd(xx,y)]=∑x,y,z∈[1,m],z|y,y|xdp[i][x][y][z]×vgcd(xx,z) d p [ i + 1 ] [ x x ] [ gcd ( x x , x ) ] [ gcd ( x x , y ) ] = ∑ x , y , z ∈ [ 1 , m ] , z | y , y | x d p [ i ] [ x ] [ y ] [ z ] × v gcd ( x x , z )其中 xx,x,y,z x x , x , y , z 在 ai a i 非 0 0 时的取值限制同上, 表示 a a 能整除 ,在 z|y,y|x z | y , y | x 的限制下, x,y,z x , y , z 所有合法的状态可以减少到 1471 1471 种,因此时间复杂度为 O(1471×nm) O ( 1471 × n m ) ,其中 gcd gcd 的计算可以 O(m2logm) O ( m 2 log m ) 通过预处理将复杂度降低到 O(1) O ( 1 ) 。
最后是 dp d p 的初始值的确定:
dp[3][x][gcd(x,y)][gcd(x,y,z)]=∑x,y,z∈[1,m]1 d p [ 3 ] [ x ] [ gcd ( x , y ) ] [ gcd ( x , y , z ) ] = ∑ x , y , z ∈ [ 1 , m ] 1其中 x,y,z x , y , z 的取值限制同上, dp d p 初始值即表示取前 3 3 位为 的合法状态数,系数为 1 1 ,初始化方式为暴力,时间复杂度为 。
最后将 ∑x,y,z∈[1,m],z|y,y|xdp[n][x][y][z] ∑ x , y , z ∈ [ 1 , m ] , z | y , y | x d p [ n ] [ x ] [ y ] [ z ] 乘上总方案数对 109+7 10 9 + 7 的逆元就是答案,总方案数为 mcnt m c n t ,其中 cnt c n t 为 a a 序列中 的个数。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;
#define LL long long
const int maxn = 101;
const int MOD = 1000000007;
int T, n, m, ans, cnt, inv;
int a[maxn], v[maxn];
int Gcd[maxn][maxn], dp[maxn][maxn][maxn][maxn];
int add(int a, int b) {
a += b;
if(a >= MOD) {
return a - MOD;
}
if(a < 0) {
return a + MOD;
}
return a;
}
int gcd(int x, int y) {
return y == 0? x: gcd(y, x % y);
}
void prepare_gcd() {
for(int i = 1; i < maxn; ++i) {
Gcd[0][i] = Gcd[i][0] = i;
for(int j = 1; j < maxn; ++j) {
Gcd[i][j] = gcd(i, j);
}
}
}
void Init() {
for(int i = 3; i <= n; ++i) {
for(int z = 1; z <= m; ++z) {
for(int y = z; y <= m; y += z) {
for(int x = y; x <= m; x += y) {
dp[i][x][y][z] = 0;
}
}
}
}
for(int x = 1; x <= m; ++x) {
if(a[3] != 0 && x != a[3]) {
continue;
}
for(int y = 1; y <= m; ++y) {
if(a[2] != 0 && a[2] != y) {
continue;
}
for(int z = 1; z <= m; ++z) {
if(a[1] != 0 && z != a[1]) {
continue;
}
++dp[3][x][Gcd[x][y]][Gcd[Gcd[x][y]][z]];
}
}
}
}
LL fast_pow(LL res, LL n) {
LL ans;
for(ans = 1; n != 0; n >>= 1) {
if((n & 1) == 1) {
ans = (ans * res) % MOD;
}
res = (res * res) % MOD;
}
return ans;
}
int main() {
#ifdef Dmaxiya
freopen("test.txt", "r", stdin);
// freopen("test1.out", "w", stdout);
#endif // Dmaxiya
ios::sync_with_stdio(false);
prepare_gcd();
scanf("%d", &T);
while(T--) {
cnt = 0;
scanf("%d%d", &n, &m);
for(int i = 1; i <= n; ++i) {
scanf("%d", &a[i]);
if(a[i] == 0) {
++cnt;
}
}
for(int i = 1; i <= m; ++i) {
scanf("%d", &v[i]);
}
Init();
for(int i = 3; i < n; ++i) {
for(int z = 1; z <= m; ++z) {
for(int y = z; y <= m; y += z) {
for(int x = y; x <= m; x += y) {
if(dp[i][x][y][z] == 0) {
continue;
}
for(int xx = 1; xx <= m; ++xx) {
if(a[i + 1] != 0 && a[i + 1] != xx) {
continue;
}
dp[i + 1][xx][Gcd[xx][x]][Gcd[xx][y]] = add(dp[i + 1][xx][Gcd[xx][x]][Gcd[xx][y]], (LL)dp[i][x][y][z] * v[Gcd[xx][z]] % MOD);
}
}
}
}
}
ans = 0;
for(int k = 1; k <= m; ++k) {
for(int j = k; j <= m; j += k) {
for(int i = j; i <= m; i += j) {
ans = add(ans, dp[n][i][j][k]);
}
}
}
ans = (LL)ans * fast_pow(fast_pow(m, cnt), MOD - 2) % MOD;
printf("%d\n", ans);
}
return 0;
}
L. Visual Cube
题意
给定一个长方体的长宽高 a,b,c a , b , c ,按照样例输出长方体。
输入
第一行为一个整数 T (1≤T≤50) T ( 1 ≤ T ≤ 50 ) ,接下去 T T 行每行三个整数 。
输出
按样例输出。
样例
输入 |
---|
2 1 1 1 6 2 4 |
输出 |
题解
由于前面覆盖后面,上面覆盖下面,右边覆盖左边,所以可以依次从后往前,从下往上,从左往右输出网格。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;
#define LL long long
const int maxn = 100 + 100;
int T, a, b, c, row, col;
char str[maxn][maxn];
void Draw(int x, int y, bool back) {
for(int i = 2 * c; i >= 0; i -= 2) {
for(int j = 0; j <= 2 * a; j += 2) {
str[x + i][y + j] = '+';
if(!back) {
str[x + i - 1][y + j] = '.';
str[x + i - 1][y + j + 1] = '/';
str[x + i][y + j + 1] = '.';
}
if(i != 2 * c) {
str[x + i + 1][y + j] = '|';
}
if(j != 0) {
str[x + i][y + j - 1] = '-';
}
if(i != 2 * c && j != 0) {
str[x + i + 1][y + j - 1] = '.';
}
}
}
}
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(false);
scanf("%d", &T);
while(T--) {
scanf("%d%d%d", &a, &b, &c);
row = 2 * c + 1 + 2 * b;
col = 2 * a + 1 + 2 * b;
for(int i = 0; i < row; ++i) {
memset(str[i], '.', sizeof(char) * col);
str[i][col] = '\0';
}
for(int i = 0; i <= 2 * b; i += 2) {
int j = 2 * b - i;
bool back = false;
if(i == 0) {
back = true;
}
Draw(i, j, back);
}
for(int i = 0; i < row; ++i) {
printf("%s\n", str[i]);
}
}
return 0;
}
M. Walking Plan
题意
给一个 n n 个节点 条边的有向图,第 i i 条边的两个端点为 ,边的长度为 wi w i , q q 次询问,每次询问从节点 到 t t 至少走过 条路径的最小距离。
输入
第一行包含一个整数 T (1≤T≤10) T ( 1 ≤ T ≤ 10 ) ,接下去有 T T 组数据,每组数据第一行为两个整数 ,接下去 m m 行每行三个整数 ,接着为一个整数 q (1≤q≤105) q ( 1 ≤ q ≤ 10 5 ) ,接下去 q q 行每行三个整数 。
输出
对于每次询问,输出最短路径长度,如果无法从节点 si s i 到达 ti t i ,则输出 −1 − 1 。
样例
输入 |
---|
2 3 3 1 2 1 2 3 10 3 1 100 3 1 1 1 1 2 1 1 3 1 2 1 1 2 1 1 2 1 1 |
输出 |
111 1 11 -1 |
题解
定义 G[i][j] G [ i ] [ j ] 为从节点 i i 恰好经过 步到达节点 j j 的最短距离(即原图按输入取最小值),无法到达设为 , dis[k][i][j] d i s [ k ] [ i ] [ j ] 表示从 i i 点出发恰好经过 步到达 j j 点的最短路径,则有递推式:
定义 A[x][i][j] A [ x ] [ i ] [ j ] 表示从 i i 恰好经过 步到达 j j 的最短距离, 表示从 i i 最少经过 次到达 j j 的最短距离,这样 就可以分为 ⌊k100⌋ ⌊ k 100 ⌋ 和 k%100 k % 100 两个部分,通过枚举中转点得到答案:
ans=min(A[⌊k100⌋][i][u]+B[k%100][u][j]),u∈[1,n] a n s = min ( A [ ⌊ k 100 ⌋ ] [ i ] [ u ] + B [ k % 100 ] [ u ] [ j ] ) , u ∈ [ 1 , n ]第一部分可以直接从上面的递推式得到,第二部分可以在原图上令 G[i][i]=0,i∈[1,n] G [ i ] [ i ] = 0 , i ∈ [ 1 , n ] 后跑一遍 floyd f l o y d ,再将从 i i 恰好经过 步到达 j j 的最短距离距离 通过 G[i][j] G [ i ] [ j ] 转化为 B[x][i][j] B [ x ] [ i ] [ j ] :
B[k][i][j]=min(dis[k][i][u],G[u][j]),u∈[1,n] B [ k ] [ i ] [ j ] = min ( d i s [ k ] [ i ] [ u ] , G [ u ] [ j ] ) , u ∈ [ 1 , n ]
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;
#define LL long long
const int Size = 51;
const int maxn = 101;
int n;
struct Matrix {
int num[Size][Size];
void Init() {
for(int i = 1; i <= n; ++i) {
memset(num[i], 0x3f, sizeof(int) * (n + 1));
}
}
void Set_zero() {
for(int i = 1; i <= n; ++i) {
num[i][i] = 0;
}
}
void operator=(const Matrix &m) {
for(int i = 1; i <= n; ++i) {
memcpy(num[i], m.num[i], sizeof(int) * (n + 1));
}
}
void Combine(const Matrix &m, Matrix &ans) {
Matrix ret;
ret.Init();
for(int i = 1; i <= n; ++i) {
for(int j = 1; j <= n; ++j) {
for(int k = 1; k <= n; ++k) {
ret.num[i][j] = min(ret.num[i][j], num[i][k] + m.num[k][j]);
}
}
}
ans = ret;
}
void floyd() {
for(int k = 1; k <= n; ++k) {
for(int i = 1; i <= n; ++i) {
for(int j = 1; j <= n; ++j) {
num[i][j] = min(num[i][j], num[i][k] + num[k][j]);
}
}
}
}
};
int T, m, u, v, dis, q, k, INF, ans;
Matrix A[maxn], B[maxn], G;
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("test1.out", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(false);
memset(&INF, 0x3f, sizeof(int));
scanf("%d", &T);
while(T--) {
scanf("%d%d", &n, &m);
G.Init();
for(int i = 0; i < m; ++i) {
scanf("%d%d%d", &u, &v, &dis);
G.num[u][v] = min(G.num[u][v], dis);
}
A[0].Init();
A[0].Set_zero();
B[0].Init();
B[0].Set_zero();
for(int i = 1; i < maxn; ++i) {
B[i - 1].Combine(G, B[i]);
}
for(int i = 1; i < maxn; ++i) {
A[i - 1].Combine(B[100], A[i]);
}
G.Set_zero();
G.floyd();
for(int i = 0; i < maxn; ++i) {
G.Combine(B[i], B[i]);
}
scanf("%d", &q);
while(q--) {
scanf("%d%d%d", &u, &v, &k);
ans = INF;
for(int i = 1; i <= n; ++i) {
ans = min(ans, A[k / 100].num[u][i] + B[k % 100].num[i][v]);
}
if(ans == INF) {
printf("-1\n");
} else {
printf("%d\n", ans);
}
}
}
return 0;
}