2017省夏令营Day7
题解:首先,我们可以得到一个规律:经过2次变换后,a和b的值都分别乘2了,所以只要用快速幂就能过啦,但是,要特判n为0的情况。
代码如下:
1 #include<cstdio> 2 #include<iostream> 3 #include<cstring> 4 #define Mod 1000000007 5 using namespace std; 6 long long a,b,n,ans1,ans2; 7 long long power(long long x){ 8 long long ret=1,tmp=2; 9 while(x){ 10 if(x&1){ret*=tmp; ret%=Mod;} 11 tmp*=tmp; tmp%=Mod; x/=2; 12 } 13 return ret; 14 } 15 int main() 16 { 17 freopen("apexis.in","r",stdin); 18 freopen("apexis.out","w",stdout); 19 scanf("%lld%lld%lld",&a,&b,&n); 20 if(n==0){ 21 printf("%lld %lld",a,b); 22 return 0; 23 } 24 if(a<b) swap(a,b); 25 long long tmp=n>>1; 26 long long t=power(tmp); 27 if(n&1) ans1=t*(a+b)%Mod,ans2=t*(a-b)%Mod; 28 else ans1=t*a%Mod,ans2=t*b%Mod; 29 printf("%lld %lld",ans1,ans2); 30 return 0; 31 }
------------------------------------------------------------华丽的分割线----------------------------------------------------------------------
题解:我们记f[i]表示区间[1,i]内素数个数,我们可以用筛法筛出数据范围内的素数并顺便求f数组,然后我们暴力枚举1-maxn中的数即可。
代码如下:
1 #include<cstdio> 2 #include<iostream> 3 #define Max 200005 4 using namespace std; 5 int n,m,ans,cnt,maxn; 6 int f[Max],num[Max],prime[Max]; 7 void shaifa(){ 8 for(int i=1;i<=Max;i++) f[i]=i; 9 for(int i=2;i<=Max;i++) 10 if(!prime[i]){ 11 f[i]--; 12 for(int j=2*i;j<=Max/2;j+=i){ 13 f[j]=f[j]*(i-1)/i; 14 prime[j]=1; 15 } 16 } 17 } 18 int main() 19 { 20 freopen("quest.in","r",stdin); 21 freopen("quest.out","w",stdout); 22 shaifa(); 23 int T; scanf("%d",&T); 24 while(T--){ 25 ans=cnt=maxn=0; 26 scanf("%d%d",&n,&m); 27 for(int i=1;i<=n;i++){ 28 int x; scanf("%d",&x); 29 num[x]++; maxn=max(maxn,x); 30 } 31 for(int i=1;i<=maxn;i++){ 32 cnt=0; 33 for(int j=i;j<=maxn;j+=i) cnt+=num[j]; 34 if(cnt>=m) ans=max(ans,f[i]); 35 num[i]=0; 36 } 37 printf("%d\n",ans); 38 } 39 return 0; 40 }
------------------------------------------------------------华丽的分割线----------------------------------------------------------------------
题解(假):首先,在n、m<=100时,我们可以暴力使枚举l到r,然后用矩阵快速幂求出每个点的答案,最后加起来即可,时间复杂度约为为O(NMlogK)。
代码如下:当做矩阵快速幂的练习吧,最蠢的暴力- -
1 #include<cstdio> 2 #include<iostream> 3 #include<cstring> 4 #define Max 100001 5 #define Mod 10007 6 using namespace std; 7 int n,m,l,r; 8 long long K,cnt[Max]; 9 struct mat{long long num[2][2];}a,b,c,ret; 10 mat mult(mat x,mat y){ 11 memset(c.num,0,sizeof(c.num)); 12 for(int i=0;i<2;i++) 13 for(int j=0;j<2;j++) 14 for(int k=0;k<2;k++) 15 c.num[i][j]=(c.num[i][j]+x.num[i][k]*y.num[k][j])%Mod; 16 return c; 17 } 18 mat matmod(mat x,long long k){ 19 memset(ret.num,0,sizeof(ret.num)); 20 for(int i=0;i<2;i++) ret.num[i][i]=1; 21 while(k){ 22 if(k&1) ret=mult(ret,x); 23 k>>=1; x=mult(x,x); 24 } 25 return ret; 26 } 27 long long Get(long long k){ 28 memset(a.num,0,sizeof(a.num)); 29 memset(b.num,0,sizeof(b.num)); 30 a.num[0][1]=1; a.num[1][0]=1; a.num[1][1]=1; 31 b.num[0][0]=1; b.num[1][0]=1; 32 a=matmod(a,k-1); 33 a=mult(a,b); 34 return a.num[0][0]; 35 } 36 int main() 37 { 38 freopen("rabbit.in","r",stdin); 39 freopen("rabbit.out","w",stdout); 40 scanf("%d%d",&n,&m); 41 for(int i=1;i<=n;i++) cnt[i]=1; 42 for(int i=1;i<=m;i++){ 43 long long ans=0; 44 scanf("%d%d%lld",&l,&r,&K); 45 if(K==0){ 46 for(int j=l;j<=r;j++) 47 if(cnt[j]!=0) ans=(ans+Get(cnt[j]))%Mod; 48 printf("%lld\n",ans); 49 } 50 else for(int j=l;j<=r;j++) cnt[j]+=K; 51 } 52 return 0; 53 }
题解:那么,在n、m<=100000时我们该怎么做呢,很明显,我们可以用线段树,但是我还没调出来。。。也许明天会记得填坑吧,大概。
转载于:https://www.cnblogs.com/Beginner-/p/7226326.html
2017省夏令营Day7相关推荐
- 2017福建夏令营Day7(数论)
埃匹希斯水晶 (apexis) ⼤家都知道,卡德加是⼀个神奇的法师. 有⼀天,他发现了⼀种可以作⽤在埃匹希斯⽔晶上的魔法:在左右两个 祭坛上放⼀定量的⽔晶,然后施放⼀个法术,左边⼀堆的⽔晶数量会变成原 ...
- 计算机在生物工程中的应用 华中农业大学,华中农业大学生物信息、生物统计、农业信息工程 2017年夏令营通知...
华中农业大学生物信息.生物统计.农业信息工程 2017年夏令营通知 一.华中农业大学信息学院简介 华中农业大学信息学科发展历史悠久.80年代中后期,学校在基础课部设立了计算机教研室.1996年9月开始 ...
- 人民大学计算机夏令营试题,中国人民大学信息学院2017年夏令营机试题目(考生回忆版)...
中国人民大学信息学院2017年夏令营机试题目(考生回忆版) (6页) 本资源提供全文预览,点击全文预览即可全文预览,如果喜欢文档就下载吧,查找使用更方便哦! 11.90 积分 中国人民大学信息学院20 ...
- 彩色方块 2017信息学夏令营第二场
NKOJ 4223 彩色方块 问题描述 何老板最近在玩一款叫"彩色方块"的小游戏,游戏虽然简单,但何老板仍旧乐此不疲. 游戏中有n个彩色方块成一排,方块的颜色用字母表示,给出目标排 ...
- 理发时间 2017信息学夏令营第三场
NKOJ4231 理发时间 时间限制 : 1000MS 空间限制 : 165536 KB 问题描述 何老板开了一家理发店,店里有m名发型师. 今天,信竞班的n名同学一起同时来到何老板的店,每个同学都要 ...
- 攻防世界 sql注入
攻防世界supersql(兼buuctf随便注) 首先添加一个单引号1' --+发现出错. error 1064 : You have an error in your SQL syntax; che ...
- 浙江大学计算机保研条件_浙江大学保研有哪些步骤?
----------------------- 2020年5月更新 又到了一年一度的保研关键点,欢迎大家关注我的公众号(公众号:查尔de理想生活)留言提问,或者知乎私信咨询. ----------- ...
- 吉林大学计算机科学与技术学院推免,2017吉林大学计算机科学与技术大学生推免夏令营...
2018考研交流群586254585 为了促进国内高校优秀大学生之间的学术交流,激发学生学习热情,提供开展高水平科学研究的机会.特此举办"2017年吉林大学计算机科学与技术学院优秀大学生暑期 ...
- Loj #6077. 「2017 山东一轮集训 Day7」逆序对
Loj #6077. 「2017 山东一轮集训 Day7」逆序对 Solution 令fi,jf_{i,j}fi,j表示前iii个数产生jjj个逆序对的方案数,每次考虑把i+1i+1i+1加入,有i ...
最新文章
- 如何快速构建服务发现的高可用能力
- Ubunt pip升级到指定版本
- boost::hana::trait用法的测试程序
- mybatis默认的数据源连接池(PooledDataSource和UnPooledDataSource)
- Cewolf使用小结
- 腾讯关系型数据库达成“双百”里程碑——6大企业级MySQL特性全面解析
- 使用JBoss EAP 7的HTTP / 2
- 记录关于vs2008 和vs2015 的报错问题
- gitlab bash_如何编写Bash一线式以克隆和管理GitHub和GitLab存储库
- python做var模型_【Python金融量化】VaR系列(五):Copula模型估计组合VaR-阿里云开发者社区...
- python中center()函数的用法
- select在各个浏览器中的兼容性问题
- 计算机图形学-消隐算法
- 如何查看源代码及一些代码的解释
- Word转PDF出现“错误,未定义书签”
- 五、Lua 变量的学习
- 技术需求文档,应当这么写!
- AR增强现实的三大关键技术
- 名词解释第七十讲:基金会
- oracle学习第二天
热门文章
- python设置字符编码_python字符编码问题一则
- 405 not allowed什么意思_clock是什么意思
- 大文件上传 之 改版了的SlickUpload.HttpUploadModule(Krystalware.SlickUpload.dll)
- 必须声明标量变量 @列名
- linux安装chrome_chrome 的安装及出现密钥问题(ubuntu16.04)
- 决策树(八)--随机森林及OpenCV源码分析
- Java项目迁移到uap上_Tomcat启动,不能加载项目问题。
- python面对对象的编程语言_python面对对象编程
- C语言——指针篇(二)指针和数组之内存分配和初始化
- linux查看php执行用户,在浏览器中打开php文件时,是Linux中的哪个用户执行的?...