LOJ#2799. 「CCC 2016」生命之環
阿新 • • 發佈:2018-11-23
題意
給你一個 \(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; }