noip模拟12[简单的区间·简单的玄学·简单的填数]
noip模拟12 solutions
这次考试靠的还是比较好的，但是还是有不好的地方，
为啥嘞？？因为我觉得我排列组合好像白学了诶，文化课都忘记了
正难则反！！！！！！！！
害没关系啦，一共拿到了\(120pts\)，其实距离我的理想分数还差那么\(100pts\)
具体是这样的，第一题AC，第二题10，第三题10
下次要把知识都回忆一下，比如这次用到的欧拉定理，差一点就忘记了
noip模拟13！！！200分！！
·
·
·
T1 简单的区间
哈哈哈这个题是我这几次考试中最成功的一道了，所以我一定要好好讲一讲我的STL算法
（因为正解好像是主席树 逃）（我现在对这树那树的非常恶心）
首先我第一眼看到这个题面的时候，我直接就想到了昨天\(noip11\)那场的\(T3\)
同样是求有关于最大值的结果，还是在区间上，那必须是单调栈走起啊，好像有人还不会诶
priority_stack
这里是利用数组来实现的单调栈，相比普通的STL栈来说要好理解
sta为栈数组，tot为栈顶
每一次向栈中加入值的时候，都把栈中小于等于它的数弹出，并将这些弹出的元素右边界赋值为i-1
注意判断栈中是否还有元素！！1
弹完之后将现在插入的值的左边界赋值为sta[tot]+1;就是栈顶元素的右边一位
最后，将栈中元素全部弹出，右边界复制为n
for(re i=1;i<=n;i++){
while(tot&&a[i]>=a[sta[tot]])r[sta[tot]]=i-1,tot--;
//弹出+右边界
l[i]=sta[tot]+1;
//左边界
sta[++tot]=i;
//加入栈中
}
while(tot)r[sta[tot]]=n,tot--;
//给右边界赋值
然后我们找到左右边界之后，再回来看这个题，
求的是区间和减去最大值\(\mod k==0\)所以，既然涉及到区间求和，那前缀和是必不可少的
设前缀和数组为\(fro[]\),左边界为\(i\)，右边界为\(j\),此时区间最大值的坐标为\(x\)
那么我们有一个结论\(fro[j]-fro[i-1]-a[x]=0\)（此处省略一堆+k、%k）
既然这样我们为了好处理这个式子，我们把所有的\(fro[i]\)、\(a[i]\)对\(k\)取模
接下来我们就可以利用一个\(vector\)数组来存储每一个状态
vec[tmp]，是一个\(vector\)类型表示余数为tmp的坐标都有哪些，
在插入值的时候，我们只需要从1～n遍历一遍，然后vec[tmp].push_back(i)就好了
（因为是顺序插入，所以每一个\(vector\)中都是排好序的坐标）
然后我们在左右区间内选取较小的那一段来进行枚举（假设选取左边较小）
那么右边界的前缀和就可以由\(fro[i-1]+a[x]\)来求出，那么我们就可以直接在\(vector\)中查找在x～r这个区间内的数有多少
这个操作可以通过STL内置函数lower_bound和upper_bound来实现
如果是右区间较短，也同理，有些细节在代码里展示
AC_code
#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int N=3e5+5;
const int M=1e6+5;
int n,k;
int a[N],fro[N];
int ys[M],ji[M],cnt;
int l[N],r[N],sta[N],tot;
vector<int> vec[M];
ll ans;
void sol(int x,int l,int r){
if(x-l<r-x){
for(re i=l;i<=x;i++){
int tmp=(fro[i-1]+a[x])%k;//注意此处fro[i-1];
int lh=lower_bound(vec[tmp].begin(),vec[tmp].end(),x)-vec[tmp].begin();
int rh=upper_bound(vec[tmp].begin(),vec[tmp].end(),r)-vec[tmp].begin();
ans+=rh-lh;
}
}
else{
for(re i=x;i<=r;i++){
int tmp=(fro[i]+k-a[x]%k)%k;//a[x]%k别忘记
int lh=lower_bound(vec[tmp].begin(),vec[tmp].end(),l-1)-vec[tmp].begin();//此处l-1
int rh=upper_bound(vec[tmp].begin(),vec[tmp].end(),x-1)-vec[tmp].begin();//x-1
ans+=rh-lh;
}
}
}
signed main(){
scanf("%d%d",&n,&k);
for(re i=1;i<=n;i++){
scanf("%d",&a[i]);
fro[i]=(fro[i-1]+a[i])%k;
ys[fro[i]]++;
if(ys[fro[i]]==1)ji[++cnt]=fro[i];//记录出现过的余数
}
for(re i=1;i<=cnt;i++){
vec[ji[i]].reserve(ys[ji[i]]+5);//限制空间
}
vec[0].push_back(0);
for(re i=1;i<=n;i++)
vec[fro[i]].push_back(i);
for(re i=1;i<=n;i++){
while(tot&&a[i]>=a[sta[tot]])r[sta[tot]]=i-1,tot--;
l[i]=sta[tot]+1;
sta[++tot]=i;
}
while(tot)r[sta[tot]]=n,tot--;
for(re i=1;i<=n;i++)sol(i,l[i],r[i]);
printf("%lld",ans-n);//最后方案数减去n，因为左边界不能等于右边界
}
你是不是看到我的代码里好像有一些并没有什么用的vector操作
那些操作可能在这个题中没啥用，但是这个很重要。看看这篇博客
T2 简单的玄学
这个题就是排列组合嘛，数学课上我就一直瞎嚷嚷：正难则反！！正难则反！！！
可是到了真正用到这个做法的时候，我倒是看不出来了，真是的
就这个题统计至少两个数相同的概率，我在考场上苦思冥想，想出来了正着做的办法，然后10pts
考完一看题解，正难则反，直接1减去两两都不相同的概率不就好了？？？
看的我人都傻了，咋考场上就想不出来啊，要是想出来了，起码也得有50pts啊
两两不相同的概率是
\(\dfrac{A_{2^n}^{m} }{2^{nm}}=\dfrac{\frac{2^n!}{(2^n-m)!}}{2^{nm}}=\dfrac{\prod_{i=2^n-m+1}^{2^n}}{2^{nm}}=\dfrac{\prod_{i=2^n-m+1}^{2^n-1}}{2^{n(m-1)}}\)
这样，目前在这个式子上来看，是没有啥可以约掉的啦
题目要求我们对这分子分母进行约分，仔细观察发现，这个分母的因子只有2
所以我们只需要找到分母中有多少个\(2\)就\(OK\)啦
现在有一个结论：\(2^n-a\)和\(a\)中，因子2的个数是一致的，
因为它们两个加起来是2的幂，所以我们可以知道，这两个数同时除以相同个数的2后，
一定相加是偶数，所以只能是都是奇数或者都是偶数，都是奇数那就符合这个结论
都是偶数，那就继续除，除到是奇数为止
所以我们要寻找\(2^n-m+1\)~~~~~\(2^n-1\)这些数相乘的结果中的2的个数，
等价于找\(1～m-1\)中的2的个数，这里确实有一个\(O(logm)\)的算法，并且极其好想
O(logm) get sum of 2
这里主要是利用2的性质
m个数中，有m/2个数可以整除2，那么就有这么多个2
同理，有m/4个数可以整除4,那么就再有这么多2
........
然后就这样一直除下去，得到fz，
for(ll i=2ll;i<m;i*=2ll)
fz+=(m-1)/i,fz%=(mod-1);
那么此时我们发现如果\(m>mod\)的话，那就一定有一个数是mod的倍数，那分子就变成0咯
但是此时的答案并不是1,因为%mod之后为0,并不代表原来为0
所以这样的话我们就计算分母就好啦，别忘记除掉约去的2
那要是m<mod的话，暴力计算分子，输出答案
答案是1减去你算出来的值啊啊啊啊啊啊啊啊啊啊
AC_code
#include<bits/stdc++.h>
using namespace std;
#define re register int
#define ll long long
const int mod=1e6+3;
ll n,m;
ll fz;
ll a,b;
ll ksm(ll x,ll y){
ll ret=1;
while(y){
if(y&1)ret=ret*x%mod;
x=x*x%mod;
y>>=1;
}
return ret;
}
signed main(){
scanf("%lld%lld",&n,&m);
for(ll i=2ll;i<m;i*=2ll)
fz+=(m-1)/i,fz%=(mod-1);
if(m>mod){
ll tmp=(n%(mod-1))*((m-1)%(mod-1))%(mod-1);
tmp=(tmp-fz+mod-1)%(mod-1);
ll bas=ksm(2ll,tmp);
printf("%lld %lld",bas,bas);
return 0;
}
a=1;
for(ll i=1;i<m;i++)a=a*(ksm(2ll,n%(mod-1))+mod-i)%mod;
a=a*ksm(ksm(2ll,fz),mod-2)%mod;
b=ksm(2ll,((n%(mod-1)*(m-1)%(mod-1))%(mod-1)+mod-1-fz)%(mod-1));
printf("%lld %lld",(b+mod-a)%mod,b);
}
T3 简单的填数
这个题考场上是真没有思路好吧，真的惨暴了
但是我在考场上确实想到了要正着走一遍，倒着走一遍，只是这个过程有点小小的复杂
然后我在考场上直接放弃了，对着我那第一题的对拍程序高兴去了
这个题主要是维护一个up和dw二元组(x,len)，存储这个位置目前最大的值及其个数，最小值及其个数
所以我们就有了一个转移，
当len==2时，up更新
当len==5时，dw更新
但是注意，在原序列上如果a[i]上有值的话：
如果当前值和a[i]相同的话，那就不需要做任何更改
如果不相等的话，那就直接赋值成a[i],
up.len=2(因为我要让up越大越好，所以把他赋值成要更改的边界)
dw.len=1(要最小嘛)
然后就剩下倒序走一遍然后得到答案序列啦，这是最难的啦
我是记录了一个now表示从上一个标记的a[i]开始5个一组的数，目前是多少，num表示现在的now有几个
然后我们的\(ans[i]=\min({now},{up[i].x})\)
如果遇到下一个标记好的a[i],就\(now=a[i]，num=0\)就行了
最后倒序输出一下这个序列就好啦
记得判断的时候有两种情况，一个是up<a[i],dw>a[i];
还有一个就是，最后一个长度不够2，输出-1;
AC_code
#include<bits/stdc++.h>
using namespace std;
#define re register int
const int N=2e5+5;
int n,a[N];
struct node{
int x,len;
}up[N],dw[N];
int ans[N];
signed main(){
scanf("%d",&n);
for(re i=1;i<=n;i++)scanf("%d",&a[i]);
up[1].x=1;up[1].len=1;
dw[1].x=1;up[1].len=1;
for(re i=2;i<=n;i++){
up[i].x=up[i-1].x;up[i].len=up[i-1].len+1;
if(up[i].len>2)up[i].x++,up[i].len=1;
dw[i].x=dw[i-1].x;dw[i].len=dw[i-1].len+1;
if(dw[i].len>5)dw[i].x++,dw[i].len=1;
if(a[i]){
if(up[i].x<a[i]||dw[i].x>a[i]){
printf("-1");return 0;
}
if(up[i].x>a[i])up[i].x=a[i],up[i].len=2;
if(dw[i].x<a[i])dw[i].x=a[i],dw[i].len=1;
}
}
if(up[n].len==1)up[n].x--,up[n].len=up[n-1].len+1;
if(up[n].x<dw[n].x){
printf("-1");return 0;
}
printf("%d\n",ans[n]=up[n].x);
int now=ans[n],num=1;
for(re i=n-1;i>=1;i--){
if(a[i]){
if(now>a[i])now=a[i],num=0;
}
if(num==5)now--,num=0;
num++;
ans[i]=min(now,up[i].x);
}
for(re i=1;i<=n;i++)printf("%d ",ans[i]);
}
完事啦
