博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
Jzoj5449【NOIP2017提高A组冲刺11.4】Pacifist
阅读量:5094 次
发布时间:2019-06-13

本文共 1219 字,大约阅读时间需要 4 分钟。

papyrus 喜欢谜题... 来解一道如何?
在你面前有一个被加密了的数组,其原数组是一个等差序列,你面前的则是将原数组中的所有数字都对m 取模再打乱后而得到的新数组
papyrus 给你出的谜题就是还原出原等差序列a

保证数据有解,而且因为papyrus 喜欢质数,所以他给你出的谜题中的m 一定是质数

今天题目还挺好玩的,虽然第二题因为数据太多OJ最后挂了全场爆0

但是还是把第一题切到了80分,后面卡卡常数就过去了

题意:给出一个模意义下的等差序列叫你求原序列

方法非常之多:有对称法,平方法,暴力枚举(雾)

我是用的是枚举+数论法

我们发现首项f一定在序列中,假设我们知道了首项那么就会有

nf+dn(n-1)/2=S (Mod m){S=Σai}

移项得到 n(n-1)d=2(S-nf) (Mod m)

这样就可以用一个扩展gcd求一下线性模方程就好了

让后考虑求出来的d是否合法

我当时考虑是使用随机选一些i并判断(id+f)%m是否在a中,如果是就接受这个答案

但是这样是有问题的,我被一个很强的数据卡了

没办法只好改为暴力枚举每一个i

但是这样似乎是可行的,因为很多情况下,可以直接break,再加上把set换成lowerbound就可以卡过

#include
#include
#include
#include
#include
#include
#define LL long longusing namespace std;LL exgcd(LL a,LL b,LL& x,LL& y){ if(b){ LL r=exgcd(b,a%b,y,x); y-=x*(a/b); return r; } else { x=1; y=0; return a; }}LL mod(LL a,LL b,LL p){ LL x,y,r=exgcd(a,p,x,y); x=(x+p)%p; return x*(b/r);}LL inv(LL a,LL p){ LL x,y,r=exgcd(a,p,x,y); return (x+p)%p;}LL in;int n,s[100010],p,b=0;int cal(LL f){ LL x,y,r=exgcd(n*(n-1ll)%p*in%p,p,x,y); return (x*((b-f*n%p+p)%p)%p+p)%p;}bool ok(LL d,int f){ for(int i=1;i
>1))&&mxd>d&&ok(d,s[i])) if(d

转载于:https://www.cnblogs.com/Extended-Ash/p/7800613.html

你可能感兴趣的文章
Python基础-数据类型
查看>>
Gym - 100989G (二分法)
查看>>
第三周学习进度条
查看>>
spark 参数调优
查看>>
精彩博文收集目录索引(程序猿就是我)
查看>>
搭建日志环境并配置显示DDL语句
查看>>
SelectorChapek(选择器的快速创建)
查看>>
使用Java合并图片、修改DPI
查看>>
JPush极光推送 Java调用服务器端API开发
查看>>
Virtualbox 无法加载usb解决办法
查看>>
Linux常用命令
查看>>
ubuntu 下LAMP服务器环境搭建
查看>>
centos7 安装SVN
查看>>
【ZH奶酪】如何用sklearn计算中文文本TF-IDF?
查看>>
USACO Shaping Regions,难题,离散化,矩形切割,逆序染色
查看>>
iis 还原配置
查看>>
设计模式——门面模式
查看>>
自己动手打造工具系列之自动刷新简历
查看>>
Sqlserver2005附加数据库为只读的解决方法
查看>>
[BZOJ 1296] 粉刷匠
查看>>