【NOIP模拟】T1+T2+T3

T1

解析：

  单调栈。
  分析题意，我们首先发现一个合法区间的因数一定是这个区间中最小的。其次发现本质是求以每个数作为最小质因子，最大的区间，最后答案就是这些区间长度的最大值。
  用单调栈计算每个点的左右端点，具体来说如果当前这个数模栈顶的元素不为$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

解析：

  斜率优化。
  这个坑以后再填吧$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;
}