bzoj4447[Scoi2015]小凸解密码

4447: [Scoi2015]小凸解密码

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 150  Solved: 58
[Submit][Status][Discuss]

Description

小凸得到了一个密码盘，密码盘被等分成N个扇形，每个扇形上有一个数字(0～9)，和一个符号（“+”或"*")

密码盘解密的方法如下：

首先，选择一个位置开始，顺时针地将数字和符号分别记在数组A和数组C巾

解密的方法如下

B0=A0

当x>0时：

若Cx为“+”，Bx=(Ax+Ax-1)%10，注意：x-1是下标值

若Cx为“*”，Bx= (Ax×Ax-1)%10,注意：x-1是下标值

操作完成后，可以得到一个长度为n的数组B，然后以B0为起点将B数组顺时针写成一个环，解密就完成

了，称得到的环为答案环。

现在小凸得到了一份指令表，指令表上有2种操作。

一种指令是修改操作，即改变原来密码盘上一个位置的数字和符号。

另一种指令是询问操作，具体如下：

首先从指令给出的位置开始完成解密，得到答案环。

答案环上会有一些0连在一起，将这些连在一起的0称为零区间，找出其中距离B0最远的那个零区间，输

出这个距离。

零区问和B0的距离定义为：零区问内所有0到B0距离中的最小值。

Input

第1行包含2个整数n,m，代表密码盘大小和指令个数

接下来n行，每行包含1个整数和1个字符，按顺时针顺序给出了密码盘上的数组和符号

接下来m行，依次给出指令

每行第1个整数代表指令类型

若第1个墼数为1，代表本行对应指令为修改操作，之后依次有2个整数pos，num和1个字符opt，分别

代表修改的位置，以及修改后该位置的数字和字符

若第1个整数为2，代表本行对应指令位询问操作，之后有1个整数pos，代表本次操作中解密的开始位置

密码盘上的位置标号为0到n-l

数据保证合法，即数据中0≤pos<N，0≤num≤9，opt为“+”或“*”

Output

对于每个询问操作1行，输出答案，若答案环上没有0，输出-1

Sample Input

5 8
0 *
0 *
0 *
0 *
0 *
2 0
1 0 1 +
1 2 1 +
2 3
1 1 1 +
1 3 1 +
1 4 1 +
2 4

Sample Output

0
2
-1

HINT

第1个询问，答案环为[0,0,0,0,0]，仅有1个零区间，且B0在其中，所以距离是0

对于第2个询问，答案环为[0,0,1,0,l]，有2个零区间，(0，1)和B0距离是o，(3，3)和B0距离是2，故答案为2

对于第3个询问，答案环为[1，2，2，2，2]，没有零区间，答案是-1

对于100%数据，5 <=n，m≤10^5

 

 

 

UPDATE：

尴尬地发现其实方法有点小BUG，WA不是因为写错了。

 

错误原因是这样的：

如果区间端点为0，左右区间拼成环后构成新的跨越左右区间的0区间，这时候直接取$min(dis_l,dis_r)$ 不一定是答案，因为有可能左右区间的第二远的0区间比上面取min得到的答案更优，举一个例子：

权值 0 0 0 0 0 1 1 1  0   0   0   1   0   0

位置 1 2 3 4 5 6 8 9 10 11 12 13 14 15

假设查询的是5位置，那么取min得到的答案就是0，而右边区间第二远的0区间的距离是4

 

所以如果出现上述的区间端点为0的情况，就需要查询新的左右区间的最远0区间的距离，即删去跨越左右区间的0区间，再做一次查询最远0区间的操作。

这样做以后虽然复杂度没变，但是常数扩大了很多，不知道能否过掉，有极大可能TLE 几组数据，所以代码就懒得改了

 

 

 

发现了更好的办法

做线段树的时候的时候可以直接维护维护离左右端点最远的和次远0区间的距离

 

下面的blog就是这种做法，可以参照

https://www.cnblogs.com/f321dd/p/6403097.html

 

 

 

以下是原方法

 

 

脑残的 二分+线段树

化环成链，倍增区间

查询的时候分成左右两段区间查询

二分查询每个点最左边的0在哪个位置，最右边的0在哪个位置，用线段树check

找到最左和最右的0之后，线段树查询当前0所在的区间的另一端点，距离为dis

如果两区间不为查询点的端点都为0，那么查询到的最远的0都是端点了。

说明左右区间拼成环之后，可以构成一个新的跨越了左右区间的0区间，答案取$min{(dis_l,dis_r)}$  否则取$max{(dis_l,dis_r)}$

WA + TLE
WA估计是写的时候出了点小错误，不想调试啦。
TLE这东西加点常数优化应该是可以卡过的
好像还有set作法貌似挺简单，直接存全为0的区间即可

 #include<bits/stdc++.h>
 #define ls u<<1
 #define rs ls|1
 #define N 200010
 using namespace std;
 int n,m,a[N],b[N],pd[N<<],lz[N<<],lx[N<<],rx[N<<];char s[N];
 void pushup(int u,int l,int r){
     int mid=(l+r)>>;
     pd[u]=pd[ls]|pd[rs];
     lx[u]=lx[ls];rx[u]=rx[rs];
     if(lx[u]==mid-l+)lx[u]+=lx[rs];
     if(rx[u]==r-mid)rx[u]+=rx[ls];
 }
 void build(int u,int l,int r){
     if(l==r){
         if(!b[l])pd[u]=lx[u]=rx[u]=;
         else pd[u]=lx[u]=rx[u]=;
         return;
     }
     int mid=(l+r)>>;
     build(ls,l,mid);
     build(rs,mid+,r);
     pushup(u,l,r);
 }
 void update(int u,int l,int r,int p){
     if(l==r){
         if(!b[l])pd[u]=lx[u]=rx[u]=;
         else pd[u]=lx[u]=rx[u]=;
         return;
     }
     int mid=(l+r)>>;
     if(p<=mid)update(ls,l,mid,p);
     else update(rs,mid+,r,p);
     pushup(u,l,r);
 }

 int query(int u,int L,int R,int l,int r){
     if(l<=L&&R<=r)return pd[u];
     int mid=(L+R)>>,ret=;
     if(l<=mid)ret|=query(ls,L,mid,l,r);
     if(r>mid)ret|=query(rs,mid+,R,l,r);
     return ret;
 }
 int asklx(int u,int L,int R,int l,int r){
     if(l==L&&R==r)return lx[u];
     int mid=(L+R)>>,ret=;
     if(r<=mid)return asklx(ls,L,mid,l,r);
     if(l>mid)return asklx(rs,mid+,R,l,r);
     ret+=asklx(ls,L,mid,l,mid);
     if(!ret)return ;
     if(ret==mid-l+)ret+=asklx(rs,mid+,R,mid+,r);
     return ret;
 }
 int askrx(int u,int L,int R,int l,int r){
     if(l==L&&R==r)return rx[u];
     int mid=(L+R)>>,ret=;
     if(r<=mid)return askrx(ls,L,mid,l,r);
     if(l>mid)return askrx(rs,mid+,R,l,r);
     ret+=askrx(rs,mid+,R,mid+,r);
     if(!ret)return ;
     if(ret==r-mid)ret+=askrx(ls,L,mid,l,mid);
     return ret;
 }
 inline int solve(int x){
     static int t1,t2,x1,x2,L,R,l,r;
     int mid=(+n)>>;
     if(x<mid)x+=n;
     L=x-mid+;R=x+mid;
     if(!((+n)&))R--;
     t1=t2=;
     l=x,r=R;
     while(l<=r){
         mid=(l+r)>>;
         if(query(,,n<<,mid,R))t2=mid,l=mid+;
         else r=mid-;
     }
     l=L,r=x;
     while(l<=r){
         mid=(l+r)>>;
         if(query(,,n<<,L,mid))t1=mid,r=mid-;
         else l=mid+;
     }
     if(!t1&&!t2)return -;
     if(t1==x&&t2==x)return ;
     x1=t1+asklx(,,n<<,t1,x)-;
     x2=t2-askrx(,,n<<,x,t2)+;
     if(x1==x&&x2==x)return ;
     if(t1==L&&t2==R)return min(x2-x,x-x1);
     return max(x2-x,x-x1);
 }
 inline void change(int p,int op){
     static int k;k=p-;if(!k)k=n;
     if(op)b[p]=(a[k]+a[p])%;
     else b[p]=1ll*a[k]*a[p]%;
     if(p>n)b[p-n]=b[p];
     else b[p+n]=b[p];
 }
 int main(){
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;i++)
     scanf("%d %c",&a[i],&s[i]);
     for(int i=;i<=n;i++){
         int j=i-;if(!j)j=n;
         if(s[i]=='+')b[i]=(a[i]+a[j])%;
         else b[i]=1ll*a[i]*a[j]%;
         b[i+n]=b[i];s[i+n]=s[i];a[i+n]=a[i];
     }
     build(,,n<<);
     int x,y,z,p;char op;
     while(m--){
         scanf("%d%d",&x,&y);++y;
         if(x==){
             scanf("%d %c",&z,&op);
             a[y]=a[y+n]=z;s[y]=s[y+n]=op;
             change(y,s[y]=='+'?:);
             update(,,n<<,y);
             update(,,n<<,y+n);
             change(y+,s[y+]=='+'?:);
             update(,,n<<,y+);
             p=y+>n?y+-n:y++n;
             update(,,n<<,p);
         }
         else printf("%d\n",solve(y));
     }
     return ;
 }