#include <bits/stdc++.h> using namespace std; typedef long long ll; const int N=(1<<21); const ll INF=4e18; /*struct Tree { int mn[2*N]; void update(int x, int a) { x+=N; mn[x]=a; while(x>0) { x/=2; mn[x]=min(mn[2*x], mn[2*x+1]); } } int query(int xa, int xb) { xa+=N; xb+=N; if(xa==xb) return mn[xa]; int res=min(mn[xa], mn[xb]); while(xa/2!=xb/2) { if(xa%2==0) res=min(res, mn[xa+1]); if(xb%2==1) res=min(res, mn[xb-1]); xa/=2; xb/=2; } //printf("?\n"); return res; } }; Tree t;*/ int a[N], d[N], num[N]; char c[N]; vector<int> v[N]; ll result(int n, int m) { for(int i=0; i<m; ++i) { d[i]=c[i*(ll)n%m]; if(d[i]=='W') d[i]=1; else d[i]=-1; num[i*(ll)n%m]=i; } for(int i=1; i<m; ++i) d[i]+=d[i-1]; for(int i=m; i<2*m; ++i) d[i]=d[i-m]+d[m-1]; for(int i=2*m; i>0; --i) d[i]=d[i-1]; d[0]=0; for(int i=0; i<2*m; ++i) if(abs(d[i])<=m) v[d[i]+m].push_back(i); //for(int i=0; i<2*m; ++i) t.update(i, d[i]); int mn=0; for(int i=1; i<=m; ++i) mn=min(d[i], mn); ll res=INF; for(int i=0; i<n; ++i) { //printf("i=%d\n", i); ll lres=0; //Check if bankrupt in a short time if(d[num[i%m]]-a[i]+m>=0&&v[d[num[i%m]]-a[i]+m].size()>0&&v[d[num[i%m]]-a[i]+m].back()>num[i%m]) { lres=*lower_bound(v[d[num[i%m]]-a[i]+m].begin(), v[d[num[i%m]]-a[i]+m].end(), num[i%m])-num[i%m]; } else { //printf("Here!\n"); lres=m-num[i%m]; //if(i==2) printf("i=%d, a[i]=%d, n=%d, m=%d, lres=%lld\n", i, a[i], n, m, lres); a[i]+=d[m]-d[num[i%m]]; if(d[m]>=0) lres=INF; else { if(a[i]<=-mn) { lres+=v[m-a[i]][0]; } else { lres+=m*(ll)((a[i]+mn)/(-d[m])); a[i]=(a[i]+mn-d[m]*(ll)m)%(-d[m])-mn; if(a[i]>-mn) { a[i]-=-d[m]; lres+=m; } lres+=v[m-a[i]][0]; } } } if(lres!=INF) lres=max(lres-1, 0LL)*n+i+1; //if(lres==2) printf("i=%d, a[i]=%d, n=%d, m=%d\n", i, a[i], n, m); res=min(res, lres); } for(int i=0; i<n; ++i) a[i]=0; for(int i=0; i<=m; ++i) d[i]=num[i]=c[i]=0; for(int i=0; i<=2*m; ++i) v[i].clear(); return res; } int ga[N]; char gc[N]; int main() { //freopen("tests/big0.in", "r", stdin); int n, m; scanf("%d", &n); for(int i=0; i<n; ++i) scanf("%d", &ga[i]); scanf("%d\n%s", &m, gc); int gcd=__gcd(n, m); ll res=INF; for(int k=0; k<gcd; ++k) { for(int i=0; i<n/gcd; ++i) a[i]=ga[k+gcd*i]; for(int i=0; i<m/gcd; ++i) c[i]=gc[k+gcd*i]; ll tmp=result(n/gcd, m/gcd); if(tmp!=INF) res=min(res, gcd*(tmp-1)+k+1); } if(res<INF) printf("%lld\n", res); else printf("-1\n"); return 0; }
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 | #include <bits/stdc++.h> using namespace std; typedef long long ll; const int N=(1<<21); const ll INF=4e18; /*struct Tree { int mn[2*N]; void update(int x, int a) { x+=N; mn[x]=a; while(x>0) { x/=2; mn[x]=min(mn[2*x], mn[2*x+1]); } } int query(int xa, int xb) { xa+=N; xb+=N; if(xa==xb) return mn[xa]; int res=min(mn[xa], mn[xb]); while(xa/2!=xb/2) { if(xa%2==0) res=min(res, mn[xa+1]); if(xb%2==1) res=min(res, mn[xb-1]); xa/=2; xb/=2; } //printf("?\n"); return res; } }; Tree t;*/ int a[N], d[N], num[N]; char c[N]; vector<int> v[N]; ll result(int n, int m) { for(int i=0; i<m; ++i) { d[i]=c[i*(ll)n%m]; if(d[i]=='W') d[i]=1; else d[i]=-1; num[i*(ll)n%m]=i; } for(int i=1; i<m; ++i) d[i]+=d[i-1]; for(int i=m; i<2*m; ++i) d[i]=d[i-m]+d[m-1]; for(int i=2*m; i>0; --i) d[i]=d[i-1]; d[0]=0; for(int i=0; i<2*m; ++i) if(abs(d[i])<=m) v[d[i]+m].push_back(i); //for(int i=0; i<2*m; ++i) t.update(i, d[i]); int mn=0; for(int i=1; i<=m; ++i) mn=min(d[i], mn); ll res=INF; for(int i=0; i<n; ++i) { //printf("i=%d\n", i); ll lres=0; //Check if bankrupt in a short time if(d[num[i%m]]-a[i]+m>=0&&v[d[num[i%m]]-a[i]+m].size()>0&&v[d[num[i%m]]-a[i]+m].back()>num[i%m]) { lres=*lower_bound(v[d[num[i%m]]-a[i]+m].begin(), v[d[num[i%m]]-a[i]+m].end(), num[i%m])-num[i%m]; } else { //printf("Here!\n"); lres=m-num[i%m]; //if(i==2) printf("i=%d, a[i]=%d, n=%d, m=%d, lres=%lld\n", i, a[i], n, m, lres); a[i]+=d[m]-d[num[i%m]]; if(d[m]>=0) lres=INF; else { if(a[i]<=-mn) { lres+=v[m-a[i]][0]; } else { lres+=m*(ll)((a[i]+mn)/(-d[m])); a[i]=(a[i]+mn-d[m]*(ll)m)%(-d[m])-mn; if(a[i]>-mn) { a[i]-=-d[m]; lres+=m; } lres+=v[m-a[i]][0]; } } } if(lres!=INF) lres=max(lres-1, 0LL)*n+i+1; //if(lres==2) printf("i=%d, a[i]=%d, n=%d, m=%d\n", i, a[i], n, m); res=min(res, lres); } for(int i=0; i<n; ++i) a[i]=0; for(int i=0; i<=m; ++i) d[i]=num[i]=c[i]=0; for(int i=0; i<=2*m; ++i) v[i].clear(); return res; } int ga[N]; char gc[N]; int main() { //freopen("tests/big0.in", "r", stdin); int n, m; scanf("%d", &n); for(int i=0; i<n; ++i) scanf("%d", &ga[i]); scanf("%d\n%s", &m, gc); int gcd=__gcd(n, m); ll res=INF; for(int k=0; k<gcd; ++k) { for(int i=0; i<n/gcd; ++i) a[i]=ga[k+gcd*i]; for(int i=0; i<m/gcd; ++i) c[i]=gc[k+gcd*i]; ll tmp=result(n/gcd, m/gcd); if(tmp!=INF) res=min(res, gcd*(tmp-1)+k+1); } if(res<INF) printf("%lld\n", res); else printf("-1\n"); return 0; } |