【Codeforces】 CF79D Password

题目链接

CF方向
Luogu方向

题目解法

看到区间异或,一个经典的套路是做差分,我们即在 l l l 处异或一次,在 r + 1 r+1 r+1 处异或一次,然后前缀和起来
于是我们可以将问题转化成:有一个序列初始全 0 0 0,每次可以把相隔 a i a_i ai 的数都 ⊕ 1 \oplus 1 1,求最少将其变成一个状态的步数
考虑 k k k 的范围很小,所以为 1 1 1 的地方一共只有 2 k 2k 2k
这里有一个非常重要的 t r i c k trick trick:在异或操作中,如果需要把 x , y x,y x,y 同时异或 1 1 1,其他不变,每次可以同时修改相隔 a i a_i ai 的位置的异或值,那么这个问题等价于建出图来从 x x x y y y 的最短路
然后发现直接状压跑最短路即可,时间复杂度 O ( 2 k k 2 ) O(2^kk^2) O(2kk2)
不难优化成 O ( 2 k k ) O(2^kk) O(2kk),但我直接 998 m s 998ms 998ms O ( 2 k k 2 ) O(2^kk^2) O(2kk2) 的做法艹过去了,就懒得改了
O ( 2 k k 2 ) O(2^kk^2) O(2kk2) 的代码:

#include 
using namespace std;
const int N=10100,M=2000100;
int n,m,k,a[110],x[30],dis[N];
int f[(1<<20)+100],D[30][30];
int e[M],ne[M],h[N],idx;
inline int read(){
    int FF=0,RR=1;
    char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if(ch=='-') RR=-1;
    for(;isdigit(ch);ch=getchar()) FF=(FF<<1)+(FF<<3)+ch-48;
    return FF*RR;
}
queue<int> que;
void bfs(int S){
    memset(dis,0x3f,sizeof(dis));
    que.push(S),dis[S]=0;
    while(!que.empty()){
        int u=que.front();que.pop();
        for(int i=h[u];~i;i=ne[i]) if(dis[u]+1<dis[e[i]])
            dis[e[i]]=dis[u]+1,que.push(e[i]);
    }
}
void add(int x,int y){ e[idx]=y,ne[idx]=h[x],h[x]=idx++;}
int main(){
    n=read(),k=read(),m=read();
    for(int i=0;i<k;i++) x[i]=read();
    for(int i=1;i<=m;i++) a[i]=read();
    for(int i=0;i<k;i++) x[i+k]=x[i]+1;
    memset(h,-1,sizeof(h));
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n-a[i]+1;j++)
            add(j,j+a[i]),add(j+a[i],j);
    for(int i=0;i<k<<1;i++){
        bfs(x[i]);
        for(int j=0;j<k<<1;j++) D[i][j]=dis[x[j]];
    }
    memset(f,0x3f,sizeof(f));
    f[0]=0;
    for(int S=0;S<1<<(k<<1);S++)
        for(int i=0;i<k<<1;i++) if(S>>i&1)
            for(int j=0;j<k<<1;j++) if(S>>j&1) if(i!=j)
                f[S]=min(f[S],f[S^(1<<i)^(1<<j)]+D[i][j]);
    printf("%d\n",f[(1<<(k<<1))-1]>1e9?-1:f[(1<<(k<<1))-1]);
    fprintf(stderr,"%d ms\n",int(1e3*clock()/CLOCKS_PER_SEC));
    return 0;
}

你可能感兴趣的:(Codeforces,算法)