T1
解析:
单调栈。
分析题意,我们首先发现一个合法区间的因数一定是这个区间中最小的。其次发现本质是求以每个数作为最小质因子,最大的区间,最后答案就是这些区间长度的最大值。
用单调栈计算每个点的左右端点,具体来说如果当前这个数模栈顶的元素不为
0
0
0,即不是栈顶的元素不是这个数的因数,那么就弹出栈顶知道合法,并且在弹出的过程中更新弹出元素的左端点/右端点。详见代码。
代码:
#include <bits/stdc++.h>
using namespace std;
const int Max=500005;
int n,m,top,ans,sum,tot;
int v[Max],l[Max],r[Max],a[Max],d[Max];
struct shu{int id,x;}p[Max];
inline int get_int()
{
int x=0;char c;
for(c=getchar();!isdigit(c);c=getchar());
for(;isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+(c^48);
return x;
}
inline void print(int x)
{
if(x>9) print(x/10);
putchar('0'+x%10);
}
int main()
{
n=get_int();
for(int i=1;i<=n;i++) a[i]=get_int();
for(int i=1;i<=n;i++)
{
while(top&&a[i]%p[top].x) r[p[top--].id]=i-1;
p[++top].x=a[i],p[top].id=i;
}
while(top) r[p[top--].id]=n;
for(int i=n;i;i--)
{
while(top&&a[i]%p[top].x) l[p[top--].id]=i+1;
p[++top].x=a[i],p[top].id=i;
}
while(top) l[p[top--].id]=1;
for(int i=1;i<=n;i++) ans=max(ans,r[i]-l[i]);
for(int i=1;i<=n;i++)
if(r[i]-l[i]==ans&&!v[l[i]]) sum++,v[l[i]]=1,d[++tot]=l[i]; //v去重
print(sum),putchar(' '),print(ans),putchar('\n');
for(int i=1;i<=tot;i++) print(d[i]),putchar(' ');
return 0;
}
T2
解析:
最小生成树?
先来看幅图:
我们发现当直径不断增加的时候,能通过的缝隙会不断减少,当缝隙减少到将通道全部隔断的时候就不能通过了,这时候加入的边即为答案。
所以可以用类似最小生成树的思想,点两两之间连边,边长位距离。特别的,对于点与两界限的情况,建立两个超级点,分别表示上界与下界,进行连边。做最小生成树,当发现两超级点连通时,说明有一个地方已经全部隔断,答案即为最后加入的这条边。
代码:
#include <bits/stdc++.h>
using namespace std;
const int Max=505;
int n,fa[Max],tot;
double ans,m;
struct shu{double x,y;}p[Max];
struct kru{double x,y,z;}e[Max*Max+Max*2];
inline char getc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int get_int()
{
int x=0,f=1;char c;
for(c=getc();(!isdigit(c))&&(c!='-');c=getc());
if(c=='-') f=-1,c=getc();
for(;isdigit(c);c=getc()) x=(x<<3)+(x<<1)+(c^48);
return x*f;
}
inline int get(int v){return fa[v]==v?v:fa[v]=get(fa[v]);}
inline bool comp(const kru &a,const kru &b){return a.z<b.z;}
inline double calc(int i,int j){return sqrt((p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y));}
int main()
{
n=get_int(),m=get_int(),fa[n+1]=n+1;
for(int i=1;i<=n;i++)
{
fa[i]=i,p[i].x=get_int(),p[i].y=get_int();
e[++tot].x=0,e[tot].y=i,e[tot].z=p[i].y;
e[++tot].x=n+1,e[tot].y=i,e[tot].z=m-p[i].y;
}
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++)
e[++tot].x=i,e[tot].y=j,e[tot].z=calc(i,j);
sort(e+1,e+tot+1,comp);
for(int i=1;i<=tot;i++)
{
if(get(e[i].x)==get(e[i].y)) continue;
fa[get(e[i].y)]=e[i].x;
if(get(n+1)==get(0)) {ans=e[i].z;break;}
}
printf("%.3f",ans);
return 0;
}
T3
解析:
斜率优化。
这个坑以后再填吧
T
T
TT
TT。当然玄学广搜也是能过的。。。
代码(广搜):
#include<iostream>
#include<cstdio>
#include<queue>
#include<cmath>
#include<cstring>
#include<cctype>
using namespace std;
inline int read(){
int x=0,f=1;char c=getchar();
while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f==1?x:-x;
}
int n,m,a[1005][1005],e[1000005],f[1000005],t=0,ans[1005][1005];
struct node{
int x,y,r;
};
queue<node>q;
int dx[8]={-1,0,1,0,-1,1,1,-1},dy[8]={0,1,0,-1,1,1,-1,-1};//八个方向!!!!!
void work1(){
char c;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++){
c=getchar();
while(c!='0'&&c!='1')c=getchar();
if(c^48){
e[++t]=i;f[t]=j;
q.push(node{i,j,t});
}
}
while(!q.empty()){
node u=q.front();q.pop();
for(int i=0;i<8;i++){
int xx=u.x+dx[i],yy=u.y+dy[i];
if(xx<1||xx>n||yy<1||yy>m)continue;
int x1=(xx-e[u.r])*(xx-e[u.r])+(yy-f[u.r])*(yy-f[u.r])+1;
if(!ans[xx][yy]||ans[xx][yy]>x1){
ans[xx][yy]=x1;
q.push(node{xx,yy,u.r});
}
}
}
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++)
printf("%d ",ans[i][j]-1);
printf("\n");
}
return;
}
int main(){
n=read();m=read();
work1();
return 0;
}