T 1 T_1 T1——mecho(3070)
Description:
有一个
n
⋅
m
n\cdot m
n⋅m的字符矩阵,其中
′
M
′
'M'
′M′为起点,
′
D
′
'D'
′D′为终点,
′
T
′
'T'
′T′为障碍,
′
H
′
'H'
′H′为敌人起点,敌人会在每一时刻扩散到周围4格,而你没一时刻最多可以移动
S
S
S格,求最迟可以在哪个时刻出发且安全到终点。
n
,
m
≤
800
n,m\le 800
n,m≤800
Solution:
- 很显然我们可以二分答案,而关键就在check(
实际就是模拟题) - 我们先用 b f s bfs bfs求出敌人到达其它点的最小时间
- 再就是从起点模拟,注意比较符号,题中是人先走,敌人再走
- 这样复杂度为 Θ ( n 2 log ( n 2 ) ) \Theta(n^2\log (n^2)) Θ(n2log(n2))
Code:
#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define SREP(i,f,t) for(int i=(f),i##_end_=(t);i<i##_end_;++i)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;--i)
#define ll long long
template<class T>inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template<class T>inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
template<class T>inline void Rd(T &x){
x=0;char c;
while((c=getchar())<48);
do x=(x<<1)+(x<<3)+(c^48);
while((c=getchar())>47);
}
const int N=802,INF=0x3f3f3f3f;
int dx[]={1,0,-1,0};
int dy[]={0,1,0,-1};
int n,S;
char mp[N][N];
int t[N][N];
struct node{
int x,y;
}M;
queue<node>Q;
int vis[N][N];
bool check_xy(int x,int y){return x<1 || x>n || y<1 || y>n;}
bool check(int tim) {
if(tim>=t[M.x][M.y]) return 0;
while(!Q.empty()) Q.pop();
memset(vis,INF,sizeof(vis));
vis[M.x][M.y]=tim*S;
Q.push((node){M.x,M.y});
while(!Q.empty()) {
node now=Q.front();Q.pop();
if(mp[now.x][now.y]=='D') return 1;
REP(i,0,3) {
int nx=now.x+dx[i];
int ny=now.y+dy[i];
if(check_xy(nx,ny))continue;
if(mp[nx][ny]=='T')continue;
if(vis[nx][ny]==INF && (vis[now.x][now.y]+1)/S<t[nx][ny]) {
vis[nx][ny]=vis[now.x][now.y]+1;
Q.push((node){nx,ny});
}
}
}
return 0;
}
int main() {
// freopen("mecho.in","r",stdin);
// freopen("mecho.out","w",stdout);
scanf("%d%d",&n,&S);
REP(i,1,n) scanf("%s",mp[i]+1);
memset(t,INF,sizeof(t));
REP(i,1,n) REP(j,1,n) {
if(mp[i][j]=='H') t[i][j]=0,Q.push((node){i,j});;
if(mp[i][j]=='M') M=(node){i,j};
}
while(!Q.empty()){
node now=Q.front();Q.pop();
REP(i,0,3){
int nx=now.x+dx[i];
int ny=now.y+dy[i];
if(check_xy(nx,ny))continue;
if(mp[nx][ny]=='T' || mp[nx][ny]=='D')continue;
if(chkmin(t[nx][ny],t[now.x][now.y]+1)) Q.push((node){nx,ny});
}
}
int L=0,R=n*n,ans=-1;
while(L<=R){
int mid=(L+R)>>1;
if(check(mid)) ans=mid,L=mid+1;
else R=mid-1;
}
printf("%d\n",ans);
return 0;
}
T 2 T_2 T2——quality(3071)
Description:
有一个
n
⋅
m
n\cdot m
n⋅m的矩阵,每个格子中的数用
1
1
1~
n
⋅
m
n\cdot m
n⋅m之间不同的数表示,求在一个
h
⋅
w
h \cdot w
h⋅w的字矩阵中的中位数最小,保证
h
⋅
w
h \cdot w
h⋅w为奇数。
n
,
m
≤
1000
n,m\le 1000
n,m≤1000
Solution:
- 因为中位数是固定的一个数,那么我们很容易想到去找一个 h ⋅ w h \cdot w h⋅w的字矩阵中第 K K K大的数最小
- 因而,我们想到可以二分答案,对于check,我们可以二维前缀和出 ( x 1 , y 1 ) (x1,y1) (x1,y1)到 ( x 2 , y 2 ) (x2,y2) (x2,y2)中不大于 m i d mid mid的数的个数。
- 这样复杂度为 Θ ( n m log ( n m ) ) \Theta(nm\log (nm)) Θ(nmlog(nm))
Code:
#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define SREP(i,f,t) for(int i=(f),i##_end_=(t);i<i##_end_;++i)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;--i)
#define ll long long
template<class T>inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template<class T>inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
template<class T>inline void Rd(T &x){
x=0;char c;
while((c=getchar())<48);
do x=(x<<1)+(x<<3)+(c^48);
while((c=getchar())>47);
}
const int N=1002;
int n,m,h,w,Mid;
int A[N][N];
struct p30{
int tmp[N],tot;
void solve(){
int ans=0x3f3f3f3f;
REP(x1,1,n-h+1) REP(y1,1,m-w+1){
tot=0;
SREP(x2,x1,x1+h) SREP(y2,y1,y1+w) tmp[++tot]=A[x2][y2];
sort(tmp+1,tmp+1+tot);
chkmin(ans,tmp[Mid]);
}
printf("%d\n",ans);
}
}p1;
struct p100{
int sum[N][N];
int Sum(int x1,int y1,int x2,int y2){
return sum[x2][y2]-sum[x2][y1-1]-sum[x1-1][y2]+sum[x1-1][y1-1];
}
bool check(int k){
memset(sum,0,sizeof(sum));
REP(i,1,n) REP(j,1,m) sum[i][j]=(A[i][j]<=k)+sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
REP(i,1,n-h+1) REP(j,1,m-w+1) if(Sum(i,j,i+h-1,j+w-1)>=Mid) return 1;
return 0;
}
void solve(){
int L=1,R=n*m,ans=0;
while(L<=R){
int mid=(L+R)>>1;
if(check(mid)) ans=mid,R=mid-1;
else L=mid+1;
}
printf("%d\n",ans);
}
}p2;
int main(){
// freopen("quality.in","r",stdin);
// freopen("quality.out","w",stdout);
Rd(n),Rd(m),Rd(h),Rd(w);
Mid=(h*w)/2+1;
REP(i,1,n) REP(j,1,m) Rd(A[i][j]);
if(n<=30 && m<=30)p1.solve();
else p2.solve();
return 0;
}
T 3 T_3 T3——salesman(3072)
Description:
有一条类似
x
x
x轴的河上,有
n
n
n个点,每个点有开放的时间
t
i
t_i
ti,坐标
p
i
p_i
pi,权值
w
i
w_i
wi。
初始在坐标
S
S
S处,你可以去其它点后并会来,而逆流而上每个单位的花费为
U
U
U,顺溜而下每个单位的花费为
D
D
D。求最后获得的权值最大为多少。
n
,
t
i
≤
500000
,
S
,
p
i
≤
500001
,
w
i
≤
4000
,
D
≤
U
≤
10
n,t_i \le 500000,S,p_i\le 500001,w_i\le4000,D\le U\le10
n,ti≤500000,S,pi≤500001,wi≤4000,D≤U≤10
Solution:
- 一道比较传统的题
- 首先是一个 Θ ( n 2 ) \Theta(n^2) Θ(n2)的 d p dp dp,这个都没有问题
- 那么发现这个逆流而上和顺流而下可以分类讨论一下,这样就可以通过线段树(BIT)来维护区间的逆流而上的和顺流而下的最大值即可。
- 这样复杂度为 Θ ( n log n ) \Theta(n\log n) Θ(nlogn)
Code:
#include<bits/stdc++.h>
using namespace std;
#define REP(i,f,t) for(int i=(f),i##_end_=(t);i<=i##_end_;++i)
#define SREP(i,f,t) for(int i=(f),i##_end_=(t);i<i##_end_;++i)
#define DREP(i,f,t) for(int i=(f),i##_end_=(t);i>=i##_end_;--i)
#define ll long long
template<class T>inline bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template<class T>inline bool chkmax(T &x,T y){return x<y?x=y,1:0;}
template<class T>inline void Rd(T &x){
x=0;char c;
while((c=getchar())<48);
do x=(x<<1)+(x<<3)+(c^48);
while((c=getchar())>47);
}
const int N=500005;
int n,U,D,S;
struct node{
int t,p,w;
bool operator<(const node &_)const{
return t!=_.t?t<_.t:p<_.p;
}
}A[N];
int dp[N];
struct p45{
int dis(int s,int t){
if(s>=t) return (s-t)*U;
else return (t-s)*D;
}
void solve(){
REP(i,1,n) dp[i]=A[i].w-dis(S,A[i].p);
int ans=0;
REP(i,1,n) {
SREP(j,1,i){
if(A[i].t==A[j].t)continue;
chkmax(dp[i],dp[j]+A[i].w-dis(A[j].p,A[i].p));
}
chkmax(ans,dp[i]-dis(A[i].p,S));
}
printf("%d\n",ans);
}
}p1;
struct p100{
int mx[2][N<<2];
void build(int p,int L,int R) {
if(L==R){
if (L==S) mx[0][p]=D*S,mx[1][p]=-U*S;
else mx[0][p]=mx[1][p]=-0x3f3f3f3f;
return;
}
int mid=(L+R)/2;
build(p<<1,L,mid);
build(p<<1|1,mid+1,R);
mx[0][p]=max(mx[0][p<<1],mx[0][p<<1|1]);
mx[1][p]=max(mx[1][p<<1],mx[1][p<<1|1]);
}
void update(int p,int L,int R,int k) {
if(L==R){
mx[0][p]=dp[k]+D*A[k].p;
mx[1][p]=dp[k]-U*A[k].p;
return;
}
int mid=(L+R)/2;
if(A[k].p<=mid)update(p<<1,L,mid,k);
else update(p<<1|1,mid+1,R,k);
mx[0][p]=max(mx[0][p<<1],mx[0][p<<1|1]);
mx[1][p]=max(mx[1][p<<1],mx[1][p<<1|1]);
}
int query(int f,int p,int L,int R,int l,int r) {
if(l<=L&&R<=r) return mx[f][p];
int mid=(L+R)/2,res=-0x3f3f3f3f;
if(l<=mid) res=query(f,p<<1,L,mid,l,r);
if(r>mid) chkmax(res,query(f,p<<1|1,mid+1,R,l,r));
return res;
}
int f[N],g[N];
void solve(){
A[++n]=(node){A[n-1].t+1,S,0};
int m=S;
REP(i,1,n) chkmax(m,A[i].p);
build(1,1,m);
for(int i=1,j;i<=n;i=j+1){
for(j=i;j<n && A[j+1].t==A[i].t;++j);
REP(k,i,j) dp[k]=max(query(0,1,1,m,1,A[k].p)-D*A[k].p,query(1,1,1,m,A[k].p,m)+U*A[k].p)+A[k].w;
f[i]=dp[i];
REP(k,i+1,j) f[k]=max(dp[k],f[k-1]-D*(A[k].p-A[k-1].p)+A[k].w);
g[j]=dp[j];
DREP(k,j-1,i) g[k]=max(dp[k],g[k+1]-U*(A[k+1].p-A[k].p)+A[k].w);
REP(k,i,j) chkmax(dp[k],max(f[k],g[k])),update(1,1,m,k);
}
printf("%d\n",dp[n]);
}
}p2;
int main(){
// freopen("salesman.in","r",stdin);
// freopen("salesman.out","w",stdout);
Rd(n),Rd(U),Rd(D),Rd(S);
REP(i,1,n) Rd(A[i].t),Rd(A[i].p),Rd(A[i].w);
sort(A+1,A+1+n);
if(n<=5000)p1.solve();
else p2.solve();
return 0;
}