[考试反思]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...)

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

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

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

#include<cstdio>
#include<iostream>
using namespace std;
#define S 1234
#define N 100001
int lst[100005],n,m,x[130005],Ans[101][130005],bl[130005];
int main(){
    scanf("%d%d",&n,&m);
    const int bs=(n-1)/S+1;
    for(int i=1;i<=n;++i)scanf("%d",&x[i]);
    for(int i=1;i<=bs;++i){
        for(int j=1;j<=N;++j)lst[j]=0;
        for(int j=(i-1)*S+1;j<=i*S&&j<=n;++j)lst[x[j]]=x[j],bl[j]=i;
        for(int j=1;j<=N;++j)lst[j]=max(lst[j-1],lst[j]);
        for(int j=2;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-1]%j);
        }
    }
    for(int t=1,l,r,q;t<=m;++t){
        scanf("%d%d%d",&l,&r,&q);int ans=0;
        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]-1)*S+1;i<=r;++i)ans=max(ans,x[i]%q);
            for(int i=bl[l]+1;i<bl[r];++i)ans=max(ans,Ans[i][q]);
        }printf("%d
",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[6001];
vector<short>pos[6001];
int mod(int p){return p>=1000000007?p-1000000007:p;}
struct ps{int x,y;friend bool operator<(ps a,ps b){return a.y>b.y;}}p[6005];
bool com(ps a,ps b){return a.x<b.x;}
int dp[6001],n,ans;
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;++i)scanf("%d%d",&p[i].x,&p[i].y);
    sort(p+1,p+1+n,com);
    for(int i=1;i<=n;++i)p[i].x=i;
    sort(p+1,p+1+n);
    for(int i=1;i<=n;++i) num[i].resize(i-1),pos[i].resize(i-1);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j)dp[j]=0;
        for(int j=0;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=1;j<i;++j)dp[p[j].x]++;
        for(int j=1;j<p[i].x;++j)dp[j]=mod(dp[j]+dp[j-1]);
        for(int j=n;j>p[i].x;--j)dp[j]=mod(dp[j]+dp[j+1]);
        ans=mod(ans+mod(dp[p[i].x-1]+dp[p[i].x+1]));
        for(int j=i+1;j<=n;++j)num[j][j-i-1]=dp[p[j].x],pos[j][j-i-1]=p[i].x;
    }
    printf("%d
",mod(ans+n));
}

思路积累：

• resize
• 前缀和

T3：画作

大神题。

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

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

复杂度$O(n^5)$

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

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

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

#include<cstdio>
#include<iostream>
using namespace std;
const int xx[]={1,0,0,-1},yy[]={0,1,-1,0};
#define tx X+xx[p]
#define ty Y+yy[p]
char s[55][55];int o[55][55],n,m,dt[55][55],qx[105][2555],qy[105][2555],qh[105],qt[105],ans=123;
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;++i)scanf("%s",s[i]+1);
    for(int x=1;x<=n;++x)for(int y=1;y<=m;++y)dt[x][y]=1000;
    for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)o[i][j]=s[i][j]-'0';
    for(int i=1;i<=n;++i)for(int j=1;j<=m;++j){
        dt[i][j]=0;qh[0]=qt[0]=1;qx[0][1]=i;qy[0][1]=j;
        for(int k=1;k<=100;++k)qh[k]=1,qt[k]=0;
        for(int k=0;k<=100;++k)for(int l=qh[k];l<=qt[k];++l){
            int X=qx[k][l],Y=qy[k][l];
            for(int p=0;p<=3;++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=0;
        for(int x=1;x<=n;++x)for(int y=1;y<=m;++y)mx=max(mx,dt[x][y]),dt[x][y]=1000;
        mx+=mx&1^o[i][j];
        ans=min(ans,mx);
    }printf("%d
",ans);
}