1. 程式人生 > >LOJ#2799. 「CCC 2016」生命之環

LOJ#2799. 「CCC 2016」生命之環

題意

給你一個 \(n\)\(\rm 01\) 組成的環,每次操作之後每個位置為1當且僅當他的左右恰好有1個1.輸出進行 \(T\) 次操作之後的環。

\(n\leq 10^5, T\leq 10^{15}\).

分析

  • 通過1~4步之內模擬可以得到結論:一個位置能夠在 \(2^k\) 的操作之後為1當且僅當他的往左往右的 \(2^k\) 個位置的異或值為1.

  • 將數字拆成若干個 \(2^k\) 進行操作即可。

  • 總時間複雜度為 \(O(nlogT)\)

程式碼

#include<bits/stdc++.h>
using namespace std;
#define go(u) for(int i=head[u],v=e[i].to;i;i=e[i].lst,v=e[i].to)
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define pb push_back
typedef long long LL;
inline int gi(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-48;ch=getchar();}
    return x*f;
}
template<typename T>inline bool Max(T &a,T b){return a<b?a=b,1:0;}
template<typename T>inline bool Min(T &a,T b){return b<a?a=b,1:0;}
const int N=1e5 +7;
LL n,T;
LL a[N],b[N];
char s[N];
int main(){
    scanf("%lld%lld",&n,&T);
    scanf("%s",s);
    rep(i,0,n-1) a[i]=s[i]-'0';
    for(int k=61;~k;--k)if(T&(1ll<<k)){
        memset(b,0,sizeof b);
        rep(i,0,n-1){
            b[((i+(1ll<<k))%n+n)%n]^=a[i];
            b[((i-(1ll<<k))%n+n)%n]^=a[i];
        }
        memcpy(a,b,sizeof a);
    }
    rep(i,0,n-1) printf("%d",a[i]);
    puts("");
    return 0;
}