ICPC数学题——补题(更新中)

Maths

2020 沈阳 I. Rise of Shadows

• Rise of Shadows

• 思路：

  1. 首先 $v_H=\frac{2\pi }{HM},v_M=\frac{2\pi }{M}$，每一次时针和分针重合的间隔时间计算为
     $$\begin{gathered} \Delta T(v_M-v_H)=2\pi \\ \Delta T=\frac{HM}{H-1} \end{gathered}$$
     每次在两针重合前后需要偏转到 $\alpha =\frac{2\pi A}{HM}$ 角度，此时需要时间为
     $$\begin{gathered} |t^{\prime }|(v_M-v_H)=\alpha \\ |t^{\prime }|=\frac{A}{H-1} \end{gathered}$$
     于是每个整数时刻(“integal moments”)即为所有满足 $\left(\frac{kHM-A}{H-1},\frac{kHM+A}{H-1}\right),k\in [0,H-1]$ 内的范围在 $[0,HM)$ (这个左闭右开的问题可以将两端情况合并到一个重合点)内的整数。

  2. 需要将问题转化为求 $-A\le t(H-1)(\mathrm{m}\mathrm{o}\mathrm{d}HM)\le A，t\in [0,HM)$，由于要求整数，便可以把范围边界约到“整数”(可取整)。又因为若 $a,b,c$ 为任意$3$个整数，$m$ 为正整数，且 $(m,c)=1$，则当$ac\equiv bc(\mathrm{m}\mathrm{o}\mathrm{d}m)$时，有$a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)$。所以为了满足互质，两边同除 $g=\gcd (H-1,HM)$，得到
     $$\frac{-A}{g}\le t\ast \frac{H-1}{g}(\mathrm{m}\mathrm{o}\mathrm{d}\frac{HM}{g})\le \frac{A}{g},t\in [0,\frac{HM}{g}]$$
     此时正负半轴都只有 $1,2,...,\frac{A}{g}$ 个解，算上重合时的情况，共有 $2\ast \frac{A}{g}+1$ 个解。还原到 $[0,HM)$ 上便有 $g\ast (2\ast \frac{A}{g}+1)$ 个解。

  3. 注意当 $A=\frac{HM}{2}$ 时，两个相邻重合情况前后重叠，于是 $[0,HM)$ 内的所有整数时刻都满足条件，需要特判。

  #include<bits/stdc++.h>
  using namespace std;
  #define rep(i, a, b) for(int i = (a); i <= (b); i++)
  #define per(i, a, b) for(int i = (a); i >= (b); i--)
  #define ll long long
  #define db double
  #define VI vector<int>
  #define PII pair<int, int>
  #define Pi 3.141592653589793
  const int INF = 0x7fffffff;
  const int N = 1e5 + 5;
  const db eps = 1e-10;
  ll H, M, A;
  db h, m, a;
  int main(){
      cin >> H >> M >> A;
      h = H, m = M, a = A;
      if(A * 2 == H * M) cout << H * M << endl;
      else{
          int x = __gcd(H * M, H - 1);
          cout << x * (2 * (A / x) + 1) << endl;
      }
  }
  /*
  5 5 4
  9
  3 5 1
  3
  */
  

2020 上海 I. Sky Garden

• Sky Garden

• 思路：由于需要走最短距离，当圆上两点角度相差大于 $\theta$ 时，经过原点会更优。计算得到
  $$\begin{gathered} \theta r<2r \\ \therefore \theta <2\approx 144° \end{gathered}$$
  于是可以固定一个轴线方向的点，最后 $\ast 2m$ 即可。我们统计一圈中角度相差在 $\theta$ 内的点一边有 $A$ 个，两边就有 $2A$ 个。分成四种情况讨论：

  1. 当角度相差 $\theta <2$ 时，两点之间路程为 $abs(R-r)+\min (R,r)\theta$，其中弧上路程 $\min (R,r)\theta$ 需要走 $2\ast (1+2+...+A)$ 个(因为有左右两边)，半径差路程 $abs(R-r)$ 需要走 $2A$ 个。但是，当 $r=R$，即在同一个圆上时， 两点被重复算了两遍；而在不同圆上，由于我们两次循环半径都不同，枚举的点没有重复。于是在同一圆上需要 $/2$ 。
  2. 当角度相差 $\theta >2$ 时，两点需要经过原点共有 $2m-2A-1$ 个点。同理当 $r=R$ 时需要 $/2$ 。
  3. 当角度相差 $\theta =0$ 时，不同半径之间只需要加上半径差

  最后每圈有 $2m$ 个点，而且当 $m\ne 1$ 时，由于原点存在，还需要加上所有圆上点到原点的距离即可。

  #include<bits/stdc++.h>
  #define rep(i, a, b) for(int i = a; i <= b; i++)
  #define db double
  #define ll long long
  #define Pi 3.141592653589793
  using namespace std;
  int n, m, A;
  db ans = 0, per;
  int main() {
      cin >> n >> m;
      A = floor(2.0 * m / Pi);
      per = Pi / (db)m;
      rep(R, 1, n){
          rep(r, R, n){
              db tmp = 0;
              //theta < 2
              tmp += min(r, R) * per * (A * (1 + A) / 2.0) * 2;
              //theta > 2
              tmp += (R + r) * (2 * m -2 * A - 1);
              if(r == R) tmp /= 2.0;  //只有在同圆内才有重复计算
              
              //theta == 0
              tmp += abs(R - r);
              //theta < 2
              tmp += abs(R - r) * 2 * A;
              ans += tmp;
          }
      }
      ans *= (2 * m);
      if(m != 1) ans += (n * (n + 1) / 2.0) * (2 * m); //原点
      printf("%.9lf", ans);
  }
  /*
  2 3
  175.4159265359
  2 4
  302.2477796076938
  2 5
  468.4955592153876
  */