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;
}