1.CF103D
对输入数据分块
A
C
C
o
d
e
:
AC \ Code:
AC Code:
#include<bits/stdc++.h>
#define maxn 300005
#define S 500
#define LL long long
#define pb push_back
using namespace std;
int n,q;
int w[maxn],a[maxn],b[maxn];
LL ans[maxn],s[maxn];
vector<int>G[S];
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&w[i]);
scanf("%d",&q);
for(int i=1;i<=q;i++){
scanf("%d%d",&a[i],&b[i]);
if(b[i]<S) G[b[i]].pb(i);
else{
for(int j=a[i];j<=n;j+=b[i])
ans[i] += w[j];
}
}
for(int i=1;i<S;i++)
if(!G[i].empty()){
for(int j=n;j>=1;j--)
s[j] = w[j] + (j+i<=n ? s[j+i] : 0);
for(int u:G[i])
ans[u] = s[a[u]];
}
for(int i=1;i<=q;i++)
printf("%I64d\n",ans[i]);
}
2.BZOJ 2054
长度为
n
n
n的线段,
m
m
m次区间颜色覆盖,求最后的线段颜色的样子。
因为只有一次询问,所以反过来每次加入一个区间然后删除这个区间,用并查集维护每一个点下一个没有被染色的地方,然后暴力跳就行了,每个点只会被删除一次。
O
(
n
α
)
O(n\alpha)
O(nα)
被告知只使用路径压缩的并查集是
O
(
log
n
)
O(\log n)
O(logn)的,得加上按秩合并
3.CF 480 E
好题
A
C
C
o
d
e
AC\ Code
AC Code
#include<bits/stdc++.h>
#define maxn 2005
using namespace std;
int n,m,q,ans;
char s[maxn][maxn];
int x[maxn],y[maxn],ret[maxn],f[maxn][maxn],l[maxn][maxn],r[maxn][maxn];
bool check(int L,int y){
deque<int>q1,q2;
for(int i=0;i<n;i++){
for(;!q1.empty() && l[q1.back()][y]<=l[i][y];q1.pop_back());
for(;!q2.empty() && r[q2.back()][y]>=r[i][y];q2.pop_back());
for(;!q1.empty() && q1.front()<=i-L;q1.pop_front());
for(;!q2.empty() && q2.front()<=i-L;q2.pop_front());
q1.push_back(i) , q2.push_back(i);
if(i>=L-1 && r[q2.front()][y] - l[q1.front()][y] + 1 >= L) return 1;
}
return 0;
}
int main(){
scanf("%d%d%d",&n,&m,&q);
for(int i=0;i<n;i++)
scanf("%s",s[i]);
for(int i=1;i<=q;i++){
scanf("%d%d",&x[i],&y[i]);x[i]--,y[i]--;
s[x[i]][y[i]] = 'X';
}
for(int i=0;i<n;i++){
for(int j=0;j<m;j++){
if(s[i][j] == '.')
f[i][j] = min(min(i?f[i-1][j]:0 , j?f[i][j-1]:0) , (i&&j)?f[i-1][j-1]:0) + 1,
l[i][j] = j?l[i][j-1]:0;
else
l[i][j] = j+1;
ans = max(ans , f[i][j]);
}
for(int j=m-1;j>=0;j--)
if(s[i][j] == '.')
r[i][j] = j+1<m?r[i][j+1]:m-1;
else
r[i][j] = j-1;
}
ret[q] = ans;
for(int i=q;i>=1;){
s[x[i]][y[i]] = '.';
for(int j=0;j<m;j++)
if(s[x[i]][j] == '.')
l[x[i]][j] = j?l[x[i]][j-1]:0;
else
l[x[i]][j] = j+1;
for(int j=m-1;j>=0;j--)
if(s[x[i]][j] == '.')
r[x[i]][j] = j+1<m?r[x[i]][j+1]:m-1;
else
r[x[i]][j] = j-1;
for(;check(ans+1,y[i]);ans++);
ret[--i] = ans;
}
for(int i=1;i<=q;i++) printf("%d\n",ret[i]);
}
4.ZJOI 2016 旅行者
分治最短路。
A C C o d e AC\ Code AC Code
#include<bits/stdc++.h>
#define maxn 100005
#define pii pair<int,int>
#define mp make_pair
using namespace std;
int n,m,Q,u1[maxn],u2[maxn],v1[maxn],v2[maxn],c[maxn],q1[maxn],q2[maxn],ans[maxn];
int info[maxn],Prev[maxn<<1],to[maxn<<1],cst[maxn<<1],cnt_e;
inline int id(int a,int b){ return (a-1)*m+b; }
void Node(int u,int v,int c){ Prev[++cnt_e]=info[u],info[u]=cnt_e,to[cnt_e]=v,cst[cnt_e]=c; }
int vis[maxn],tim,dis[maxn];
priority_queue<pii,vector<pii >,greater<pii > >q;
void Solve(int l1,int l2,int r1,int r2,int l,int r){
if(l1>r1 || l2>r2 || l>r) return;
int L1=0,L2=0;
if(r1-l1>r2-l2){
int mid = (l1 + r1) >> 1;
for(int i=l2;i<=r2;i++){
for(int a=l1;a<=r1;a++)
for(int b=l2;b<=r2;b++)
{
if(a!=mid || b!=i)
dis[id(a,b)] = 0x3f3f3f3f;
else
dis[id(a,b)] = 0,
q.push(mp(dis[id(a,b)],id(a,b)));
vis[id(a,b)] = 1;
}
for(;!q.empty();){
pii now = q.top();q.pop();
if(!vis[now.second]) continue;
vis[now.second] = 0;
for(int j=info[now.second];j;j=Prev[j])
if(vis[to[j]] && dis[to[j]] > dis[now.second] + cst[j]){
dis[to[j]] = dis[now.second] + cst[j];
q.push(mp(dis[to[j]] , to[j]));
}
}
for(int j=l;j<=r;j++)
ans[c[j]] = min(ans[c[j]] , dis[id(u1[c[j]],u2[c[j]])] + dis[id(v1[c[j]],v2[c[j]])]);
}
for(int i=l;i<=r;i++)
if(u1[c[i]]<mid && v1[c[i]]<mid) q1[L1++] = c[i];
else if(u1[c[i]]>mid && v1[c[i]]>mid) q2[L2++] = c[i];
for(int i=l;i<l+L1;i++) c[i] = q1[i-l];
for(int i=l+L1;i<l+L1+L2;i++) c[i] = q2[i-l-L1];
Solve(l1,l2,mid-1,r2,l,l+L1-1),Solve(mid+1,l2,r1,r2,l+L1,l+L1+L2-1);
}
else{
int mid = (l2 + r2) >> 1;
for(int i=l1;i<=r1;i++){
for(int a=l1;a<=r1;a++)
for(int b=l2;b<=r2;b++){
if(a!=i || b!=mid)
dis[id(a,b)] = 0x3f3f3f3f;
else
dis[id(a,b)] = 0,
q.push(mp(dis[id(a,b)],id(a,b)));
vis[id(a,b)] = 1;
}
for(;!q.empty();){
pii now = q.top();q.pop();
if(!vis[now.second]) continue;
vis[now.second] = 0;
for(int j=info[now.second];j;j=Prev[j])
if(vis[to[j]] && dis[to[j]] > dis[now.second] + cst[j]){
dis[to[j]] = dis[now.second] + cst[j];
q.push(mp(dis[to[j]] , to[j]));
}
}
for(int j=l;j<=r;j++)
ans[c[j]] = min(ans[c[j]] , dis[id(u1[c[j]],u2[c[j]])] + dis[id(v1[c[j]],v2[c[j]])]);
}
for(int i=l;i<=r;i++)
if(u2[c[i]]<mid && v2[c[i]]<mid) q1[L1++] = c[i];
else if(u2[c[i]]>mid && v2[c[i]]>mid) q2[L2++] = c[i];
for(int i=l;i<l+L1;i++) c[i] = q1[i-l];
for(int i=l+L1;i<l+L1+L2;i++) c[i] = q2[i-l-L1];
Solve(l1,l2,r1,mid-1,l,l+L1-1),Solve(l1,mid+1,r1,r2,l+L1,l+L1+L2-1);
}
}
int main(){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
for(int j=1,a;j<m;j++){
scanf("%d",&a);
Node(id(i,j),id(i,j+1),a);
Node(id(i,j+1),id(i,j),a);
}
for(int i=1,a;i<n;i++)
for(int j=1;j<=m;j++){
scanf("%d",&a);
Node(id(i,j),id(i+1,j),a);
Node(id(i+1,j),id(i,j),a);
}
scanf("%d",&Q);
for(int i=1;i<=Q;i++){
scanf("%d%d%d%d",&u1[i],&u2[i],&v1[i],&v2[i]);
c[i] = i;
}
memset(ans,0x3f,sizeof ans);
Solve(1,1,n,m,1,Q);
for(int i=1;i<=Q;i++)
printf("%d\n",ans[i]);
}