Codeforces Beta Round #2

最新推荐文章于 2022-09-19 15:21:11 发布

原创最新推荐文章于 2022-09-19 15:21:11 发布 · 644 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#up #ini #struct #c

本文解析了一场编程比赛中三道题目的解决方案，包括通过遍历寻找最大值的策略、使用动态规划解决路径选择问题以减少乘积后的0的数量，以及采用模拟退火算法来逼近最优解。

我对第二场的唯一感想是……B题写得真tmd丑- -|||

A题：遍历一次找最大值，再遍历一次找第一个超过m的ren。不会用map，表示鸭梨甚大……

#include <stdio.h> #include <string.h> typedef struct { char name[100]; int val; }win; win ord[1005]; win a[1005]; char ans[1005][100]; int main() { int i,j,n,up,cnt,s; scanf("%d",&n); up=0; for (i=0;i<n;i++) { scanf("%s%d",ord[i].name,&ord[i].val); } up=0; for (i=0;i<n;i++) { for (j=0;j<up;j++) { if (strcmp(ord[i].name,a[j].name)==0) break; } if (j==up) { strcpy(a[up].name,ord[i].name); a[up].val=ord[i].val; up++; } else { a[j].val+=ord[i].val; } } s=-100000000; cnt=0; for (i=0;i<up;i++) { if (a[i].val>s) { s=a[i].val; } } for (i=0;i<up;i++) { if (a[i].val==s) { strcpy(ans[cnt],a[i].name); cnt++; } } if (cnt==1) { printf("%s",ans[0]); return 0; } for (i=0;i<up;i++) { a[i].val=0; } for (i=0;i<n;i++) { for (j=0;j<up;j++) { if (strcmp(a[j].name,ord[i].name)==0) break; } a[j].val+=ord[i].val; if (a[j].val>=s) { for (j=0;j<cnt;j++) { if (strcmp(ans[j],ord[i].name)==0) break; } if (j<cnt) break; } } printf("%s/n",ans[j]); return 0; }

B题……先写了190+行，结果debug又加了四五十行，代码实在难看……

题目说找一条路径使得将路径上所有数字乘起来后缀的0最少，那么dp两次，第一次找2最少的，第二次找5最少的。注意路径上可能出现0，所以如果有0另外算，然后dp时将0跳过。

代码如下

#include <stdio.h> #include <string.h> #define INF 999999999 int dp2[1005][1005]; int dp5[1005][1005]; int map[1005][1005]; int two[1005][1005]; int five[1005][1005]; char ret[3005]; int Get2(int t) { int s=0; while(t%2==0) { t/=2; s++; } return s; } int Get5(int t) { int s=0; while(t%5==0) { t/=5; s++; } return s; } int Min(int x,int y) { return x<y?x:y; } int main() { int i,j,n,ans,x,y,t,p,q,up; scanf("%d",&n); t=INF; for (i=0;i<n;i++) { for (j=0;j<n;j++) { scanf("%d",&map[i][j]); if (map[i][j]==0) { t=1; x=i; y=j; } else { two[i][j]=Get2(map[i][j]); five[i][j]=Get5(map[i][j]); } } } memset(dp2,0,sizeof(dp2)); memset(dp5,0,sizeof(dp5)); for (i=0;i<n;i++) { for (j=0;j<n;j++) { if (map[i][j]==0) { dp2[i][j]=INF; dp5[i][j]=INF; } else if (i>0 && j>0) { if (dp2[i-1][j]>dp2[i][j-1] || dp2[i-1][j]==dp2[i][j-1] && dp5[i-1][j]>dp5[i][j-1]) { dp2[i][j]=dp2[i][j-1]+two[i][j]; dp5[i][j]=dp5[i][j-1]+five[i][j]; } else { dp2[i][j]=dp2[i-1][j]+two[i][j]; dp5[i][j]=dp5[i-1][j]+five[i][j]; } } else if (i>0) { dp2[i][j]=dp2[i-1][j]+two[i][j]; dp5[i][j]=dp5[i-1][j]+five[i][j]; } else if (j>0) { dp2[i][j]=dp2[i][j-1]+two[i][j]; dp5[i][j]=dp5[i][j-1]+five[i][j]; } else { dp2[i][j]=two[i][j]; dp5[i][j]=five[i][j]; } } } ans=Min(dp2[n-1][n-1],dp5[n-1][n-1]); if (ans<t) { p=n-1; q=n-1; up=0; while(p!=0 || q!=0) { // printf("%d %d/n",p,q); if (q==0) { p--; ret[up++]='D'; } else if (p==0) { q--; ret[up++]='R'; } else if (dp2[p-1][q]+two[p][q]==dp2[p][q] && dp5[p-1][q]+five[p][q]==dp5[p][q]) { p--; ret[up++]='D'; } else { q--; ret[up++]='R'; } } } memset(dp2,0,sizeof(dp2)); memset(dp5,0,sizeof(dp5)); for (i=0;i<n;i++) { for (j=0;j<n;j++) { if (map[i][j]==0) { dp2[i][j]=INF; dp5[i][j]=INF; } else if (i>0 && j>0) { if (dp5[i-1][j]>dp5[i][j-1] || dp5[i-1][j]==dp5[i][j-1] && dp2[i-1][j]>dp2[i][j-1]) { dp2[i][j]=dp2[i][j-1]+two[i][j]; dp5[i][j]=dp5[i][j-1]+five[i][j]; } else { dp2[i][j]=dp2[i-1][j]+two[i][j]; dp5[i][j]=dp5[i-1][j]+five[i][j]; } } else if (i>0) { dp2[i][j]=dp2[i-1][j]+two[i][j]; dp5[i][j]=dp5[i-1][j]+five[i][j]; } else if (j>0) { dp2[i][j]=dp2[i][j-1]+two[i][j]; dp5[i][j]=dp5[i][j-1]+five[i][j]; } else { dp2[i][j]=two[i][j]; dp5[i][j]=five[i][j]; } } } if (ans<=Min(dp2[n-1][n-1],dp5[n-1][n-1]) && ans<=t) { printf("%d/n",ans); for (i=up-1;i>=0;i--) { printf("%c",ret[i]); } printf("/n"); return 0; } if (t<ans && t<=Min(dp2[n-1][n-1],dp5[n-1][n-1])) { printf("%d/n",t); for (i=0;i<x;i++) { printf("D"); } for (i=0;i<y;i++) { printf("R"); } for (i=x+1;i<n;i++) { printf("D"); } for (i=y+1;i<n;i++) { printf("R"); } printf("/n"); return 0; } printf("%d/n",Min(dp2[n-1][n-1],dp5[n-1][n-1])); p=n-1; q=n-1; up=0; while(p!=0 || q!=0) { if (q==0) { p--; ret[up++]='D'; } else if (p==0) { q--; ret[up++]='R'; } else if (dp2[p-1][q]+two[p][q]==dp2[p][q] && dp5[p-1][q]+five[p][q]==dp5[p][q]) { p--; ret[up++]='D'; } else { q--; ret[up++]='R'; } } for (i=up-1;i>=0;i--) { printf("%c",ret[i]); } printf("/n"); return 0; }

C题

模拟退火……因为sin(A/2)=R/d，要算关于角A之间的差就要让A/2变成A，因此需要平方……原来用fabs显然过不了……

#include <stdio.h> #include <math.h> #define eps 1e-5 typedef struct { double x,y,r; }circle; int dx[4]={1,0,-1,0}; int dy[4]={0,1,0,-1}; circle a[3]; double sqr(double t) { return t*t; } double Dist(double x1,double y1,double x2,double y2) { return sqrt(sqr(x1-x2)+sqr(y1-y2)); } double f(double x,double y) { int i; double ret=0; double sub[3]; for (i=0;i<3;i++) { sub[i]=Dist(x,y,a[i].x,a[i].y)/a[i].r; } ret+=sqr(sub[1]-sub[0])+sqr(sub[2]-sub[1])+sqr(sub[2]-sub[0]); return ret; } int main() { int i,j,n,p; double min,step,xx,yy,tag; xx=0; yy=0; for (i=0;i<3;i++) { scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].r); xx+=a[i].x; yy+=a[i].y; } step=1; xx/=3; yy/=3; while(step>eps) { // printf("%lf....%lf/n",xx,yy); min=f(xx,yy); p=-1; for (i=0;i<4;i++) { tag=f(xx+dx[i]*step,yy+dy[i]*step); if (tag<min) { min=tag; p=i; } } if (p==-1) step/=2; else { xx=xx+dx[p]*step; yy=yy+dy[p]*step; } } if (f(xx,yy)<=eps) { printf("%.5lf %.5lf/n",xx,yy); } return 0; }