CSUST2020集训队选拔赛题解

最新推荐文章于 2021-03-13 19:43:22 发布
最新推荐文章于 2021-03-13 19:43:22 发布
阅读量409 收藏 2
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hddddh/article/details/105834868
版权

E.恶心心的题

题意：
给一个序列 ai，q次询问，求每次LCM(al…ar,x)的值，对p取模。
思路：

  1. 先对每个数都唯一分解吧，考虑一下怎么求多个数的 lcm;举个例子
    2 ^ 3 * 3 ^ 1* 5 ^ 7
    2 ^ 2 * 3 ^ 2 * 5 ^ 3
    2 ^ 1 *3 ^ 3 * 5 ^ 2
    这三个数的 lcm 就是 2^3 * 3^3 *5^7 ,也就是各个质数指数的最大值了。因为总共只有60多个质数，可以用 60棵线段树维护最大值（还可以用 rmq）。

线段树做法：

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+20;
int isprime[301];
int pri[70];
int a[maxn][65],sum[maxn<<2][65],k,s[65];
ll poww(int a,int b){//快速幂
       long long ans=1;
       while(b>0){
            if(b&1) ans=ans*a;
            a=a*a;
            b>>=1;
       }
      return ans;
}
void sieve(){
    for(int i=2;i<=300;i++){
        if(isprime[i]==0){
            for(int j=i*i;j<=300;j+=i){
                isprime[j]=1;
            }
            pri[++k]=i;
        }
    }
}
void pushup(int rt,int x){
     sum[rt][x]=max(sum[rt<<1][x],sum[rt<<1|1][x]);
}
void build(int l,int r,int rt,int x){
     if(l==r){
        sum[rt][x]=a[l][x];
        return ;
     }
     int m=(l+r)>>1;
     build(l,m,rt<<1,x);
     build(m+1,r,rt<<1|1,x);
     pushup(rt,x);
}
void query(int L,int R,int l,int r,int rt){
    if(L<=l&&R>=r){
        for(int i=1;i<=62;i++){
            s[i]=max(s[i],sum[rt][i]);
        }
        return ;
    }
    int m=(l+r)>>1;
    if(L<=m) query(L,R,l,m,rt<<1);
    if(R>m)  query(L,R,m+1,r,rt<<1|1);

}
int main (){
    int p,q,x,n;
    cin>>n>>p>>q;
    sieve();
    for(int i=1;i<=n;i++){
        scanf("%d",&x);
        for(int j=1;j<=62&&x!=1;j++){
            int cnt=0;
                while(x%pri[j]==0){
                     cnt++;
                    x/=pri[j];
                }
                a[i][j]=cnt;
        }
    }
   for(int i=1;i<=62;i++){
      build(1,n,1,i);
   }
    while(q--){
        int l,r,x,ans;
        ll anss=1;
        scanf("%d%d%d",&l,&r,&x);
        memset(s,0,sizeof(s));
        query(l,r,1,n,1);
        for(int j=1;j<=62;j++){
            int cnt=0;
                while(x%pri[j]==0){
                    cnt++;
                    x/=pri[j];
                }
         ans=max(s[j],cnt);
         anss=anss*poww(pri[j],ans)%p;
        }
        printf ("%lld\n",anss*x%p);
    }
    return 0;
}

st表：

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+20;
int isprime[300];
int pri[70],n,ans[70],k;
short int a[65][maxn],dpma[65][maxn][17];
const int mod=1e9+7;
long long poww(long long a,long long b){//快速幂
       long long ans=1;
       while(b>0){
            if(b&1) ans=ans*a;
            a=a*a;
            b>>=1;
       }
      return ans;
}
void sieve(){
    for(int i=2;i<=300;i++){
        if(isprime[i]==0){
            for(int j=i*i;j<=300;j+=i){
                isprime[j]=1;
            }
            pri[++k]=i;
        }
    }
}
void init(){
     for(int i=1;i<=62;i++){
         for(int j=1;j<=n;j++){
            dpma[i][j][0]=a[i][j];
         }
     }
     for(int x=1;x<=62;x++){
        for(int j=1;(1<<j)<=n;j++){
           for(int i=1;i+(1<<j)-1<=n;i++){
            dpma[x][i][j]=max(dpma[x][i][j-1],dpma[x][i+(1<<(j-1))][j-1]);
           }
        }
     }
}
int qa(int x ,int l, int r){
    int k=log((r-l+1)*1.0)/log(2.0);
    return max(dpma[x][l][k],dpma[x][r-(1<<k)+1][k]);
}
int main (){
    int p,q,x;
    cin>>n>>p>>q;
    sieve();
    for(int i=1;i<=n;i++){
        scanf("%d",&x);
        for(int j=1;j<=62&&x!=1;j++){
            if(x%pri[j]==0){
                while(x%pri[j]==0){
                    a[j][i]++;
                    x/=pri[j];
                }
            }
        }
    }
    init();
    while(q--){
        int l,r,x,ans;
        ll anss=1;
        scanf("%d%d%d",&l,&r,&x);
         for(int j=1;j<=62;j++){
            int cnt=0;
            if(x%pri[j]==0){
                while(x%pri[j]==0){
                    cnt++;
                    x/=pri[j];
                }
            }
            int t=qa(j,l,r);
         ans=max(t,cnt);
         anss=anss*poww(pri[j],ans)%p;
        }
        printf ("%lld\n",anss*x%p);
    }
    return 0;
}

需要注意的：

  1. 线段树查询的时候，不能像st表一样查询60次，这样会t,仔细想一下线段树查询的实质，其实就是找出所有需要比较的节点的下标，如果查询60次，就会重复这个步骤。所以我们只要查询一次，每到达一个区间就维护一下60个数的最大值，用一个数组记录一下。
  2. st表会卡内存，所以用了 short int.
  3. 最后X的范围是 1e9 ，还能分解出300以外的质数，所以最后还要乘 x 分解出300以内的质数后的值。

ps:虽然hdw聚聚每次的题坑点都挺多的，但写完收获也很大，吹爆hdwdl.

I.摸鱼的tomjobs2

题意：
给n个数，每一个连续区间的按位与，作为一个交叉值，求有多少个不同的交叉值。
思路：

  1. 考虑以每一个数为左端点，考虑每一个二进制位，要想使得这个二进制位发生改变，就要找到右边第一个不为1的二进制位作为右端点，然后记录这一段的值。这个方法必须先用 st表预处理 区间值，至于为什么能用 st 表，聚聚告诉我这是个可覆盖的区间，豁然开朗。然后找右边第一个不为1的位置，虽然也可以预处理出来，但貌似不太会，，我用了 前缀和+二分。
  2. 考虑dp的思路，每次都往右扩展一位，把每次能到达的状态塞进 set 然后不断地转移。但set中的元素为什么不会超过 62个 本菜鸡始终没想明白。

法1：

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+7;
map<ll,int >mp;
inline ll read(){
    char c=getchar();
    ll f=1,x=0;
    while(c<'0'||c>'9'){
        if(c=='-')
            f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        x=(x<<1)+(x<<3)+(c^'0');
        c=getchar();
    }
    return x*f;
}
ll a[maxn],dp[maxn][32];
int n,pre[64][maxn];
void init(){
    for(int i=1;i<=n;i++) dp[i][0]=a[i];
    for(int j=1;(1<<j)<=n;j++){
        for(int i=1;i+(1<<j)-1<=n;i++){
            dp[i][j]=dp[i][j-1]&dp[i+(1<<(j-1))][j-1];
        }
    }
}
ll qi(int l,int r){
    int k=log((r-l+1)*1.0)/log(2.0);
    return dp[l][k]&dp[r - (1 << k) + 1][k];
}
int main (){
    int ans=0;
    n=read();
    for(int i=1; i<=n; i++){
        a[i]=read();
        if(mp[a[i]]==0){
            mp[a[i]]=1;
            ans++;
        }
    }
    init();
    for(int i=0;i<61;i++){
        for(int j=1;j<=n;j++){
            pre[i][j]=pre[i][j-1]+(((a[j]>>i)&1)==0);
        }
    }
    for(int i=1; i<=n; i++){
        for(int j=0;j<61;j++){
            int l=i,r=n,mid,t=i-1;
            while(l<=r){
                mid=(l+r)/2;
                if(pre[j][mid]-pre[j][i-1]==0) t=mid,l=mid+1;
                else r=mid-1;
            }
            if(t<n) t++;
            ll x=qi(i,t);
            if(mp[x]==0) mp[x]=1,ans++;
        }
    }
    printf ("%d\n",ans);

}

法2:

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<set>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+7;

string s1,s2;
set < ll > dp[maxn];
ll a[maxn];
map<ll ,int >mp;
int main (){
    int n,ans=0;
    cin>>n;
    for(int i=1;i<=n;i++){
        scanf("%lld",&a[i]);
        dp[i].insert(a[i]);
        if(!mp[a[i]]){
            mp[a[i]]=1;
            ans++;
        }
    }
    for(int i=1;i<n;i++){
        set<ll>::iterator it;
        for(it=dp[i].begin();it!=dp[i].end();it++){
            dp[i+1].insert(*it&a[i+1]);
            if(!mp[*it&a[i+1]]){
                mp[*it&a[i+1]]=1;
                ans++;
            }
        }
    }
    printf ("%d\n",ans);
}

F.打扑克牌

题意：
给你一个长度为n的数字串，你可以任意打乱顺序，求有多少个不同的数字串可以被m整除。（n<=15,m<=50）

思路：

  1. 一开始想的裸的 dfs（状压没入门，可怜），交了一发 t 了，仔细想想 15! 不T 就怪了。然后又想到数字只能是 0-9，只要记录 0-9数字的数量，再进行 dfs 会有一些优化，比如 11111，就只会搜索一次，而原来要搜索 5! 次相同状态。然而还是会 t .
  2. 还是别挣扎了，还是需要记录状态呀，不然会 t 傻的。
  3. 状压dp 就是把状态用二进制压缩作为dp的状态。例如

大概是这样的

dp[3] 011  //代表你已经选了 第一个数和第二个数。
dp[2] 010  //代表选了第二个数
dp[1] 001  //代表选了第一个数

dp[3] 可以由 dp[2]和dp[1]转移
     dp[2|(1<<0)]+=dp[2];
     dp[1|(1<<1)]+=dp[1];
     1.由 dp[2]转移相当于第一次选了第二个数第二次选第一个数。
     2.由 dp[1]转移相当于第一次选了第一个数第二次选第二个数。

这个题需要再加一个模数的状态。
大概就是这样了

for(int i=0; i<(1<<n); i++) {
        for(int j=0; j<n; j++){
            if((i>>j)&1) continue;  //表示第j个数已经选过了
            for(int k=0; k<m; k++){
                dp[i|(1<<j)][(k*10+a[j])%m]+=dp[i][k];
            }
        }
    }

最后的答案就是 dp[(1<<n)-1][0],因为会出现重复的数，所以最后还要除以每个数重复次数的阶乘。

#include<iostream>
#include<cstdio>
#include<map>
#include<math.h>
#include<queue>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn=1e5+7;

string s;
ll m,a[20],dp[1<<16][60],n,cnt[20];
ll num[20];
int main ()
{
    cin>>s>>m;
    n=s.size(),num[0]=1,dp[0][0]=1;
    for(int i=0; i<n; i++)
        a[i]=s[i]-'0',cnt[a[i]]++,num[i+1]=num[i]*(i+1);
    for(int i=0; i<(1<<n); i++)
    {
        for(int j=0; j<n; j++)
        {
            if((i>>j)&1)
                continue;
            for(int k=0; k<m; k++)
            {
                dp[i|(1<<j)][(k*10+a[j])%m]+=dp[i][k];
            }
        }
    }
    ll ans=dp[(1<<n)-1][0];
    for(int i=0;i<=9;i++){
        ans/=num[cnt[i]];
    }
    printf ("%lld\n",ans);

}

最后放上一张偷来的图
在这里插入图片描述
也算是第一次写 状压dp 吧

stduy_ing
关注
ACM 训练大纲（CSUST_ACM）
bo-jwolf的专栏
04-29 2449
ACM 训练大纲 Changsha University of Science & Technology July 31, 2012 1 推荐题库 • http://ace.delos.com/usaco/ 美国的OI 题库，如果是刚入门的新手，可以尝试先把它刷通，能够学到几乎全部的基础算法极其优化， 全部的题解及标程还有题目翻译可以baidu 一个叫N
2015 CSUST校赛 - 超级快速幂【费马小定理】+【快速幂取模】
Fighting!
05-02 2425
 超级快速幂 Time Limit: 3000/1000 MS(Java/Others)   Memory Limit:65536/65536 K (Java/Others) Description   请计算：    a^(b^c)mod(1e9+7) Input   多组数据，第一行为一个整数T，代表数据组数，接下来每行
CSUST 集训队选拔赛题解
airen6935的博客
03-06 219
选拔赛的题解，~~~ 题目链接：请点击 A题 素数筛 + 线段树（树状数组） 先用素数筛打表，然后线段树更新，遍历求出值，O（1）查询即可 AC代码： /*num数组 是把记录 数是否存在 存在即为1。 总共有N个数，如何判断第i+1个数到最后一个 数之间有多少个数小于第i个数呢？不妨假设 有一个区间 [1,N]，只需要判断区间[i+1,N]之 间有多少...
CF - 658 -- B.Sequential Nim【简单博弈论】
qq_45249273的博客
07-23 394
Sequential Nim 题意： 顺序给出n堆石头，有两名玩家(First and Second) 首先是First先拿第一堆石头，它可以拿至少一个，如果谁不能拿了，则另一方获胜。 思路： 简单博弈论：分析子问题 2 3 1，这肯定是第一个拿第一堆的人(First)获胜，再来分析 5 2 2 2，根据前面分析的2 3 1，2 2 2 同理也是当前先拿的人获胜，所以第一个拿第一堆的人可以计划他也是第一个拿第二堆的人，那么它可以先拿到第一堆石头的a[1]-1个，这样他也是第一个拿第二堆的，由此分析，谁先
HUT_CSUST 新生友谊赛题解CSUST版本
weixin_30678821的博客
12-22 104
所有的题目百度都有,问问题前先百度,这是黄金定律 所有题都给个思路:具体代码都可百度 A: 把素数筛出来,先构造原矩形,就是一个一个往下推,然后构造被挖掉合数的矩形 遍历原矩形找到起点和终点,然后后构造的矩形中的素数就相当于障碍物,整个问题就化成一个迷宫问题了 用BFS一搜就能得出最短距离 B: 具体证明: http://www.doc88.com/p-792551521...
CSUST2021集训队选拔赛题解
the best 的博客
03-13 521
本场难度： 简单：A,B,E,H,K 中等：C,D,F,G,I 签到题又双叒叕来了 线段树模板题，在建树时，把奇数位的数设为正数，偶数位置的数设为负数，在查询区间和的时候，判断lll 位置是正还是负，如果为负，只需要把区间和乘上-1即可。 代码： #include<bits/stdc++.h> #define fi first #define se second #define debug cout<<"I AM HERE"<<endl; using namespace
get-csust-info-helper
06-02
"get-csust-info-helper" 是一个专门为长沙理工大学学生设计的工具，主要功能是帮助学生查询校园信息，包括课程表和图书馆的借阅信息。这个工具最初由用户yinbro开发，并在Gitcafe上进行了分享。后来，该项目被其他...
CSUST_:长沙理工大学——校园客户端
05-26
"CSUST_"，即长沙理工大学的校园客户端，旨在提供一个便捷、全面的服务平台，帮助师生更好地融入校园生活。这款客户端是用Java编程语言开发的，充分体现了Java在跨平台应用开发中的优势。 首先，我们来探讨一下Java...
2013ACM多校联合(1)_CSUST(标程+数据)第一次修改
03-17
标题 "2013ACM多校联合(1)_CSUST(标程+数据)第一次修改" 暗示这是一个关于ACM（国际大学生程序设计竞赛，International Collegiate Programming Contest）的比赛资料集，包含了2013年3月17日某场赛事的题目解析、...
spark_3_2_0-master-3.2.3-1.el7.noarch.rpm
09-16
Ambari+Bigtop 一站式编译和部署解决方案 https://gitee.com/tt-bigdata/ambari-env
浙大城市学院在河南2021-2024各专业最低录取分数及位次表.pdf
09-16
全国各大学在河北2021-2024年各专业最低录取分数及录取位次数据，高考志愿必备参考数据
第4周玩转案例分析.pdf
09-16
第4周玩转案例分析
基于MATLAB的教室人数统计系统源代码+使用说明，带有丰富的人机交互GUI界面
09-16
基于MATLAB的教室人数统计，带有丰富的人机交互GUI界面。实现教室人数的计数统计，进而统计出勤率，上课认真听讲的比例。原先预设总人数，未出勤或者书本遮住脸部、趴着玩手机等必然无法采集到人脸，从而对应数量缺失，得出出勤率或听讲率。该课题采用肤色原理进行人脸定位和人脸分割。本课题中，为了界面GUI的美观，刻意将分割出的人脸单独显示在GUI对应的axes里
java-ssm+jsp药品销售网站系统实现源码(项目源码-说明文档)
09-16
药品销售系统的主要使用者分为管理员、用户和医生，实现功能包括管理员；管理员个人中心、用户管理、医生管理、药品信息管理、药品分类管理、订单配送管理、系统管理、订单管理，用户；个人中心、订单配送管理、订单管理，医生：个人中心、药品信息管理、药品分类管理，前台首页：首页、药品信息、我的、跳转到后台、购物车、客服等功能。 项目关键技术 开发工具：IDEA 、Eclipse 编程语言: Java 数据库: MySQL5.7+ 后端技术：ssm 前端技术：jsp 关键技术：jsp、spring、ssm、MYSQL、MAVEN 数据库工具：Navicat、SQLyog
AI健身体能测试之基于paddlehub实现引体向上计数个数统计源码+模型+视频例子+视频结果文件.zip
最新发布
09-16
使用环境： paddlepaddle==2.4.2 paddlehub==2.4.0 理论上安装最新paddle环境也行。 实现的细节参考博文：https://blog.csdn.net/FL1623863129/article/details/142307043
CPA《财务成本管理》刘正兵 专题班 资本成本 债务成本的估计+加权平均资本成本.pdf
09-16
CPA《财务成本管理》刘正兵 专题班 资本成本 债务成本的估计+加权平均资本成本.pdf
CPA《财务成本管理》刘正兵 专题班 风险与报酬 资本资产定价模型.pdf
09-16
CPA《财务成本管理》刘正兵 专题班 风险与报酬 资本资产定价模型.pdf
智能优化特征选择-基于FOA果蝇优化的特征选择分类算法（matlab程序）
09-16
智能优化特征选择|基于FOA果蝇优化的特征选择分类算法（matlab程序） 1.程序功能已完成调试，用户可以通过一键操作生成图形和评价指标。 2.数据输入以Excel格式保存，只需更换文件，即可运行以获得个人化的实验结果。 3.代码中包含详细注释，具有较强的可读性，特别适合初学者和新手。 4.在实际数据集上的效果可能较差，需要对模型参数进行微调。 CSDN：机器不会学习CL 智能优化特征选择|基于FOA果蝇优化的特征选择分类算法（matlab程序） 1.程序功能已完成调试，用户可以通过一键操作生成图形和评价指标。 2.数据输入以Excel格式保存，只需更换文件，即可运行以获得个人化的实验结果。 3.代码中包含详细注释，具有较强的可读性，特别适合初学者和新手。 4.在实际数据集上的效果可能较差，需要对模型参数进行微调。 CSDN：机器不会学习CL
gee在首页该如何导航至代码块
10-06
1. 打开浏览器，并输入Gee的官方网站地址：https://gee.csust.edu.cn/。 2. 在首页上寻找代码块导航的相关按钮或链接。一般来说，这些导航选项会以明显的方式展示在首页的某个位置，比如菜单栏、导航栏或页面的底部...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值