[考试反思]0921csp-s模拟测试49：困顿

太弱。还是太弱。

拉不开分差，离第一机房分数线估计还是300多分。

但是，还是要骂：XX出题人。

部分分非常少且没有意义，T1基本只有0/纯暴力20/100三个档，

T2正解是n²但是n³一分不给，还要卡n²的空间，T3的n⁵有51分但没有任何提示，除了送了与正解根本就无关的7分。

而且造的数据还出了锅，这个出题人啊。。。

出题人是有多痛恨人类。

开考看题。

T1上来先看错了题以为是xor，然后就是trie裸题了，然后直接开始打。

很快地打出来发现过不了样例，然后发现是mod。

重新看题，全不会，傻眼了。

想再找一道最简单的题先做稳住，但是无从下手。

刚开始T2有一个n³的思路，发现没有设分，就绝望了。

想了一个半小时多吧，想出了n²，然后发现卡空间。

太弱了没有直接想到正解，zkt就知道直接按x排序，我和脸和toot都是按y排的。

然后大脸秒掉了空间优化，toot放弃了空间优化拿80分走了。

我以为只有50分于是不想放弃，但是剩下两题还没有打，于是走了。

回去弄了一个T1的随机化+测试点分治暴力，然后忘记取模，本来40分的代码炸成0了。

因为没有打正解所以就没有打对拍。。。

然后T3顺手把特殊性质7分拿下，直接放弃。

回T2，对拍，然后空间优化。

思路还可以，可以节约25%的空间，但是忘了vector会翻倍导致空间浪费，于是结果就是和toot没优化一样，80分。

暴力也要打对拍。

对于卡空间的题要精打细算，注意vector的2倍空间。

get黑科技：resize。

还是没有能把所有能拿下的分拿下，还是弱啊。。。

T1：养花

分块。预处理块中对于100001的所有询问的最大答案。

预处理的方法类似于lower_bound但是可以记在数组上就去掉了不必要的log。

开数组lst[i]记录下这个块里小于等于i的值里的最大值。

那么对于每一个询问就可以搞了：mod k的最大值就是max(lst[k-1],lst[2k-1]-k,lst[3k-1]-2k...)

而计算这个式子的复杂度是调和级数$\sum\limits_{i=1}^{n} \frac{n}{i} = ln n$

所以总的复杂度就是$O(\sqrt{n} \times 100001 \times ln 100001 + m \sqrt{n})$

然而其实不是$\sqrt{n}$，会T的，块要开的稍大一点（1300左右）跑的能快一些

 #include<cstdio>
 #include<iostream>
 using namespace std;
 #define S 1234
 #define N 100001
 int lst[],n,m,x[],Ans[][],bl[];
 int main(){
     scanf("%d%d",&n,&m);
     const int bs=(n-)/S+;
     for(int i=;i<=n;++i)scanf("%d",&x[i]);
     for(int i=;i<=bs;++i){
         for(int j=;j<=N;++j)lst[j]=;
         for(int j=(i-)*S+;j<=i*S&&j<=n;++j)lst[x[j]]=x[j],bl[j]=i;
         for(int j=;j<=N;++j)lst[j]=max(lst[j-],lst[j]);
         for(int j=;j<=N;++j)for(int k=j;;k+=j){
             if(k>N){Ans[i][j]=max(Ans[i][j],lst[N]%j);break;}
             Ans[i][j]=max(Ans[i][j],lst[k-]%j);
         }
     }
     for(int t=,l,r,q;t<=m;++t){
         scanf("%d%d%d",&l,&r,&q);int ans=;
         if(bl[l]==bl[r])for(int i=l;i<=r;++i)ans=max(ans,x[i]%q);
         else{
             for(int i=l;i<=bl[l]*S;++i)ans=max(ans,x[i]%q);
             for(int i=(bl[r]-)*S+;i<=r;++i)ans=max(ans,x[i]%q);
             for(int i=bl[l]+;i<bl[r];++i)ans=max(ans,Ans[i][q]);
         }printf("%d\n",ans);
     }
 }

思路积累：

• 分块不一定都是暴力。
• 预处理所有询问。

T2:折射

dp式子都不一样。我说我的按y排序后的n²空间的了。（因为另一种我并想不到）

先把x离散化，然后按y排序。我是dp[i][j]表示选到了第i个装置，上上个选的装置的x是j。

先考虑n³，emm还是不说了吧真的是暴力啊，加个if就行。

然后其实我们发现是一个区间加，那么前缀和一下复杂度就对了。

这时候的dp[i][x[j]]+=dp[j][x[i]];

我们发现对于每一个i只从前面的i-1项转移而来，那么总的就是1/2 n²而不是n²

那么我们把要加的项目提前放进vector里，把dp数组滚动，就不会再炸空间了。

 #include<cstdio>
 #include<algorithm>
 #include<vector>
 using namespace std;
 vector<int>num[];
 vector<short>pos[];
 int mod(int p){return p>=?p-:p;}
 struct ps{int x,y;friend bool operator<(ps a,ps b){return a.y>b.y;}}p[];
 bool com(ps a,ps b){return a.x<b.x;}
 int dp[],n,ans;
 int main(){
     scanf("%d",&n);
     for(int i=;i<=n;++i)scanf("%d%d",&p[i].x,&p[i].y);
     sort(p+,p++n,com);
     for(int i=;i<=n;++i)p[i].x=i;
     sort(p+,p++n);
     for(int i=;i<=n;++i) num[i].resize(i-),pos[i].resize(i-);
     for(int i=;i<=n;++i){
         for(int j=;j<=n;++j)dp[j]=;
         for(int j=;j<num[i].size();++j)dp[pos[i][j]]=mod(dp[pos[i][j]]+num[i][j]);
         num[i].clear();pos[i].clear();
         for(int j=;j<i;++j)dp[p[j].x]++;
         for(int j=;j<p[i].x;++j)dp[j]=mod(dp[j]+dp[j-]);
         for(int j=n;j>p[i].x;--j)dp[j]=mod(dp[j]+dp[j+]);
         ans=mod(ans+mod(dp[p[i].x-]+dp[p[i].x+]));
         for(int j=i+;j<=n;++j)num[j][j-i-]=dp[p[j].x],pos[j][j-i-]=p[i].x;
     }
     printf("%d\n",mod(ans+n));
 }

思路积累：

• resize
• 前缀和

T3：画作

大神题。

有一个结论就是最优方案中一定有一种的所有操作的交集不是空集。

那么我们只要抓住某一个点向四周扩展，把所有同色区域取反就好了，知道全部区域都是白色。

复杂度$O(n^5)$

同色区域被多次扩展了，我们只在乎最远的点几次能扩展几次。

同色连0边，异色连1边，跑bfs即可。

因为懒得打双端队列所以开了100个队列。。。（因为曼哈顿距离不会超过100）

 #include<cstdio>
 #include<iostream>
 using namespace std;
 const int xx[]={,,,-},yy[]={,,-,};
 #define tx X+xx[p]
 #define ty Y+yy[p]
 char s[][];int o[][],n,m,dt[][],qx[][],qy[][],qh[],qt[],ans=;
 int main(){
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;++i)scanf("%s",s[i]+);
     for(int x=;x<=n;++x)for(int y=;y<=m;++y)dt[x][y]=;
     for(int i=;i<=n;++i)for(int j=;j<=m;++j)o[i][j]=s[i][j]-'';
     for(int i=;i<=n;++i)for(int j=;j<=m;++j){
         dt[i][j]=;qh[]=qt[]=;qx[][]=i;qy[][]=j;
         for(int k=;k<=;++k)qh[k]=,qt[k]=;
         for(int k=;k<=;++k)for(int l=qh[k];l<=qt[k];++l){
             int X=qx[k][l],Y=qy[k][l];
             for(int p=;p<=;++p)if(dt[tx][ty]>dt[X][Y]+(o[X][Y]^o[tx][ty]))
                 dt[tx][ty]=dt[X][Y]+(o[X][Y]^o[tx][ty]),
                 qx[dt[tx][ty]][++qt[dt[tx][ty]]]=tx,
                 qy[dt[tx][ty]][qt[dt[tx][ty]]]=ty;
         }
         int mx=;
         for(int x=;x<=n;++x)for(int y=;y<=m;++y)mx=max(mx,dt[x][y]),dt[x][y]=;
         mx+=mx&^o[i][j];
         ans=min(ans,mx);
     }printf("%d\n",ans);
 }