【解题报告】2021牛客寒假算法基础集训营 6

最新推荐文章于 2022-01-25 19:15:51 发布

溢流眼泪

最新推荐文章于 2022-01-25 19:15:51 发布

阅读量266

点赞数 2

分类专栏：【解题报告 / 补题报告】文章标签：算法

本文链接：https://blog.csdn.net/weixin_45775438/article/details/114062116

版权

【解题报告 / 补题报告】专栏收录该内容

32 篇文章 6 订阅

订阅专栏

【解题报告】2021牛客寒假算法基础集训营 6

前面的话
$\mathfrak{A}$ ：回文括号序列计数 | QD思维
$\mathfrak{B}$ ：系数 | 数学
$\mathfrak{C}$ ：末三位 | QD数学
$\mathfrak{D}$ ：划数 | QD数学
$\mathfrak{E}$ ：网格 | $d p$
$\mathfrak{F}$ ：组合数问题 | 数学 + 打表
$\mathfrak{G}$ ：机器人 | 状压 $d p$ + 高精度
$\mathfrak{H}$ ：动态最小生成树
$\mathfrak{I}$ ：贪吃蛇 | QD $B F S$
$\mathfrak{J}$ ：天空之城 | $K r u s k a l$ $M S T$

前面的话

2021牛客寒假算法基础集训营 6
这场真的是 $\lceil$ 命途多舛 $\rfloor$ ， $A$ 坑题， $H$ 暴力能过，签到题还挺多…
以后我签到题写个 $Q D$ 吧，就只给个题意和思路吧，不然长度太长不容易翻阅。

$\mathfrak{A}$ ：回文括号序列计数 | QD思维

【题意】
问长度为 $n$ 的回文合法括号序列的个数
【范围】
样例组数 $T\le 10^6$
$0\le n\le 10^9$
【思路】赛内
问：以下哪个括号串是回文的：
A： $(() ())$
B： $(() (()$
C： $(()) (($
答案是 $\tiny{c}$ ， $A$ 是镜像串 $Symmetry\ String$ 。
那就送分题了。第一个位置为 $)$ 或者 $($ 都会导致结果不合法，只有 $n = 0$ 时才为 $1$ （空串）
顺带一提，如果题目问的是镜像串的数量，那么首先判断 $n$ 必须为偶数，然后 $n = n / 2$ ，最终答案为 $C_n^{\lfloor n/2\rfloor}$

$\mathfrak{B}$ ：系数 | 数学

【题意】
问你 $x^2+x+1)^n$ 的第 $k$ 项系数取模 $3$ 是多少。
【范围】
样例组数 $T\le 10^4$
$0\le n\le 10^{15}$
$0\le k\le 2\times n$
【思路】题解 + 赛后
如果直接去做的话，是一个 $O (n k)$ 的 $d p$ ，式子如下：
$\begin{cases} 1& j=0 或者 j=2i\\ dp[i-1][j-1]+dp[i-1][j]&j=1\\ dp[i-1][j-2]+dp[i-1][j-1]&j=2i-1\\ dp[i-1][j-2]+dp[i-1][j-1]+dp[i-1][j]&Others \end{cases}$
~~应该没写错吧~~
这是一个不规则三角阵，貌似没法直接去做…有题解看出来是一个分形，进行具体分析的 $\%\%\%$
我们注意到一个细节： $输出取模 3 的系数$ ，原式子为 $x^2+x+1)^n$
考虑同余，我们知道 $(x^2+x+1)^n\equiv (x^2-2x+1)^n\equiv (x-1)^{2n}\pmod 3$
然后直接二项式，答案就得到了：
$1)^{2n-k}C_{2n}^k$
因为模数为 $3$ 是一个质数，直接用 $L u c a s$ 定理做就行了。
【代码】
时间复杂度： $T\log n$
$32 / 1000 M s$

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
const int MAX = 1e5+50;

ll MOD;
ll qpow(ll a,ll n){/* */ll res = 1LL;while(n){if(n&1)res=res*a%MOD;a=a*a%MOD;n>>=1;}return res;}
ll inv(ll a){/* */return qpow(a,MOD-2);}

ll fac[MAX],ivfac[MAX];
void init(int n){
    fac[0] = 1;
    for(int i = 1;i <= n;++i)fac[i] = fac[i-1] * i % MOD;
    ivfac[n] = inv(fac[n]);
    for(int i = n-1;i >= 0;--i)ivfac[i] = ivfac[i+1]*(i+1)%MOD;
}
ll C(ll n,ll m){
    if(m > n)return 0;
    return fac[n] * ivfac[m] % MOD * ivfac[n - m] % MOD;
}
ll Lucas(ll n,ll m){
    if(m == 0)return 1;
    return C(n % MOD,m % MOD) * Lucas(n / MOD, m / MOD) % MOD;
}
int main()
{
    ll n,k;
    int T;cin >> T;
    while(T--){
        cin >> n >> k;
        MOD = 3;
        init(MOD - 1);
        int sgn = 1;
        if((2*n-k)&1)sgn = -1;
        ll res = (Lucas(2*n,k) * sgn % MOD + MOD) % MOD;
        cout << res << endl;
    }
    return 0;
}

$\mathfrak{C}$ ：末三位 | QD数学

【题意】
$5^n$ 的末三位为多少
【范围】
样例组数 $T\le 10^6$
$0\le n\le 10^9$
【思路】
对于任意 $x$ ， $x^n$ 的末三位都有循环节的。 $x = 5$ 的循环节如下：
$005,025,125,625,125,625\cdots$

$\mathfrak{D}$ ：划数 | QD数学

【题意】
有 $n$ 个数，每次划去两个数 $x ， y$ ，增加一个数 $(x+y)\%11$
划到只剩两个数时，一个数为 $c n t$ ，问你另个一数为几。
【范围】
$cnt\ge11$
$2\le n\le 15000$
$1\le num[i]\le 10^6$
【思路】
剩下的那个 $c n t$ 肯定是没划过的。
$n = 2$ ，剩下的就是另一个数字。
$n > 2$ ，剩下的就是其他所有划过的数字。等于其他数字的总和取模 $11$ 。

$\mathfrak{E}$ ：网格 | $d p$

【题意】
$n\times m$ 的网格，每个位置有一个权值 $a_{i,j}$ ，每个位置要从上下左右方向选择互相垂直的两个。
定义 $w (x) = x + p o p c n t (x)$ ，其中 $p o p c n t (x)$ 表示 $x$ 的二进制表示中 $1$ 的个数。
如果两个相邻的位置 $x_1,y_1)$ 和 $x_2,y_2)$ 互相位与对方选择的方向之中，那么答案增加 $w(a_{x_1,y_1}\oplus a_{x_2,y_2})$ ，其中 $\oplus$ 表示二进制按位异或。
问你答案最大值为多少。
【范围】
$1\le n,m\le 1000$
$0\le a_{i,j}<1024$
【思路】题解
赛内的时候感觉：方向好复杂！ $w (x)$ 好复杂！贡献还要异或！一堆给人不可做的感觉…
但是我们把要求 $\lceil$ 抽丝剥茧 $\rfloor$ 之后，就没有那么复杂了。
方向选择的要求是上下左右方向选择互相垂直的两个，我们改成：向左和向右选择一个方向，向上和向下选择一个方向。
先考虑左右方向。那么对答案的贡献只能是一行里面某两个相邻位置对答案有贡献。
这不就是最简单的 $d p$ 吗？设 $d p [i] [0]$ 表示当前行第 $i$ 列选择方向向左， $d p [i] [1]$ 表示选择方向向右。
$\begin{aligned} dp[i][1]&=\max\{dp[i-1][0],dp[i-1][1]\}\\ dp[i][0]&=\max\{dp[i-1][0],dp[i-1][1]+w(贡献)\} \end{aligned}$
这个贡献，我们只要预处理好每个数字的 $p o p c n t [i]$ 即可。
对于上下方向，和该转移式子差别不大。
【代码】
时间复杂度： $O (n m)$
$19 / 3000 M s$

const int MAX = 1e3+50;

int pop[MAX];
int aa[MAX][MAX];
ll dp[MAX][2];
int lowbit(int x){
    return (x & -x);
}
void init(int n){
    pop[0] = 0;
    for(int i = 1;i <= n;++i){
        pop[i] = pop[i-lowbit(i)] + 1;
    }
}
int main()
{
    init(1024);
    int n,m;n = read();m = read();
    for(int i = 1;i <= n;++i)
        for(int j = 1;j <= m;++j)
            aa[i][j] = read();
    ll res = 0;
    for(int i = 1;i <= n;++i){
        dp[1][0] = dp[1][1] = 0;
        for(int j = 2;j <= m;++j){
            ll tmp = aa[i][j] ^ aa[i][j-1];
            ll duo = tmp + pop[tmp];
            dp[j][1] = max(dp[j-1][0],dp[j-1][1]);
            dp[j][0] = max(dp[j-1][0],dp[j-1][1] + duo);
        }
        res += dp[m][0];
    }
    for(int i = 1;i <= m;++i){
        dp[1][0] = dp[1][1] = 0;
        for(int j = 2;j <= n;++j){
            ll tmp = aa[j][i] ^ aa[j-1][i];
            ll duo = tmp + pop[tmp];
            dp[j][1] = max(dp[j-1][0],dp[j-1][1]);
            dp[j][0] = max(dp[j-1][0],dp[j-1][1] + duo);
        }
        res += dp[n][0];
    }
    printf("%lld",res);
    return 0;
}

$\mathfrak{F}$ ：组合数问题 | 数学 + 打表

【题意】
$4 ∣ n$ ，求计算 $C_n^0+C_n^4+C_n^8+\cdots+C_n^n$
【范围】
$1\le n\le 10^{18}$
【思路】赛内
前几项打表，得出式子： $t = n / 4$
$O E I S$ 之后得出答案：
$ans=(-4)^t\times \frac{1}{2}+16^n\times \frac{1}{4}$
~~虽然有点偷懒~~

$\mathfrak{G}$ ：机器人 | 状压 $d p$ + 高精度

【题意】
$n$ 个机器，第 $i$ 个入读一个数字 $x$ 就会输出 $a_ix+b_i$
你一开始有一个数字 $x$ ，每个机器只用一次，最终最大值为多少。
【范围】
$1\le a_i,b_i,x,n\le 20$
【思路】赛内
$20^{20}$ 爆 $l l$ 了，我最后用了 $J a v a$ 大数写，虽然也可以 $_\_int128$
一个简单的状压 $d p$ ，假设现在的机器人集合为 $i$ ， $d p [i]$ 表示选择完这个集合的机器人，输出的最大值
$dp[i]=\max\{dp[i-j]\times a_j+b_j\}\qquad i\&j=j,j=lowbit(j)$
很好懂的吧…
【代码】
时间复杂度： $O(2^{n}\times n)$
$2645 / 6000 M s$

import java.math.BigInteger;
import java.util.Scanner;

public class Main {
    static BigInteger[] A = new BigInteger[30];
    static BigInteger[] B = new BigInteger[30];
    static BigInteger[] dp = new BigInteger[1100000];
    static int[] pw = new int[30];
    public static void main(String[] args){
        Scanner read = new Scanner(System.in);
        int n = read.nextInt();
        BigInteger x = read.nextBigInteger();
        int ed = 1;
        for(int i = 0;i < n;++i){
            A[i] = read.nextBigInteger();
            B[i] = read.nextBigInteger();
            pw[i] = ed;
            ed = ed * 2;
        }
        ed--;
        dp[0] = x;
        for(int i = 1;i <= ed;++i){
            dp[i] = BigInteger.ZERO;
            for(int j = 0;j < n;++j){
                if((i & pw[j]) > 0){
                    int re = i - pw[j];
                    BigInteger tmp = dp[re].multiply(A[j]);
                    tmp = tmp.add(B[j]);
                    dp[i] = dp[i].max(tmp);
                }
            }
        }
        System.out.println(dp[ed]);
    }
}

【补充】
这题也可以贪心排序，不过这个数据，出题人的本意是让我们做状压 $d p$ 的。

$\mathfrak{H}$ ：动态最小生成树

赛后看别人代码：暴力能过？暴力能过！暴力能过…
正解显然我不会。

$\mathfrak{I}$ ：贪吃蛇 | QD $B F S$

【题意】
走迷宫，蛇每走一格，就长长一格。
【范围】
$1\le n,m\le 100$
【思路】赛内
走迷宫，就是无法走回头路，还是 $B F S$ 就能搞定。
【补充】
这题无法到达输出 $- 1$

$\mathfrak{J}$ ：天空之城 | $K r u s k a l$ $M S T$

【题意】
有 $n$ 个城市， $q$ 条边。每个边的权值表示走该路的花费时间。
给定一个起点 $s$ 。你希望走过所有城市，且花费时间最短。
你走过的边再走就会不花时间。
【范围】
$1\le n\le 5000$
$1\le q\le 200000$
$1\le val\le 10^9$
【思路】赛内
城市名字字符串转成编号。
走过的边再走不会花时间？那就是一个 $M S T$ 的裸题。
写个并查集然后用 $K r u s k a l$ 即可。
【代码】
时间复杂度： $O(q\log q)$
$425 / 5000 M s$

/*
 _            __   __          _          _
| |           \ \ / /         | |        (_)
| |__  _   _   \ V /__ _ _ __ | |     ___ _
| '_ \| | | |   \ // _` | '_ \| |    / _ \ |
| |_) | |_| |   | | (_| | | | | |___|  __/ |
|_.__/ \__, |   \_/\__,_|_| |_\_____/\___|_|
        __/ |
       |___/
*/
const int MAX = 2e5+50;

unordered_map<string,int>M;
struct node{
    int ta,tb;
    ll w;
    bool operator <(const node &ND)const{
        return w < ND.w;
    }
}aa[MAX];
int fa[MAX];
int find_fa(int x){
    if(x == fa[x])return x;
    return fa[x] = find_fa(fa[x]);
}
bool sam(int x,int y){
    int fx = find_fa(x);
    int fy = find_fa(y);
    return fx == fy;
}
void add(int x,int y){
    int fx = find_fa(x);
    int fy = find_fa(y);
    if(fx != fy){
        fa[fx] = fy;
    }
}

int main()
{
    int n,m;
    while(~scanf("%d%d",&n,&m)){
        int cnt = 0;
        M.clear();
        for(int i = 1;i <= n;++i){
            fa[i] = i;
        }
        string st;cin >> st;
        for(int i = 1;i <= m;++i){
            string ta,tb;ll q;
            cin >> ta >> tb >> q;
            int sa,sb;
            if(!M[ta])M[ta] = ++cnt;
            if(!M[tb])M[tb] = ++cnt;
            sa = M[ta];
            sb = M[tb];
            aa[i].ta = sa;
            aa[i].tb = sb;
            aa[i].w = q;
        }
        sort(aa+1,aa+1+m);
        ll res = 0;
        for(int i = 1;i <= m;++i){
            int ta = aa[i].ta;
            int tb = aa[i].tb;
            if(!sam(ta,tb)){
                add(ta,tb);
                res += aa[i].w;
            }
        }
        bool can = true;
        for(int i = 2;i <= n;++i){
            if(!sam(i,1)){
                can = false;
                break;
            }
        }
        if(!can)puts("No!");
        else printf("%lld\n",res);
    }
    return 0;
}

【补充】
这题无法到达输出 $N o!$

溢流眼泪

关注

2
点赞
踩
0

收藏

觉得还不错? 一键收藏
1
评论
【解题报告】2021牛客寒假算法基础集训营 6

【解题报告】2021牛客寒假算法基础集训营 6前面的话A\mathfrak{A}A ：回文括号序列计数 | QD思维B\mathfrak{B}B ：系数 | 数学C\mathfrak{C}C ：末三位 | QD数学D\mathfrak{D}D ：划数 | QD数学E\mathfrak{E}E ：网格 | dpdpdpF\mathfrak{F}F ：组合数问题 | 数学 + 打表G\mathfrak{G}G ：机器人 | 状压 dpdpdp + 高精度H\mathfrak{H}H ：动态最小生成树I\m
复制链接

扫一扫