3/28 并查集+最短路+爆搜dfs(记忆化,方向控制)

P2658 汽车拉力比赛
这题的数据好像有问题，思路明明一样，代码的实现稍微有点区别，就会有样例报re。真的浪费了好长时间

#include <bits/stdc++.h>

using namespace std;
int h[501][501],f[250001],a[250001],cnt;
int ans,l,r,n,m,mid;
int dx[4]={-1,0,1,0};
int dy[4]={0,1,0,-1};
int id(int x,int y)
{
    return (x-1)*m+y;
}
int r_find(int r)
{
    if(r==f[r]) return f[r];
    f[r]=r_find(f[r]);
    return f[r];
}
void add(int x,int y){f[r_find(x)]=r_find(y);return;}

bool check(int high)
{
	for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            f[id(i,j)]=id(i,j);
	for(int i=1;i<=n;i++)
	 for(int j=1;j<=m;j++)
	  for(int k=0;k<4;k++)
	  {
	  	int nx=i+dx[k],ny=j+dy[k];
	 	if(nx<1||ny<1||nx>n||ny>m)
            continue;
	 	if(abs(h[i][j]-h[nx][ny])>high)
            continue;
	 	add(id(i,j),id(nx,ny));
      }
	for(int i=1;i<cnt;i++)
	if(r_find(a[i])!=r_find(a[i+1]))
        return 0;
	return 1;
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            scanf("%d",&h[i][j]);
            if(h[i][j]>r) r=h[i][j];
        }
	for(int i=1;i<=n;i++)
	for(int j=1,t;j<=m;j++)
	{
	    scanf("%d",&t);if(t)
        a[++cnt]=id(i,j);
    }
	while(l<=r)
	{
		mid=(l+r)>>1;
		if(check(mid))
		{
			ans=mid;
			r=mid-1;
		}
		else
            l=mid+1;
	}
	printf("%d",ans);
}

P2937 [USACO09JAN]Laserphones S
开一个三维数组控制方向。
记忆化：如果已知到达这个方向的点的最小次数，则不必在进行搜索了。
方向的控制:刚开始的方向特殊处理，无需累加进去。

#include <bits/stdc++.h>

using namespace std;
const int inf=0x3f3f3f3f;
const int N=105;
int n,m,f[N][N][5],ans=inf;
char mp[N][N];
bool vis[N][N];
int dx[4]={-1,0,1,0};
int dy[4]={0,1,0,-1};
void dfs(int x,int y,int num,int dir) //dir表示方向
{
    if(dir!=-1&&f[x][y][dir]<=num)
        return;
    if(dir!=-1) f[x][y][dir]=num;
    if(mp[x][y]=='C')
    {
        ans=min(ans,num);return;
    }
    if(dir!=-1)
    {
        int nx=x+dx[dir],ny=y+dy[dir];
        if(!vis[nx][ny]&&mp[nx][ny]!='*'&&nx>=1&&nx<=n&&ny>=1&&ny<=m)
        {
            vis[nx][ny]=1;
            dfs(nx,ny,num,dir);
            vis[nx][ny]=0;
        }

    }
    for(int i=0;i<4;i++)
    {
        if(i==dir) continue;
        int nx=x+dx[i],ny=y+dy[i];
        if(!vis[nx][ny]&&mp[nx][ny]!='*'&&nx>=1&&nx<=n&&ny>=1&&ny<=m)
        {
            vis[nx][ny]=1;
            dfs(nx,ny,num+(dir!=-1),i);
            vis[nx][ny]=0;
        }
    }

}
int main()
{
    memset(f,inf,sizeof(f));
	scanf("%d%d",&m,&n);
	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++)
        cin>>mp[i][j];
	for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
        {
            if(mp[i][j]=='C')
            {
                vis[i][j]=1;mp[i][j]='*';
                dfs(i,j,0,-1);
                cout<<ans<<endl;
                return 0;
            }
        }

	return 0;
}

G. Gliding
不知道为啥过不了，可能被卡精度了。

#include <bits/stdc++.h>
using namespace std;
const double inf=1e20;
const int N=5e4+6;
struct point
{
    double x,y;double v;
}a[N];
struct Edge
{
    int to,nxt;double dis;
}e[N];
int head[N],n,cnt,sx,sy,tx,ty;
double vp,vf,vh;
double dist[N];
bool vis[N];
void add(int from,int to,double dis)
{
    e[++cnt].to=to;
    e[cnt].dis=dis;
    e[cnt].nxt=head[from];
    head[from]=cnt;
}
bool cmp(point p1,point p2)
{
    return p1.v<p2.v;
}
struct node
{
    int pos;double dis;
    bool operator <(const node &x)const
    {
        return x.dis<dis;
    }
};
void dijistra()
{
    for(int i=0;i<=n;i++)
        dist[i]=inf,vis[i]=0;
    priority_queue<node>q;
    dist[n]=0;
    q.push(node{n,0});
    while(!q.empty())
    {
        node cur=q.top();q.pop();
        int u=cur.pos;
        if(vis[u]) continue;
        vis[u]=1;
        for(int i=head[u];~i;i=e[i].nxt)
        {
            int v=e[i].to;
            if(dist[v]>dist[u]+e[i].dis)
            {
                dist[v]=dist[u]+e[i].dis;
                q.push(node{v,dist[v]});
            }
        }
    }
}
void init()
{
    cnt=0;
    memset(head,-1,sizeof(head));
    memset(a,0,sizeof(a));
    memset(e,0,sizeof(e));
}

int main()
{
    ios::sync_with_stdio(0); cin.tie(0);
    cout << fixed << setprecision(11);
    int t;cin>>t;
    while(t--)
    {
        init();
        cin>>sx>>sy>>tx>>ty;
        cin>>vf>>vp>>vh>>n;
        for(int i=0;i<=n;i++)
        {
            cin>>a[i].x>>a[i].y>>a[i].v;
        }
        n++;
        a[n].x=tx,a[n].y=ty,a[n].v=inf;
        sort(a,a+n,cmp);
        for(int i=n;i>=1;i--)
        {
            if(a[i].v<=vp) break;
            for(int j=i-1;j>=0;j--)
            {
                if(a[j].v<=vp) break;
                double dd=sqrt((double)(a[i].x-a[j].x)*(a[i].x-a[j].x)*1.0+(double)(a[i].y-a[j].y)*(a[i].y-a[j].y)*1.0)/(double)(vh); //时间
                dd+=dd*vp/(a[j].v-vp);
                //cout<<"ddd "<<dd<<" "<<i<<" "<<j<<endl;
                add(i,j,dd);
            }
        }
        dijistra();
        int g;
        for(int i=0;i<=n;i++)
        {
            if(a[i].x==sx&&a[i].y==sy)
            {
                g=i;break;
            }
        }
        cout<<dist[g]<<endl;
    }
	return 0;
}