Jzoj5449【NOIP2017提高A组冲刺11.4】Pacifist

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<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<time.h>
#include<set>
#define LL long long
using 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<n;++i)if(*lower_bound(s+1,s+1+n,(d*i+f)%p)!=(d*i+f)%p) return 0;return 1;
}
int main(){freopen("pacifist.in","r",stdin);freopen("pacifist.out","w",stdout);scanf("%d%d",&p,&n); const int M=p;for(int i=1;i<=n;++i){ scanf("%d",s+i);b=(b+s[i])%M; }sort(s+1,s+1+n);if(s[1]==s[2]) return 0&puts("-1");if(n==p){printf("0 1");return 0;} in=inv(2,p);int mxd=1<<30,mf;for(int i=1;i<=n;++i){int d=cal(s[i]);if(d&&(d<=(p>>1))&&mxd>d&&ok(d,s[i]))if(d<mxd){ mxd=d; mf=s[i]; }}printf("%d %d\n",mf,mxd);
}

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

Jzoj5449【NOIP2017提高A组冲刺11.4】Pacifist相关推荐

JZOJ 5441. 【NOIP2017提高A组冲刺11.1】序列
Description 给定一个1~n的排列x,每次你可以将x1~xi翻转.你需要求出将序列变为升序的最小操作次数.有多组数据. Input 第一行一个整数t表示数据组数. 每组数据第一行一个整数n, ...
jzoj_5455. 【NOIP2017提高A组冲刺11.6】拆网线
Description 企鹅国的网吧们之间由网线互相连接,形成一棵树的结构.现在由于冬天到了,供暖部门缺少燃料,于是他们决定去拆一些网线来做燃料.但是现在有K只企鹅要上网和别人联机游戏,所以他们需要把 ...
[JZOJ5459]【NOIP2017提高A组冲刺11.7】密室
Description 小X 正困在一个密室里,他希望尽快逃出密室. 密室中有N 个房间,初始时,小X 在1 号房间,而出口在N 号房间. 密室的每一个房间中可能有着一些钥匙和一些传送门,一个传送门会 ...
JZOJ 5458. 【NOIP2017提高A组冲刺11.7】质数
Description 小X 是一位热爱数学的男孩子,在茫茫的数字中,他对质数更有一种独特的情感.小X 认为,质数是一切自然数起源的地方. 在小X 的认知里,质数是除了本身和1 以外,没有其他因数的数 ...
JZOJ 5455. 【NOIP2017提高A组冲刺11.6】拆网线
Description 企鹅国的网吧们之间由网线互相连接,形成一棵树的结构.现在由于冬天到了,供暖部门缺少燃料,于是他们决定去拆一些网线来做燃料.但是现在有K只企鹅要上网和别人联机游戏,所以他们需要把 ...
JZOJ 5445. 【NOIP2017提高A组冲刺11.2】失格
Description 胆小鬼连幸福都会害怕,碰到棉花都会受伤,有时还被幸福所伤. --太宰治<人间失格> 回顾我的一生,一共有n个事件,每一个事件有一个幸福值p_i. 我想用n-1条线把 ...
JZOJ5445. 【NOIP2017提高A组冲刺11.2】失格
Description 胆小鬼连幸福都会害怕,碰到棉花都会受伤,有时还被幸福所伤. --太宰治<人间失格> 回顾我的一生,一共有n个事件,每一个事件有一个幸福值p_i. 我想用n-1条线把 ...
JZOJ 5466. 【NOIP2017提高A组冲刺11.9】玩游戏
Description 小A得了忧郁综合症,小B正在想办法开导她. 机智的小B决定陪着小A玩游戏,他从魔法的世界里变出一张无向联通图,每条边上都有边权.小B定义一条路径的权值为所有经过边中的最大权值, ...
JZOJ 5462. 【NOIP2017提高A组冲刺11.8】好文章
Description nodgd写了一篇文章,自认为这是一篇好文章.nodgd的文章由n个小写英文字母组成.文章的一个子串指的是文章中的一段连续的字母,子串的长度就是这一段的字母个数.nodgd在文 ...

Jzoj5449【NOIP2017提高A组冲刺11.4】Pacifist

Jzoj5449【NOIP2017提高A组冲刺11.4】Pacifist相关推荐

最新文章

热门文章