[cf1428H]Rotary Laser Lock

为了方便,定义0区间初始左端点所在的编号为0,其余位置顺时针依次为$[1,nm)$

考虑对0区间顺时针旋转,记$s_{i}$为0区间左端点旋转到$i$时的答案(约定$s_{i+n}=s_{i}$)

性质1:若$s_{i}-s_{i-1}=1$且$s_{i+1}-s_{i}\ne 1$(其中$0\le i<nm$),则存在一段区间以$i$为左端点

记$a_{i}$为除0以外是否存在区间包含位置$i$(存在为1),不难得到$s_{i}-s_{i-1}=a_{i+m-1}-a_{i-1}$

代入条件,由$s_{i}-s_{i-1}=1$得$a_{i+m-1}=1$且$a_{i-1}=0$,同时$a_{i+m}\ne 1$或$a_{i}\ne 0$

对后者分类讨论,以第一种情况为例,即$a_{i+m-1}=1$且$a_{i+m}\ne 1(=0)$,显然得证

同时,显然$s_{i}$中必然存在这样的$i$,记为$I$并对其余区间分别逆时针旋转一次

性质2:记$s'_{i}$为区间$j$逆时针旋转后的$s_{i}$,则$j$以$I$为左端点当且仅当$s'_{I}-s'_{I-1}\ne 1$

若$j$以$I$为左端点,注意到这会使得$a_{I-1}$变为1,进而即$s'_{I}-s'_{I-1}\ne 1$

另一方面,其需要影响$a_{I-1}$或$a_{I+m-1}$,且前者初始为0、后者初始为1,因此仅能以$I$为左端点

 

综上,考虑重复以下过程直至得到答案——

1.暴力找到$I$,同时将第0个区间移动到以$I$为左端点的位置

2.枚举所有未确定位置的区间,通过4次移动可以判定其是否以$I$为左端点并还原

另外,可以通过合理的操作顺序,使得其同时判定$s_{i}-s_{i-1}$是否为1(不为1即已经不存在以$I$为左端点的区间)

3.对于左端点在$I$上的区间,再逆时针移动$m+1$次(要求不覆盖到$I-1$)

(实际上,这个过程即不断找到连续段的开头)

进一步的,考虑上述过程所需操作次数,对每一步分别考虑:

1.对于第一步,不难发现第一次找到后至多再转一圈,共计$2nm+2(n-1)$次

2.对于第二步,一共有$n-1$轮,第$i$轮至多有$n-i$个区间,共计$4\sum_{i=1}^{n-1}(n-i)$

3.对于第三步,若第二步中最后一次不还原,每一次实际仅额外产生$m-1$次,共计$(n-1)(m-1)$次

求和后,共计需要$2n^{2}+3nm-n-m-1=25879$次,无法通过

 

对性质2拓展,每次旋转多个区间,设其构成集合$S$,则$S$中存在以$I$为左端点的区间当且仅当$s'_{I}-s'_{I-1}=1$

以此法不断进行二分,即0区间的移动次数变为$2\log n$,次数降为$2\sum_{i=1}^{n-1}(n-i)+2(n-1)\log n$

同时,由于是不断二分,确定不以$I$为左端点的区间并不需要还原(原来可能会从$I+1$移动到$I$后跳过)

另外,此时第二步还会多出$4(n-1)$次检验,第三步每一次实际额外产生$m+1$次

综上,后两步的次数分别变为$\sum_{i=1}^{n-1}(n-i)+2(n-1)(\log n+2)$和$(n-1)(m+1)$

求和后,共计需要$\frac{n^{2}}{2}+3nm+2n\log n+\frac{11n}{2}-2\log n-m-7=13107$次,可以通过

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 #define N 105
 4 vector<int>v;
 5 int n,m,L,now,ans[N];
 6 void move(int k,int p){
 7     printf("? %d %d\n",k,p);
 8     fflush(stdout);
 9     scanf("%d",&now);
10 }
11 int main(){
12     scanf("%d%d",&n,&m),L=n*m;
13     for(int i=1;i<n;i++)ans[i]=-1;
14     move(0,1);
15     int lst=now,flag=0;
16     for(int i=0;i<(L<<1);i++){
17         move(0,1);
18         if (ans[0]==n-1)continue;
19         if (lst+1==now)flag=1;
20         else{
21             if (flag){
22                 move(0,-1);
23                 while (1){
24                     if (ans[0]==n-1)break;
25                     move(0,-1),lst=now,move(0,1);
26                     if (lst+1!=now)break;
27                     lst=now,move(0,1);
28                     if (lst+1==now){
29                         move(0,-1);
30                         break;
31                     }
32                     move(0,-1),v.clear();
33                     for(int j=1;j<n;j++)
34                         if (ans[j]<0)v.push_back(j);
35                     int l=0,r=(int)v.size()-1,flag=0;
36                     while (l<r){
37                         int mid=(l+r>>1);
38                         if (!flag){
39                             for(int j=l;j<=mid;j++)move(v[j],-1);
40                         } 
41                         else{
42                             for(int j=mid+1;j<=r;j++)move(v[j],1);
43                         }
44                         move(0,-1),lst=now,move(0,1);
45                         if (lst+1==now)l=mid+1,flag=0;
46                         else r=mid,flag=1;
47                     }
48                     for(int j=0;j<=m-flag;j++)move(v[l],-1);
49                     ans[0]++,ans[v[l]]=(i+L-m-1)%L;
50                 }
51                 lst=now,move(0,1);
52             }
53             flag=(lst+1==now);
54         }
55         lst=now;
56     }
57     printf("! ");
58     for(int i=1;i<n;i++)printf("%d ",ans[i]);
59     fflush(stdout);
60     return 0;
61 }
View Code

 

posted @ 2022-01-15 11:52  PYWBKTDA  阅读(113)  评论(0编辑  收藏  举报