2019 East Central North America

2019 East Central North America

A. A Retribution!

​ 没什么东西，阅读理解题，枚举所有边，选取打标记 (nm + np)

Code：

#include <bits/stdc++.h>
using namespace std;

const int N = 1000 + 10;
struct point{ double x, y; }judge[N], tar[N], store[N];
double dist(double xa,double ya,double xb, double yb)
{
	return sqrt((xa-xb)*(xa-xb)+(ya-yb)*(ya-yb));
}
struct  Node
{
	int x,y;
	double w;
}e[N*N];
bool cmp(Node a, Node b){ return a.w < b.w; }
int visx[N], visy[N];
signed main(){
    ios_base::sync_with_stdio(false), cin.tie(0), cout.tie(0);
    int n, m, p; cin >> n >> m >> p;
    for(int i = 0; i < n; i++) cin >> judge[i].x >> judge[i].y; 
    for(int i = 0; i < m; i++) cin >> tar[i].x >> tar[i].y; 
    for(int i = 0; i < p; i++) cin >> store[i].x >> store[i].y; 
	int cnt=0;
	for(int i=0;i<n;i++)
	{
		for(int j=0; j<m; j++)
		{
			e[++cnt].x = i; e[cnt].y = j; 
			e[cnt].w = dist(judge[i].x, judge[i].y, tar[j].x, tar[j].y); 
		}
	}
	sort(e+1,e+1+cnt,cmp);
	double ans = 0.0;
	for(int i=1;i<=cnt;i++)
	{
		int x = e[i].x, y = e[i].y;
		if(!visx[x] && !visy[y])
		{
			ans += e[i].w;
			visx[x] = 1; visy[y] = 1;
		}
	}
	memset(visx,0,sizeof(visx));
	memset(visy,0,sizeof(visy));
	cnt=0;
	for(int i=0;i<n;i++)
	{
		for(int j=0; j<p; j++)
		{
			e[++cnt].x = i; e[cnt].y = j; 
			e[cnt].w = dist(judge[i].x, judge[i].y, store[j].x, store[j].y); 
		}
	}
	sort(e+1,e+1+cnt,cmp);
	for(int i=1;i<=cnt;i++)
	{
		int x = e[i].x, y = e[i].y;
		if(!visx[x] && !visy[y])
		{
			ans += e[i].w;
			visx[x] = 1; visy[y] = 1;
		}
	}
	cout << fixed << setprecision(6) << ans << endl;

	return 0;
}

E. E Just Passing Through

题意：

​ 从最西边（也就是第一列的某一点）出发，到达最东边某一点结束，过程中必须要经过p个通道（不多也不少），求权值和最小。

复盘：

​ 读完题目，是dp没毛病，而且还是三维的，f[i] [j] [k] 表示在第 i 行第 j 列 通过了 k 个通道的权值最小值，最后再遍历比较一下f [i] [m] [p] 即可。

坑点：

​ 当a[i] [j] = -1 时，是不能走的（这都知道），并且（敲重点）它旁边是没有通道的，题目讲了（肯定是字太小了，就没看到了）。就这个坑点，debug de 了一个多小时 …

Code：

#include <bits/stdc++.h>
#define int long long
using namespace std;

const int N=505;
int a[N][N], vis[N][N], ans=1e9;
int n,m,p,cnt;
int f[N][N][12];
int mn(int a,int b,int c)
{
    return min(a,min(b,c));
}

void init()
{
    for(int i=0;i<=n+1;i++)
    {
        for(int j=0;j<=m+1;j++)
        {
            for(int k=0;k<=p+1;k++)
            {
                f[i][j][k] = 1e9;
            }
        }
    }
    for(int i=0;i<=n;i++)
    {
        f[i][0][0] = 0;
    }
}

signed main()
{
    //freopen("070.in","r",stdin);
    //freopen("070.out","w",stdout);
    cin>>n>>m>>p;
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            cin>>a[i][j];
        }
    }
    for(int i=2;i<n;i++)
    {
        for(int j=2;j<m;j++)
        {
            if(a[i][j]!=-1 && a[i-1][j]!=-1 && a[i+1][j]!=-1 && a[i][j-1]!=-1 && a[i][j+1]!=-1 && a[i-1][j]>a[i][j] && a[i+1][j]>a[i][j] && a[i][j-1]<a[i][j] && a[i][j+1]<a[i][j])
            {
                vis[i][j]=2;
            }
        }
    }
    init();

    for(int j=1;j<=m;j++)
    {
        for(int i=1;i<=n;i++)
        {
            if(a[i][j]!=-1)
            {
                for(int k=0;k<=p;k++)
                {
                    if(vis[i][j]!=2){
                        //f[i][j][k] = 1e9;
                        f[i][j][k] = min(f[i][j][k], mn(f[i-1][j-1][k],f[i][j-1][k],f[i+1][j-1][k])+a[i][j]);
                        //cout<<i-1<<' '<<j-1<<' '<<k<<' '<<f[i-1][j-1][k]<<" zzz"<<endl;
                        //cout<<i<<' '<<j<<' '<<k<<' '<<f[i][j][k]<<endl;
                    }
                    else{
                        f[i][j][k+1] = min(f[i][j][k+1], mn(f[i-1][j-1][k],f[i][j-1][k],f[i+1][j-1][k])+a[i][j]);
                        //cout<<i<<' '<<j<<' '<<k+1<<' '<<f[i][j][k+1]<<endl;
                    }
                }
            }
        }
    }
    int ans=1e9;
    for(int i=1;i<=n;i++)
    {
        ans = min(ans, f[i][m][p]);
    }
    if(ans != 1e9) cout<<ans<<endl;
    else cout<<"impossible"<<endl;

    return 0;
}

G.G Out of Sorts

题意：

​ 给一组序列，问用二分查找每一个数，有多少数能查询成功 (n <= 1e6)

思路：

​ n 才1e6， 暴力枚举也才 nlogn

Code：

#include <bits/stdc++.h>
#define int long long
using namespace std;

const int N=1000005;
int x[N], ans;

signed main()
{
    int n,m,a,c;
    cin>>n>>m>>a>>c>>x[0];
    for(int i=1;i<=n;i++)
    {
        x[i] = (a*x[i-1] + c) % m;
    }
    for(int i=1;i<=n;i++)
    {
        int l=1, r=n, k=x[i];
        while(l <= r)
        {
            int md = (l+r)>>1;
            if(x[md] == k) { ans++; break; }
            if(k > x[md]) l = md+1;
            if(k < x[md]) r = md-1;
        }
    }
    cout<<ans<<endl;

    return 0;
}

​ 赛后又想了想，发现可以用 logn 的算法实现：去二分查所有的数，对于查到的数，能满足条件的，ans++。

​ 那么怎样才算满足条件？

​ 假设现在已经查到了x[mid] , 接下来要去查找区间 [l, mid-1] 和 [mid+1, r], 对于区间 [l, mid-1] 的查询，<x[mid] 才算有效的，对于区间 [mid+1, r] 的查询，>x[mid] 才算有效的，所以要更新上下限，具体操作看代码。

Code：

#include <bits/stdc++.h>
#define int long long
using namespace std;

const int N=1000005;
int x[N], ans;
void binary(int l,int r,int mn,int mx)
{
	if(l > r || mx < mn) return ;
	int md = (l+r) >> 1;
	if(x[md] > mn && x[md] < mx) ans++;
	binary(md+1, r, max(mn, x[md]), mx);
	binary(l, md-1, mn, min(mx, x[md]));
}

signed main()
{
	//freopen("1.out","w",stdout);
    int n,m,a,c;
    cin>>n>>m>>a>>c>>x[0];
    for(int i=1;i<=n;i++)
    {
        x[i] = (a*x[i-1] + c) % m;
    }
	int mx=1e12;
	binary(1,n,-1,mx);
    cout<<ans<<endl;

    return 0;
}