2021“MINIEYE杯”中国大学生算法设计超级联赛(8)(1002,1004,1006,1009)

前言

依旧是白嫖账号,只打了一些题/kk


正题


1002 Buying Snacks

题目大意

nnn个物品,每个可以买一次也可以不买,如果买需要选择1/21/21/2块钱的,然后也可以相邻两个一起买并且减少一块的花销,求恰好用掉mmm块钱的方案。

1≤n≤109,1≤m≤20000,1≤T≤51\leq n\leq 10^9,1\leq m\leq 20000,1\leq T\leq 51n109,1m20000,1T5

解题思路

fi,jf_{i,j}fi,j表示iii个物品花jjj块钱的方案那么有
fi,j=fi−1,j+fi−1,j−1+fi−1,j−2+fi−2,j−1+2×fi−2,j−2+fi−2,j−3f_{i,j}=f_{i-1,j}+f_{i-1,j-1}+f_{i-1,j-2}+f_{i-2,j-1}+2\times f_{i-2,j-2}+f_{i-2,j-3}fi,j=fi1,j+fi1,j1+fi1,j2+fi2,j1+2×fi2,j2+fi2,j3
然后化成生成函数就是
Fi(x)=(1+x+x2)Fi−1(x)+(x+2x2+x3)Fi−2(x)F_i(x)=(1+x+x^2)F_{i-1}(x)+(x+2x^2+x^3)F_{i-2}(x)Fi(x)=(1+x+x2)Fi1(x)+(x+2x2+x3)Fi2(x)
和之前一道倍增FFTFFTFFT很像,然后考虑分割位置填两个就是F1(x)=1+x+x2F_{1}(x)=1+x+x^2F1(x)=1+x+x2,然后
Fa+b(x)=Fa(x)Fb(x)+(x+2x2+x3)Fa−1(x)Fb−1(x)F_{a+b}(x)=F_{a}(x)F_b(x)+(x+2x^2+x^3)F_{a-1}(x)F_{b-1}(x)Fa+b(x)=Fa(x)Fb(x)+(x+2x2+x3)Fa1(x)Fb1(x)
上倍增FFTFFTFFT就好了。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=1<<16,P=998244353;
ll T,n,k,m,r[N],f[3][N],t[3][N],g[2][N];
void fm(ll &x){x+=x>>31&P;}
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
void NTT(ll *f,ll op){for(ll i=0;i<n;i++)if(i<r[i])swap(f[i],f[r[i]]);for(ll p=2;p<=n;p<<=1){ll tmp=power(3,(P-1)/p),len=p>>1;if(op==-1)tmp=power(tmp,P-2);for(ll k=0;k<n;k+=p){ll buf=1;for(ll i=k,tt;i<(k|len);i++){tt=1ll*buf*f[i|len]%P;fm(f[i|len]=f[i]-tt);fm(f[i]=f[i]+tt-P);buf=1ll*buf*tmp%P;}}}if(op==-1){ll invn=power(n,P-2);for(ll i=0;i<n;i++)f[i]=1ll*f[i]*invn%P;}return;
}
void print(ll x)
{if(x>9)print(x/10);putchar(x%10+'0');return;}
signed main()
{scanf("%lld",&T);while(T--){memset(f,0,sizeof(f));memset(g,0,sizeof(g));memset(t,0,sizeof(t)); scanf("%lld%lld",&m,&k);k++;n=1;while(n<(k*2))n<<=1;for(ll i=0;i<n;i++)r[i]=(r[i>>1]>>1)|((i&1)?(n>>1):0);f[0][0]=f[0][1]=f[0][2]=f[1][0]=g[0][0]=1;for(ll d=1;d<=m;d<<=1){if(m&d){for(ll j=0;j<3;j++){for(ll i=0;i<n;i++)t[j][i]=(i<k)?f[j][i]:0;NTT(t[j],1);}NTT(g[0],1);NTT(g[1],1);for(ll i=0;i<n;i++){ll b0=g[0][i],b1=g[1][i];g[0][i]=1ll*b0*t[0][i]%P;g[1][i]=1ll*b0*t[1][i]%P;t[0][i]=1ll*t[1][i]*b1%P;t[1][i]=1ll*t[2][i]*b1%P;}NTT(g[0],-1);NTT(g[1],-1);NTT(t[0],-1);NTT(t[1],-1);for(ll i=0;i<k;i++){(g[0][i+1]+=t[0][i])%=P;(g[0][i+2]+=2ll*t[0][i])%=P;(g[0][i+3]+=t[0][i])%=P;(g[1][i+1]+=t[1][i])%=P;(g[1][i+2]+=2ll*t[1][i])%=P;(g[1][i+3]+=t[1][i])%=P;}for(ll i=k;i<n;i++)g[0][i]=g[1][i]=0;}if(d*2>m)break;for(ll j=0;j<3;j++){for(ll i=0;i<n;i++)t[j][i]=(i<k)?f[j][i]:0;NTT(t[j],1);}for(ll i=0;i<n;i++){f[0][i]=1ll*t[0][i]*t[0][i]%P;f[1][i]=1ll*t[1][i]*t[0][i]%P;f[2][i]=1ll*t[1][i]*t[1][i]%P;t[0][i]=1ll*t[1][i]*t[1][i]%P;t[1][i]=1ll*t[1][i]*t[2][i]%P;t[2][i]=1ll*t[2][i]*t[2][i]%P;}for(ll j=0;j<3;j++)NTT(f[j],-1),NTT(t[j],-1);for(ll i=0;i<k-1;i++){(f[0][i+1]+=t[0][i])%=P;(f[0][i+2]+=2ll*t[0][i])%=P;(f[0][i+3]+=t[0][i])%=P;(f[1][i+1]+=t[1][i])%=P;(f[1][i+2]+=2ll*t[1][i])%=P,(f[1][i+3]+=t[1][i])%=P;(f[2][i+1]+=t[2][i])%=P;(f[2][i+2]+=2ll*t[2][i])%=P;(f[2][i+3]+=t[2][i])%=P;}for(ll i=k;i<n;i++)f[0][i]=f[1][i]=f[2][i]=0;}for(ll i=1;i<k;i++)print(g[0][i]),putchar(' ');putchar('\n');}return 0;
}

1004 Counting Stars

题目大意

nnn个数字要求支持操作

  1. 将一个区间所有数字的最低位减去
  2. 将一个区间非零数字的最高位提高一位

每次操作玩输出操作区间的和

1≤T≤5,1≤n,q≤105,1≤ai≤1091\leq T\leq 5,1\leq n,q\leq 10^5,1\leq a_i\leq 10^91T5,1n,q105,1ai109

解题思路

每个数字操作一最多执行log⁡ai\log a_ilogai次所以这个操作暴力做就好了

然后第二个操作我们把每个数字的最高位分离出来就变成了区间乘二。

代码是lemondinosaur\text{lemondinosaur}lemondinosaur爷写的


1006 GCD Game

题目大意

nnn个数字每次操作的人可以选择一个数字xxx然后找一个1≤k<x1\leq k<x1k<xxxx变为gcd(x,k)gcd(x,k)gcd(x,k)

不能操作者输,求是否先手必胜
1≤T≤100,∑n≤106,1≤ai≤1071\leq T\leq 100,\sum n\leq 10^6,1\leq a_i\leq 10^71T100,n106,1ai107

解题思路

其实就是去掉任意个质因子,所以用线性筛筛出每个数字的最小质因子,然后递推出每个数字的质因子数就好了。
时间复杂度O(ai+∑n)O(a_i+\sum n)O(ai+n)

code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e7+10;
int T,n,cnt,pri[N/10],v[N],c[N];
void Prime(){v[1]=1;for(int i=2;i<N;i++){if(!v[i])pri[++cnt]=i,v[i]=i;for(int j=1;j<=cnt&&i*pri[j]<N;j++){v[i*pri[j]]=pri[j];if(i%pri[j]==0)break;}}return;
}
int main()
{Prime();for(int i=2;i<N;i++)c[i]=c[i/v[i]]+1;scanf("%d",&T);while(T--){scanf("%d",&n);int ans=0;for(int i=1,x;i<=n;i++)scanf("%d",&x),ans=ans^c[x];if(ans)puts("Alice");else puts("Bob");}return 0;
}

1009 Singing Superstar

题目大意

给一个字符串SSSnnn次询问字符串aia_iaiSSS中出现的最大不重次数。

1≤T≤5,1≤n≤105,1≤∣ai∣≤30,∑∣S∣,∑∣ai∣≤4×1051\leq T\leq 5,1\leq n\leq 10^5,1\leq |a_i|\leq 30,\sum |S|,\sum|a_i|\leq 4\times 10^51T5,1n105,1ai30,S,ai4×105

解题思路

根据贪心的思想是能够匹配尽量匹配
把所有SSS长度不超过303030的子串建一棵TrieTrieTrie,然后记下每个子串的最后出现位置然后统计ansansans就好了。

code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ull unsigned long long
using namespace std;
const int N=1e5+10;
int n,T,l,cnt,ans[N*30],t[N*30][26],last[N*30];
char s[N],st[N];
int main()
{scanf("%d",&T);while(T--){scanf("%s",s+1);l=strlen(s+1);cnt=1;last[1]=ans[1]=0;memset(t[1],0,sizeof(t[1]));for(int i=1;i<=l;i++){int x=1;for(int j=1;j<=30;j++){if(i-j+1<1)break;int c=s[i-j+1]-'a';if(!t[x][c]){cnt++;memset(t[cnt],0,sizeof(t[cnt]));t[x][c]=cnt;last[cnt]=ans[cnt]=0;}x=t[x][c];if(i-last[x]>=j)last[x]=i,ans[x]++;}}scanf("%d",&n);while(n--){scanf("%s",st+1);int tl=strlen(st+1),x=1;for(int i=tl;i>=1;i--)x=t[x][st[i]-'a'];printf("%d\n",ans[x]);}}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/318498.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

分羊(区间dp:分治与决策单调性优化)

文章目录题目描述样例输入样例输出解析代码题目描述 给你一个数列&#xff0c;分成m段&#xff0c;每段的价值为相同数对的对数 求最小价值和 样例输入 10 2 1 2 1 2 1 2 1 2 1 2 样例输出 8 解析 容易想到区间dp 定义dp[i][j]:把前i个数分成j段的最小价值 那么枚举最后一段…

网络流专题(最大流与费用流)(一)

流量网络 • 想要将一些水从S运到T&#xff0c;必须经过一些水站&#xff0c;链接水站的是管道&#xff0c;每条管道都有它的最大能容纳的水量&#xff0c;求最多S到T能流多少流量。 基本概念 • 这是一个典型的网络流模型。我们先了解网络流的有关定义和概念。 • 若有向图G(…

eShopOnContainers 看微服务 ①:总体概览

一、简介eShopOnContainers是一个简化版的基于.NET Core和Docker等技术开发的面向微服务架构的参考应用。该参考应用是一个简化版的在线商城/电子商务微服务参考示例应用。其包含基于浏览器的Web应用、基于Xamarin的Android、IOS、Windows/UWP 移动应用&#xff0c;以及服务端应…

【模板】文艺平衡树

题目链接 题目描述 您需要写一种数据结构&#xff08;可参考题目标题&#xff09;&#xff0c;来维护一个有序数列。 其中需要提供以下操作&#xff1a;翻转一个区间&#xff0c;例如原有序序列是 543215\ 4\ 3\ 2\ 15 4 3 2 1&#xff0c;翻转区间是 [2,4][2,4][2,4] 的话&a…

P4022-[CTSC2012]熟悉的文章【广义SAM,dp,单调队列】

正题 题目链接:https://www.luogu.com.cn/problem/P4022 题目大意 给出mmm个模板串。 然后nnn次询问给出一个串SSS要求找到一个最大的LLL使得能够将SSS超过90%90\%90%的部分拿出来分后每个串都是某个模板串的子串且长度不小于LLL。 所有输入文件长度不超过 110000011000001…

猜数(二分、线段树)

题目描述 对于一个长度为n的数列给出m个描述 每一个描述给出一个区间[a,b]的最小值的x 求从第几个描述开始矛盾 解析 本题关键是一个关于矛盾的充要条件&#xff1a; 如果存在一个最小值x&#xff0c;其所在的区间的交集&#xff08;就是它真正可以存在的区间&#xff09;是…

网络流专题(最大流与费用流)例题总结

文章目录NC 106056 poj1459 Power Network题目大意&#xff1a;题解&#xff1a;NC213817 [网络流24题]最小路径覆盖问题题目&#xff1a;题解&#xff1a;例2&#xff1a;NC213818 [网络流24题]魔术球问题题目&#xff1a;题解&#xff1a;方法2&#xff1a;NC 213820 [网络流…

FHQ Treap

FHQ Treap 主要通过 merge⁡\operatorname{merge}merge 和 spilt⁡\operatorname{spilt}spilt 两个核心操作和其他辅助函数来实现。 本文用 lcxlc_xlcx​ 表示 xxx 的左儿子&#xff0c;rcxrc_xrcx​ 表示 xxx 的右儿子&#xff0c;datxdat_xdatx​ 表示 xxx 的随机优先级 me…

从.Net框架Bug的提交到修复代码成功合并到.NET CoreFX主线

从发现.NET Framework中SmtpClient的Bug并拿出解决方案&#xff0c;然后给微软开发者社区提交Bug开始&#xff0c;总共耗时一个多月&#xff0c;对Bug修复的代码最终被采纳&#xff0c;现已合并到.NET Core Libraries (CoreFX)主线中。修复记录https://github.com/dotnet/coref…

P2469-[SDOI2010]星际竞速【费用流】

正题 题目链接:https://www.luogu.com.cn/problem/P2469 题目大意 给出nnn个点mmm条边的一张无向图&#xff0c;你每次可以选择一条s−>ts->ts−>t的路径花费AsA_sAs​的代价出发&#xff0c;而且只能从编号小的点到编号大的点走&#xff0c;求最小花费使得所有点被…

周期长度和(KMP)

文章目录题目描述解析问题总结代码题目描述 解析 我们可以看到 如果A是B的周期 那么B一定可以写成&#xff1a; A1A2A1 的形式 注意到&#xff1a;A1就是KMP中的公共前后缀 要使A最大&#xff0c;要使A1最短 也就是求最短公共前后缀 这怎么求呢&#xff1f; 我们注意到&#x…

计算几何基础-1

文章目录基本概念点与向量的运算精度问题线段&#xff0c;射线和直线点积&#xff1a;夹角叉积向量的极角旋转一个向量求三角形面积直线交点点到直线距离点在直线上的投影判断两条线段是否相交点与直线的位置关系点是否在直线左侧点是否在直线上点是否在线段上点与多边形的位置…

.net core i上 K8S(四).netcore程序的pod管理,重启策略与健康检查

目录1.pod管理2.重启策略3.健康检查4.进入容器正文上一章我们已经通过yaml文件将.netcore程序跑起来了&#xff0c;但还有一下细节问题可以分享给大家。1.pod管理1.1创建podkubectl create -f netcore-pod.yaml我们创建一个netcore-pod.yaml文件&#xff0c;内容如下&#xff1…

【模板】可持久化并查集

题目链接 题目描述 给定 nnn 个集合&#xff0c;第 iii 个集合内初始状态下只有一个数&#xff0c;为 iii。 有 mmm 次操作。操作分为 333 种&#xff1a; 1 a b 合并 a,ba,ba,b 所在集合&#xff1b; 2 k 回到第 kkk 次操作&#xff08;执行三种操作中的任意一种都记为一次…

P1758-[NOI2009]管道取珠【dp】

正题 题目链接:https://www.luogu.com.cn/problem/P1758 题目大意 给出一个大小为nnn和一个大小为mmm的栈&#xff0c;每次选择一个栈弹出栈顶然后记录这个字母&#xff0c;求所有弹出序列的弹出方案的二次方和。 1≤n,m≤5001\leq n,m\leq 5001≤n,m≤500 解题思路 二次方和…

洛谷P2680:运输计划(倍增、二分、树上差分)

传送门 文章目录题目描述解析问题代码题目描述 解析 求最大值的最小值 容易想到二分 然后。。。就没有然后了。。。 看了题解 学会了一个新技能&#xff1a;树上差分 &#xff08;其实学长之前好像讲过。。。&#xff09; 一般的&#xff0c;对于一条A到B的路径&#xff0c;如…

Loj#2880-「JOISC 2014 Day3」稻草人【CDQ分治,单调栈,二分】

正题 题目链接:https://loj.ac/problem/2880 题目大意 给出平面上的nnn个点&#xff0c;然后求有多少个矩形满足 左下角和右上角各有一个点矩形之间没有其他点 1≤n≤2105,1≤xi,yi≤109,1\leq n\leq 2\times 10^5,1\leq x_i,y_i\leq 10^9,1≤n≤2105,1≤xi​,yi​≤109,保证…

计算几何基础-2

文章目录直线&#xff1a;图形&#xff1a;求垂足求两圆交点直线与圆交点多边形问题判断一个点是否在任意多边形内部Pick定理凸包求点集的凸包水平法&#xff1a;增量法&#xff1a;半平面半平面交求半平面交直线&#xff1a; struct Line{point p,v;Line(){}Line(point _p.po…

eShopOnContainers 看微服务 ②:配置 启动

一、什么是dockerDocker 是一个开源项目&#xff0c;通过把应用程序打包为可移植的、自给自足的容器&#xff08;可以运行在云端或本地&#xff09;的方式&#xff0c;实现应用程序的自动化部署。使用 Docker 的时候&#xff0c;需要创建一个应用或服务&#xff0c;然后把它和它…

字符串距离(opj )(动态规划)

题目描述2988:计算字符串距离 对于两个不同的字符串&#xff0c;我们有一套操作方法来把他们变得相同&#xff0c;具体方法为&#xff1a; 修改一个字符&#xff08;如把“a”替换为“b”&#xff09;删除一个字符&#xff08;如把“traveling”变为“travelng”&#xff09;…