7/28-7/29 期望+思维+后缀数组+ST表

最新推荐文章于 2024-07-25 13:44:24 发布
最新推荐文章于 2024-07-25 13:44:24 发布
阅读量164 收藏
点赞数 1
分类专栏： 22暑期ACM训练 文章标签： 算法 图论 c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_51934288/article/details/126047387
版权

期望的概念：
E(X+Y)=E(X)+E(Y) 因此可根据线性的可加性做题
在条件概率中，一个事件发生的概率固定

每日一题

Game (20上海ICPC热身赛)

类型：期望题，未与dp、数学结合
思路：
1.长度为n的序列中，存在多少互质的组合cnt
2.将n分成奇偶判断，若为奇数，则会产生（n-1）/2次消去，每次概率为cnt/C(n,2);
若为偶数，则会产生n/2次消去，每次概率为cnt/C(n,2)

#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)

using namespace std;
const int N=5e5+5;
const int mod=1e9+7;
int n;
int fac[N];
int fastpow(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
int getinv(int a) //求逆元
{
    return fastpow(a,mod-2)%mod;
}
int C(int n,int m)  //C(n,m)
{
    return ((fac[n]*getinv(fac[m])%mod)*(getinv(fac[n-m])%mod))%mod;
}

signed main()
{
    IOS;
    cin>>n;
    int cnt=0;
    for(int i=1;i<n;i++)
    {
        for(int j=i+1;j<=n;j++)
            if(__gcd(i,j)==1)   cnt++;
    }
    if(n%2)
    {
        int g=__gcd(cnt,n);
        cout<<cnt/g<<"/"<<n/g<<endl;
    }
    else
    {
        int g=__gcd(cnt,n-1);
        cout<<cnt/g<<"/"<<(n-1)/g<<endl;
    }
    return 0;
}

Music Game 五一集训派对day2

类型：数学+期望的结合
思路：
1.处理出任意段事件发生的概率；在处理出不同长度m幂次方的值。
2.根据时间期望的可加性，枚举不同长度发生的期望值进行累加。

#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)

using namespace std;
const int N=5e3+5;
const int mod=1e9+7;
int n,m,p[N],mul[N],pre[N][N];
int fastpow(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1) res=res*a%mod;
        a=a*a%mod;
        b>>=1;
    }
    return res;
}
int getinv(int a) //求逆元
{
    return fastpow(a,mod-2)%mod;
}

signed main()
{
    IOS;
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        cin>>p[i],mul[i]=fastpow(i,m);
    int inv=getinv(100);
    p[0]=p[n+1]=0;
    for(int i=1;i<=n;i++)   //不同区间段成功的概率
    {
        pre[i][i]=p[i]*inv%mod;
        for(int j=i+1;j<=n;j++)
            pre[i][j]=pre[i][j-1]*(p[j]*inv%mod)%mod;
    }
    int ans=0;
    for(int i=0;i<=n;i++)
    {
        for(int j=i+2;j<=n+1;j++)
        {
            ans+=pre[i+1][j-1]*mul[j-i-1]%mod*(100-p[i])%mod*inv%mod*(100-p[j])%mod*inv%mod;
            ans%=mod;
        }
    }
    cout<<ans<<endl;
    return 0;
}

后缀数组

P4051 [JSOI2007]字符加密

思路：将字符串扩大两倍，在进行后缀数组的排序。证明与题意处理方式等价：
1.abcab中ab在abcab前面，题意中ababc在abcab前面；若扩大两倍为abcababcab，则ababcab在abcababcab前面，后面多余并不影响
2.转为为了sa的裸题

#include <bits/stdc++.h>
#define int long long
#define endl '\n'

using namespace std;
const int N=1e6+5;
const int mod=1e9+7;
int n,k;
int rk[N],rk2[N];   //以i开头后缀的排名
char s[N];
int sa[N];   //表示sa[i]表示排名i的后缀的开头下标

//求解各个以i为起始下标的后缀字符串的排名
bool cmp(int i,int j)
{
    if(rk[i]!=rk[j])
        return rk[i]<rk[j];
    int ri=(i+k<=n ? rk[i+k]:-1);
    int rj=(j+k<=n ? rk[j+k]:-1);
    return ri<rj;
}
void getsa(int n,char *str)
{
    for(int i=1;i<=n;i++)
        sa[i]=i,rk[i]=s[i];  //利用ASCLL码
    for(k=1;k<=n;k*=2)
    {
        sort(sa+1,sa+1+n,cmp);
        rk2[sa[1]]=1;
        for(int i=2;i<=n;i++)
            rk2[sa[i]]=rk2[sa[i-1]]+cmp(sa[i-1],sa[i]);
        for(int i=1;i<=n;i++)
            rk[i]=rk2[i];
    }
}

int ht[N];   //维护数组rk相邻两个后缀的lcp(i-1和i的最长公共前缀)
void getht(int n,char *s)
{
    for(int i=1;i<=n;i++)
        rk[sa[i]]=i;
    int h=0;
    ht[1]=0;
    for(int i=1;i<=n;i++)
    {
        int j=sa[rk[i]-1];
        if(h>0)
            h--;
        for(;j+h<=n&&i+h<=n;h++)
            if(s[j+h]!=s[i+h])
                break;
        ht[rk[i]]=h;
    }
}
/*
int top,ll[N],rr[N],sk[N],ans;
int cal(int n,char *s)
{
    getsa(n,s);
    getht(n,s);
    top=1,sk[1]=1;
    for(int i=2;i<=n;i++)
    {
        while(top&&ht[sk[top]]>ht[i])
            rr[sk[top]]=i,top--;
        ll[i]=sk[top];
        sk[++top]=i;
    }
    while(top)
        rr[sk[top]]=n+1,top--;
    int res=0;
    for(int i=2;i<=n;i++)
        res+=(i-ll[i])*(rr[i]-i)*ht[i];
    return res;
}*/

signed main()
{
    cin>>(s+1);
    n=strlen(s+1);
    for(int i=1;i<=n;i++)
        s[i+n]=s[i];
    n*=2;
    getsa(n,s);
    for(int i=1;i<=n;i++)
        if(sa[i]<=n/2)
        cout<<s[sa[i]+n/2-1];
    cout<<endl;
	return 0;
}

思维

A. Color the Picture

思路：
1.要保证一个颜色至少三个相邻和它有相同的颜色，则保证该颜色至少能涂两列
2.若m为技术的情况，则需找到至少一种颜色能涂三列
3.n和m的数值调换的情况也要考虑到

#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)

using namespace std;
const int N=5e5+5;
const int mod=1e9+7;
int n,m,k,a[N];
int f,f1,ans;
void check(int n,int m)
{
    ans=0;
    for(int i=1;i<=k;i++)
    {
        if(a[i]/n>=2)
        {
            if(a[i]/n>=3)
                f1=1;
            ans+=a[i]/n;
        }
    }
    if(m%2)
    {
        if(ans>=m&&f1)
            f=1;
    }
    else
        if(ans>=m)
            f=1;
}
signed main()
{
    IOS;
    int t;cin>>t;
    while(t--)
    {
        f=f1=0;
        cin>>n>>m>>k;
        for(int i=1;i<=k;i++)
            cin>>a[i];
        check(n,m);
        swap(n,m);
        check(n,m);
        if(f)
            cout<<"YES"<<endl;
        else
            cout<<"NO"<<endl;
    }
    return 0;
}

ST表
1.基于倍增思想，可做到O(n*logn)的预处理，O(1)的查询
2.解决指定运算的区间询问
模板：

//ST表
int lg[N];  //记录log2N的整数值
int st[N][30];  //以i为左端点长度为2的j次方的序列所求最值
void ST(int n)
{
    lg[1]=0;
    for(int i=2;i<=n;i++)
        lg[i]=lg[i-1]+(1<<(lg[i-1]+1)==i);
    for(int i=1;i<=n;i++)
        st[i][0]=a[i];
    for(int j=1;(1<<j)<=n;j++)
    {
        for(int i=1;i+(1<<j)-1<=n;i++)
            st[i][j]=max(st[i][j-1],st[i+(1<<j-1)][j-1]);
    }
}
int sch(int l,int r)
{
    if(l>r)
        return 0;
    int k=lg[r-l+1];
    return max(st[l][k],st[r-(1<<k)+1][k]);
}

D. Rorororobot

题意：n行m列的网格中，每列下方中a[i]个格子被阻塞，可以往上下左右四个方向移动，不可出边界或碰到阻塞格子，不计次数，能否正好移动到重点。

思路：
1.先判断起点和终点间的行可否移动到同一行，每次移动需走k步
2.同理，列格子可否移动到同一列
3.每一列阻塞格子的最大值不可超过最终到达的起点值

#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
#define y1 y11
using namespace std;
const int N=5e5+5;
const int mod=1e9+7;
int n,m,a[N];
int x1,y1,x2,y2,k;

//ST表
int lg[N];  //记录log2N的整数值
int st[N][30];  //以i为左端点长度为2的j次方的序列所求最值
void ST(int n)
{
    lg[1]=0;
    for(int i=2;i<=n;i++)
        lg[i]=lg[i-1]+(1<<(lg[i-1]+1)==i);
    for(int i=1;i<=n;i++)
        st[i][0]=a[i];
    for(int j=1;(1<<j)<=n;j++)
    {
        for(int i=1;i+(1<<j)-1<=n;i++)
            st[i][j]=max(st[i][j-1],st[i+(1<<j-1)][j-1]);
    }
}
int sch(int l,int r)
{
    if(l>r)
        return 0;
    int k=lg[r-l+1];
    return max(st[l][k],st[r-(1<<k)+1][k]);
}
signed main()
{
    IOS;
    cin>>n>>m;
    for(int i=1;i<=m;i++)
        cin>>a[i];
    ST(m);
    int q;cin>>q;
    while(q--)
    {
        cin>>x1>>y1>>x2>>y2>>k;
        x1+=(n-x1)/k*k;
        x2+=(n-x2)/k*k;
        if(y1>y2)   swap(y1,y2);
        if(x1==x2&&abs(y2-y1)%k==0&&sch(y1,y2)<x1)
            cout<<"YES"<<endl;
        else
            cout<<"NO"<<endl;
    }
    return 0;
}

P3865 【模板】ST 表

ST表静态查询，ST表知识优化查询的工具，一般不会单独考，会结合LCA等算法考察

#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define IOS ios::sync_with_stdio(false),cin.tie(0),cout.tie(0)
#define y1 y11
using namespace std;
const int N=5e5+5;
const int mod=1e9+7;
int n,m,a[N];
//ST表
int lg[N];  //记录log2N的整数值
int st[N][30];  //以i为左端点长度为2的j次方的序列所求最值
void ST(int n)
{
    lg[1]=0;
    for(int i=2;i<=n;i++)
        lg[i]=lg[i-1]+(1<<(lg[i-1]+1)==i);
    for(int i=1;i<=n;i++)
        st[i][0]=a[i];
    for(int j=1;(1<<j)<=n;j++)
    {
        for(int i=1;i+(1<<j)-1<=n;i++)
            st[i][j]=max(st[i][j-1],st[i+(1<<j-1)][j-1]);
    }
}
int sch(int l,int r)
{
    if(l>r)
        return 0;
    int k=lg[r-l+1];
    return max(st[l][k],st[r-(1<<k)+1][k]);
}
signed main()
{
    IOS;
    cin>>n>>m;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    ST(n);
    while(m--)
    {
        int l,r;cin>>l>>r;
        cout<<sch(l,r)<<endl;
    }
    return 0;
}
钟钟终
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
compat-libstdc++-33-3.2.3-72.el7.x86_64/libaio-devel-0.3.105
08-05
标题中的"compat-libstdc++-33-3.2.3-72.el7.x86_64"是一个重要的软件包，它在Linux环境中主要用于兼容性支持。这个包包含的是GNU标准C++库的一个特定版本（3.2.3），它允许运行那些依赖于旧版libstdc++的程序。在...
Java版计算器支持四则混合运算.+-*/()
06-08
这个过程被称为后缀表达式（也叫逆波兰表示法）转换，它简化了运算符优先级的处理。 最后，结果输出阶段从栈中取出最终结果并显示。在Java中，可以使用`System.out.println()`打印结果。 这个Java版计算器的实现还...
POJ2217 Secretary、POJ2774 Long Long Message - 最长公共子串 - 字符串哈希+二分/后缀数组
太阳升起
08-03 223
Secretary、Long Long Message     题目大概意思： 给出两个字符串 S,TS,TS,T ，求 SSS 与 TTT 的最大公共子串的长度。其中 S,TS,TS,T 的长度均不超过 10510^5105 .   样例输入 yeshowmuchiloveyoumydearmotherreallyicannotbelieveit yeaphowmu...
luoguP5108 仰望半月的夜空 [官方？]题解 后缀数组 / 后缀树 / 后缀自动机 + 线段树 / st表 + 二分...
weixin_30622181的博客
12-16 67
仰望半月的夜空 题解 可以的话，支持一下原作吧... 这道题数据很弱..... 因此各种乱搞估计都是能过的.... 算法一 暴力长度然后判断判断，复杂度\(O(n^3)\) 期望得分15分 算法二 通过二分+\(hash\)或者等等来优化字典序比较，复杂度\(O(n^2 \log n)\)，可能要松一下 期望得分30分 ps:好吧有55分... 算法三 我们考虑字符集非常小的情况 我们猜想...
超硬核十万字！全网最全 数据结构 代码，随便秒杀老师/面试官，我说的
热门推荐
hebtu666
04-20 11万+
本文代码实现基本按照《数据结构》课本目录顺序，外加大量的复杂算法实现，一篇文章足够。能换你一个收藏了吧？
Codeforces 1264C Beautiful Mirrors with queries(线段树+期望)
sigh_的博客
12-08 286
Beautiful Mirrors with queries time limit per test：2 seconds memory limit per test：256 megabytes Problem Description Creatnx has n mirrors, numbered from 1 to n. Every day, Creatnx asks exactly one mi...
CY7C68013-PVC+XC9572 开发板protel设计硬件原理图+PCB文件.zip
02-04
PVC+XC9572 FPGA 开发板protel设计硬件原理图+PCB文件计，采用2层板设计，板子大小为165x147mm，双面布局布线，主要器件为XC9572-7PC44C，CY7C68013-PVC，CY7C1049BNV33，LM150等.Protel 99se 设计的DDB后缀项目工程...
YOLOv7源代码+权重文件
07-14
压缩包中的文件“yolov7-e6e.pt”、“yolov7-d6.pt”、“yolov7-e6.pt”、“yolov7x.pt”、“yolov7-w6.pt”和“yolov7.pt”是YOLOv7的预训练权重文件。这些文件是模型训练后的结果，包含了模型在大量数据上学习到的...
2014-许智磊-后缀数组1
08-03
后缀数组是一种在字符串处理中极其重要的数据结构，由许智磊在IOI2004国家集训队论文中介绍。它是一个一维数组，包含字符串的所有后缀按照字典顺序排序后的起始索引。后缀数组的构建是通过特定算法实现的，如O(nlogn...
算法-KMP字符串匹配
xx_xb的博客
07-23 198
next[i] 就是使子串 s[0…i] 有最长相等前后缀的前缀的最后一位的下标。总体来说解next数组和模板串匹配的过程很相似，触类旁通。
AtCoder Beginner Contest 363
2301_80105412的博客
07-21 686
AtCoder Beginner Contest 363 ABCDE
【Hot100】LeetCode—152. 乘积最大子数组
weixin_44382896的博客
07-23 252
【代码】【Hot100】LeetCode—152. 乘积最大子数组。
算法力扣刷题记录 五十八【701.二叉搜索树中的插入操作】
danaaaa的博客
07-24 839
二叉搜索树向左或向右指明方向，不用遍历整个树。
数据结构与算法--【数组2】力扣练习 || 双指针 / 移除元素 / 数组排序
qq_44578108的博客
07-25 567
简单删除元素 / 指定某元素位置，其余元素位置相对不变双指针指向数组头。快指针扫描数组并判断，慢指针收集变化后的数组元素。不创建新数组，只覆盖原数组。对数组元素排序一个指针指开头，一个指针指结尾。判断两指针指向值大小，循环结束条件first <= last；for循环，while循环均可，本人更习惯while循环。使用malloc为新数组分配空间，注意，创建了新数组。
河南萌新联赛2024第（二）场：南阳理工学院
a13997460860的博客
07-25 830
这场比赛里面出到了很多的题，都是需要很多的时间去写，主要考察到了代码的实现能力，与codeforces里面的题有些差距。codeforces打多了还是有点不习惯的。题挺好的，尤其是H题感觉很考代码实现能力与STL的使用熟练度。多写写这样的题还是对我们很有好处的。只不过前三题都太水了，没有一点难度。
算法力扣刷题记录 五十六【501.二叉搜索树中的众数】
danaaaa的博客
07-21 934
求众数方法
代码随想录算法训练营day8 | 344.反转字符串、541.反转字符串 II、卡码网：54.替换数字
最新发布
m0_55763458的博客
07-25 487
代码随想录算法训练营day8 | 344.反转字符串、541.反转字符串 II、卡码网：54.替换数字
三维点云中常用的表面重建算法
Argulo的博客
07-22 1257
表面重建算法起源于计算机视觉和计算几何学领域。早期的研究集中在从二维图像中提取三维信息。然而，随着三维扫描技术的进步，越来越多的研究转向了如何从点云数据中重建表面。三维点云表面重建技术的发展，始于20世纪90年代，主要推动力是激光扫描和结构光扫描的广泛应用。
a-b*c/(3+6)#后缀式
05-23
首先按照运算优先级将中缀表达式转换为后缀表达式： a - b * c / (3 + 6) => a b c * 3 6 + / - 然后依次按照后缀表达式进行计算： 1. 将 a 压入栈中 2. 将 b 压入栈中 3. 将 c 压入栈中 4. 弹出栈顶的两个元素 c 和 b，计算 c * b 的值得到 18，并将结果压入栈中 5. 将数字 3 压入栈中 6. 将数字 6 压入栈中 7. 弹出栈顶的两个元素 6 和 3，计算 6 + 3 的值得到 9，并将结果压入栈中 8. 弹出栈顶的两个元素 18 和 9，计算 18 / 9 的值得到 2，并将结果压入栈中 9. 弹出栈顶的两个元素 a 和 2，计算 a - 2 的值得到最终的结果。 因此，后缀表达式 a-b*c/(3+6)# 的计算结果为 a - 2。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值