2021西安邀请赛部分题解

最新推荐文章于 2024-07-13 15:34:59 发布

敲代码的欧文

最新推荐文章于 2024-07-13 15:34:59 发布

阅读量1k

点赞数

本文链接：https://blog.csdn.net/Irving0323/article/details/118391586

版权

G

题意：一棵树每个点权有一个 $w_i$ ， $z$ 可以作为 $u,v$ 的中转站的条件如下。

$z$ 在 $u,v$ 之间的路径上。
$z$ 可以是 $u,v$ 之一。
$u,v$ 需要不同。
$w_u - w_z \leqslant p$
$w_v - w_z \leqslant p$

定义点 $z$ 的承载量是 $z$ 可以作为 $u,v$ 的中转站的无序对 $(u,v)$ 的个数。

最大化 $p$ 使得每个点的承载量小于 $k$ 。

题解：二分 $p$ 。将 $w$ 升序，按顺序依次将每个点塞到按照dfs序建立的树状数组上。

有一个坑点：需要根据 $z$ 是 $u,v$ 之一， $z$ 不是 $u,v$ 之一进行分类讨论，需要考虑清楚 $z$ 是 $u,v$ 之一这个case。

H

题意：有 $n$ 种宝石，一次魔法可以把第 $i$ 种宝石变成第 $j$ 种宝石并且个数变成原来的 $a[i][j]$ 倍。每次询问最少需要多少次魔法可以把 $1$ 个第 $x$ 种宝石变成 $2^k$ 个第种宝石。

数据范围：询问次数是 $1000$ 次。 $1\leqslant n \leqslant 100,1\leqslant k \leqslant 10^5,0\leqslant a[i][j]\leqslant 10$ 。

题解：

先用一个对数trick。

考虑矩阵快速幂 $A[t]][i][j]$ 表示进行 $2^t$ 次魔法后第 $j$ 种宝石到第 $i$ 种宝石的最大翻倍数。

每次询问需要从高位到低位拼出答案。用 $A[t]$ 和向量 $res$ 相乘。 $res$ 的大小是 $n*1$ 。 $res[i]$ 表示当前次数下第 $x$ 种宝石到第 $i$ 种宝石的最大翻倍数。需要保证当前次数下第 $x$ 种宝石的翻倍数小于 $2^k$ 。拼出的二进制数是 $ans$ 。答案是 $ans+1$ 即翻倍数不小于 $2^k$ 。

复杂度 $O(1000*20*n^2)$ 。

感想：确实是个好题，贴个代码。

#include<bits/stdc++.h>
using namespace std ;
typedef double db ;
typedef vector<db> vec ;
typedef vector<vec> mat ;
mat mul(mat &A , mat &B)
{
    mat C(A.size() , vec(B[0].size() , -1e9)) ;
    for(int i = 0 ; i < (int)A.size() ; i ++)
        for(int j = 0 ; j < (int)B[0].size() ; j ++)
            for(int k = 0 ; k < (int)A[0].size() ; k ++)
                C[i][j] = max(C[i][j] , A[i][k] + B[k][j]) ;
    return C ;
}
int main()
{
    std::ios::sync_with_stdio(false) , cin.tie(0) ;
    int n , q ;
    cin >> n >> q ;
    mat a(n , vec(n)) ;
    vector<mat> A(20 , mat(n , vec(n))) ;
    for(int i = 0 ; i < n ; i ++)
        for(int j = 0 ; j < n ; j ++)
            cin >> a[i][j] ;
    for(int i = 0 ; i < n ; i ++)
        for(int j = 0 ; j < n ; j ++)
        {
            if(a[j][i] == 0)  A[0][i][j] = -1e9 ;
            else  A[0][i][j] = log2(a[j][i]) ;
            if(i == j)  A[0][i][j] = 0 ;
        }
    for(int i = 1 ; i <= 19 ; i ++)  A[i] = mul(A[i - 1] , A[i - 1]) ;
    while(q --)
    {
        int x , k ;
        cin >> x >> k ;
        x -- ;
        mat res(n , vec(1 , -1e9)) ;
        res[x][0] = 0 ;
        int ans = 0 ;
        for(int i = 19 ; i >= 0 ; i --)
        {
            mat tmp = mul(A[i] , res) ;
            if(tmp[x][0] < k)  ans += (1 << i) , res = tmp ;
        }
        cout << ans + 1 << '\n' ;
    }
    return 0 ;
    //you should actually read the stuff at the bottom before submitting or in the confusion
}
/* stuff you should look for
    * long long
    * array bounds
    * init
    * ios
    * special cases (n=1?)
    * do smth instead of nothing and stay organized
    * WRITE STUFF DOWN
    * DON'T GET STUCK ON ONE APPROACH
    * DON'T GET STUCK ON ONE PROBLEM
*/