SDOI2018 Round2 简要题解

最新推荐文章于 2019-05-01 18:06:55 发布
最新推荐文章于 2019-05-01 18:06:55 发布
阅读量750 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/wxh010910/article/details/80382638
版权

Day 1

物理实验

转一下坐标系，扫描线，用set可以维护出每一段的最低位置能看到的线段的斜率，然后再双指针扫一遍就行了。

#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
#define Debug(...) fprintf(stderr, __VA_ARGS__)

typedef long long LL;
typedef long double LD;
typedef unsigned int uint;
typedef pair <int, int> pii;
typedef unsigned long long uLL;

template <typename T> inline void Read(T &x) {
  char c = getchar();
  bool f = false;
  for (x = 0; !isdigit(c); c = getchar()) {
    if (c == '-') {
      f = true;
    }
  }
  for (; isdigit(c); c = getchar()) {
    x = x * 10 + c - '0';
  }
  if (f) {
    x = -x;
  }
}

template <typename T> inline bool CheckMax(T &a, const T &b) {
  return a < b ? a = b, true : false;
}

template <typename T> inline bool CheckMin(T &a, const T &b) {
  return a > b ? a = b, true : false;
}

const int N = 10005;
const LD eps = 1e-9;

struct Point {
  LD x, y;

  Point(LD x = 0, LD y = 0):x(x), y(y) {}

  Point operator - (const Point &b) const {
    return Point(x - b.x, y - b.y);
  }
} c, d, a[N], b[N];

struct Line {
  Point l, r;
  int i;

  Line(Point l = Point(0, 0), Point r = Point(0, 0), int i = -1):l(l), r(r), i(i) {}

  inline LD F(LD x) const {
    return ((x - l.x) * r.y + (r.x - x) * l.y) / (r.x - l.x);
  }

  inline LD S() const {
    return sqrt(1 + ((r.y - l.y) * (r.y - l.y)) / ((r.x - l.x) * (r.x - l.x)));
  }

  bool operator < (const Line &b) const {
    LD x = max(l.x, b.l.x);
    return fabs(F(x)) < fabs(b.F(x));
  }
} l[N];

LD all, ans, v[N << 1], x[N << 1];
pair <LD, int> eve[N << 1];
set <Line> s, t;
int n, m;

int Main() {
  Read(n), m = 0, ans = 0, s.clear(), t.clear();
  for (int i = 1; i <= n; ++i) {
    Read(a[i].x), Read(a[i].y), Read(b[i].x), Read(b[i].y);
  }
  Read(c.x), Read(c.y), Read(d.x), Read(d.y), Read(all), d = d - c;
  LD p, q, len = sqrt(d.x * d.x + d.y * d.y);
  for (int i = 1; i <= n; ++i) {
    a[i] = a[i] - c;
    p = (d.x * a[i].x + d.y * a[i].y) / len, q = (d.x * a[i].y - d.y * a[i].x) / len;
    a[i] = Point(p, q);
    b[i] = b[i] - c;
    p = (d.x * b[i].x + d.y * b[i].y) / len, q = (d.x * b[i].y - d.y * b[i].x) / len;
    b[i] = Point(p, q);
    if (a[i].x > b[i].x) {
      swap(a[i], b[i]);
    }
    l[i] = Line(a[i], b[i], i);
    eve[i - 1 << 1] = mp(a[i].x, i), eve[i - 1 << 1 | 1] = mp(b[i].x, -i);
  }
  sort(eve, eve + (n << 1));
  for (int i = 0; i < n << 1; ++i) {
    if (!i || abs(eve[i].X - x[m - 1]) > eps) {
      x[m] = eve[i].X;
      if (m) {
        v[m - 1] = 0;
        if (!s.empty()) {
          v[m - 1] += s.begin()->S();
        }
        if (!t.empty()) {
          v[m - 1] += t.begin()->S();
        }
      }
      m++;
    }
    if (eve[i].Y > 0) {
      int p = eve[i].Y;
      if (a[p].y > 0) {
        s.insert(l[p]);
      } else {
        t.insert(l[p]);
      }
    } else {
      int p = -eve[i].Y;
      if (a[p].y > 0) {
        s.erase(l[p]);
      } else {
        t.erase(l[p]);
      }
    }
  }
  LD cur;
  cur = 0;
  for (int i = 0, j = 0; i < m; ++i) {
    for (; j + 1 < m && x[j + 1] < x[i] + all; ++j) {
      cur += (x[j + 1] - x[j]) * v[j];
    }
    CheckMax(ans, cur + (x[i] + all - x[j]) * (j == m - 1 ? 0 : v[j]));
    cur -= (x[i + 1] - x[i]) * v[i];
  }
  cur = 0;
  for (int i = m - 1, j = m - 1; ~i; --i) {
    for (; j && x[j - 1] > x[i] - all; --j) {
      cur += (x[j] - x[j - 1]) * v[j - 1];
    }
    CheckMax(ans, cur + (x[j] + all - x[i]) * (!j ? 0 : v[j - 1]));
    if (i) {
      cur -= (x[i] - x[i - 1]) * v[i - 1];
    }
  }
  printf("%.9Lf\n", ans);
  return 0;
}

int main() {
#ifdef wxh010910
  freopen("d.in", "r", stdin);
#endif
  int T;
  for (Read(T); T; --T) {
    Main();
  }
#ifdef wxh010910
  Debug("My Time: %.3lfms\n", (double)clock() / CLOCKS_PER_SEC);
#endif
  return 0;
}
战略游戏

建出圆方树，询问就是问虚树上圆点个数。

#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
#define Debug(...) fprintf(stderr, __VA_ARGS__)

typedef long long LL;
typedef long double LD;
typedef unsigned int uint;
typedef pair <int, int> pii;
typedef unsigned long long uLL;

template <typename T> inline void Read(T &x) {
  char c = getchar();
  bool f = false;
  for (x = 0; !isdigit(c); c = getchar()) {
    if (c == '-') {
      f = true;
    }
  }
  for (; isdigit(c); c = getchar()) {
    x = x * 10 + c - '0';
  }
  if (f) {
    x = -x;
  }
}

template <typename T> inline bool CheckMax(T &a, const T &b) {
  return a < b ? a = b, true : false;
}

template <typename T> inline bool CheckMin(T &a, const T &b) {
  return a > b ? a = b, true : false;
}

const int N = 200005;
const int M = 20;

int n, m, q, ans, cnt, tim, top, dep[N], dfn[N], low[N], seq[N], siz[N], sta[N], sum[N], f[M][N];
vector <int> adj[N], adv[N];

inline void DFS(int x, int p) {
  dfn[x] = low[x] = ++tim, sta[++top] = x;
  bool f = false;
  for (auto y : adj[x]) {
    if (y != p || f) {
      if (!dfn[y]) {
        DFS(y, x), CheckMin(low[x], low[y]);
        if (dfn[x] <= low[y]) {
          ++m, adv[x].pb(m);
          for (int k = 0; k != y; adv[m].pb(k = sta[top--]));
        }
      } else {
        CheckMin(low[x], dfn[y]);
      }
    } else {
      f = true;
    }
  }
}

inline void DFS(int x) {
  dfn[x] = ++tim, siz[x] = 1, sum[x] = sum[f[0][x]] + (x <= n);
  for (int i = 1; i < M; ++i) {
    f[i][x] = f[i - 1][f[i - 1][x]];
  }
  for (auto y : adv[x]) {
    dep[y] = dep[x] + 1, f[0][y] = x, DFS(y), siz[x] += siz[y];
  }
}

inline int LCA(int x, int y) {
  if (dep[x] < dep[y]) {
    swap(x, y);
  }
  if (dep[x] > dep[y]) {
    for (int i = 0; i < M; ++i) {
      if (dep[x] - dep[y] >> i & 1) {
        x = f[i][x];
      }
    }
  }
  if (x == y) {
    return x;
  }
  for (int i = M - 1; ~i; --i) {
    if (f[i][x] != f[i][y]) {
      x = f[i][x], y = f[i][y];
    }
  }
  return f[0][x];
}

int Main() {
  Read(n), Read(m), tim = top = 0;
  for (int i = 1, x, y; i <= m; ++i) {
    Read(x), Read(y);
    adj[x].pb(y), adj[y].pb(x);
  }
  m = n, DFS(1, 0), tim = 0, DFS(1);
  for (Read(q); q; --q) {
    Read(cnt), ans = -cnt;
    for (int i = 1; i <= cnt; ++i) {
      Read(seq[i]);
    }
    sort(seq + 1, seq + cnt + 1, [&](const int &x, const int &y) {
      return dfn[x] < dfn[y];
    });
    for (int i = cnt; i > 1; --i) {
      seq[++cnt] = LCA(seq[i], seq[i - 1]);
    }
    sort(seq + 1, seq + cnt + 1, [&](const int &x, const int &y) {
      return dfn[x] < dfn[y];
    });
    cnt = unique(seq + 1, seq + cnt + 1) - seq - 1, top = 0;
    for (int i = 1; i <= cnt; ++i) {
      for (; top && !(dfn[sta[top]] <= dfn[seq[i]] && dfn[sta[top]] + siz[sta[top]] >= dfn[seq[i]] + siz[seq[i]]); --top);
      if (top) {
        ans += sum[f[0][seq[i]]] - sum[sta[top]];
      }
      sta[++top] = seq[i];
    }
    for (int i = 1; i <= cnt; ++i) {
      if (seq[i] <= n) {
        ++ans;
      }
    }
    printf("%d\n", ans);
  }
  for (int i = 1; i <= m; ++i) {
    adj[i].clear(), adv[i].clear(), dfn[i] = 0;
  }
  return 0;
}

int main() {
#ifdef wxh010910
  freopen("d.in", "r", stdin);
#endif
  int T;
  for (Read(T); T; --T) {
    Main();
  }
#ifdef wxh010910
  Debug("My Time: %.3lfms\n", (double)clock() / CLOCKS_PER_SEC);
#endif
  return 0;
}
反回文串

枚举串的周期 d d ，计算最短周期等于自身长度 d 的串的个数记为 f f ，不难发现 f(n)=k⌈n2⌉−∑d|nf(d) ，答案是 ∑d|nf(d)(d1+[dmod2=0]) ∑ d | n f ( d ) ( d 1 + [ d mod 2 = 0 ] ) 。计算 f f 的时候可以发现它是一个类似积性函数的东西，推一推就可以算了。

#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
#define Debug(...) fprintf(stderr, __VA_ARGS__)

typedef long long LL;
typedef long double LD;
typedef unsigned int uint;
typedef pair <int, int> pii;
typedef unsigned long long uLL;

template <typename T> inline void Read(T &x) {
  char c = getchar();
  bool f = false;
  for (x = 0; !isdigit(c); c = getchar()) {
    if (c == '-') {
      f = true;
    }
  }
  for (; isdigit(c); c = getchar()) {
    x = x * 10 + c - '0';
  }
  if (f) {
    x = -x;
  }
}

template <typename T> inline bool CheckMax(T &a, const T &b) {
  return a < b ? a = b, true : false;
}

template <typename T> inline bool CheckMin(T &a, const T &b) {
  return a > b ? a = b, true : false;
}

const int N = 200000;
const int M = 20;

LL n, k, a[N], pri[N];
int m, cnt, mod, b[N];

inline int Rand() {
#ifdef wxh010910
  return rand() << 15 | rand();
#else
  return rand();
#endif
}

inline LL RandLL() {
  return (LL)Rand() << 30 | Rand();
}

inline LL Qul(LL x, LL y, LL mod) {
  return ((x * y - (LL)(((LD)x * y + 0.5) / mod) * mod) % mod + mod) % mod;
}

inline LL Qow(LL x, LL y, LL mod) {
  LL r = 1;
  for (; y; y >>= 1, x = Qul(x, x, mod)) {
    if (y & 1) {
      r = Qul(r, x, mod);
    }
  }
  return r;
}

inline bool MillerRabin(LL n) {
  if (n == 2) {
    return true;
  }
  if (n < 2 || !(n & 1)) {
    return false;
  }
  LL a, x, y, u = n - 1;
  int t = 0;
  for (; !(u & 1); u >>= 1, ++t);
  for (int i = 0; i < 10; ++i) {
    a = RandLL() % (n - 1) + 1, x = Qow(a, u, n);
    for (int j = 0; j < t; ++j, x = y) {
      if ((y = Qul(x, x, n)) == 1 && x != 1 && x != n - 1) {
        return false;
      }
    }
    if (x != 1) {
      return false;
    }
  }
  return true;
}

inline LL PollardRho(LL a, LL c) {
  LL k = 2, x = RandLL() % a, y = x, p = 1;
  for (LL i = 1; p == 1; ++i) {
    x = (Qul(x, x, a) + c) % a, p = __gcd(abs(x - y), a);
    if (i == k) {
      y = x, k <<= 1;
    }
  }
  return p;
}

inline void Divide(LL x) {
  if (x == 1) {
    return ;
  }
  if (MillerRabin(x)) {
    pri[cnt++] = x;
    return ;
  }
  LL p = x;
  while (p >= x) {
    p = PollardRho(x, RandLL() % (x - 1));
  }
  Divide(p), Divide(x / p);
}

inline int Qow(int x, LL y) {
  int r = 1;
  for (; y; y >>= 1, x = 1LL * x * x % mod) {
    if (y & 1) {
      r = 1LL * r * x % mod;
    }
  }
  return r;
}

inline int DFS(int x, LL cur, int val) {
  if (x == m) {
    if (cur & 1) {
      return 1LL * cur % mod * Qow(k, cur + 1 >> 1) % mod * val % mod;
    } else {
      return 1LL * (cur >> 1) % mod * Qow(k, cur >> 1) % mod * val % mod;
    }
  }
  LL ans = 0, cof = 1LL * val * ((1 - a[x]) % mod + mod) % mod;
  if (a[x] & 1) {
    ans += DFS(x + 1, cur, cof);
  }
  for (int i = 1; i < b[x]; ++i) {
    ans += DFS(x + 1, cur *= a[x], cof);
  }
  ans += DFS(x + 1, cur *= a[x], val);
  return ans % mod;
}

int Main() {
  Read(n), Read(k), Read(mod), k %= mod, cnt = 0, Divide(n);
  sort(pri, pri + cnt), m = 0;
  for (int i = 0; i < cnt; ++i) {
    if (!i || pri[i] != pri[i - 1]) {
      a[m++] = pri[i], b[m - 1] = 0;
    }
    ++b[m - 1];
  }
  printf("%d\n", DFS(0, 1, 1));
  return 0;
}

int main() {
#ifdef wxh010910
  freopen("d.in", "r", stdin);
#endif
  int T;
  for (Read(T); T; --T) {
    Main();
  }
#ifdef wxh010910
  Debug("My Time: %.3lfms\n", (double)clock() / CLOCKS_PER_SEC);
#endif
  return 0;
}

Day 2

原题识别

求出欧拉序，对于 (x,y) ，我们在 (dfnx,dfny) ( d f n x , d f n y ) 和 (endx,endy) ( e n d x , e n d y ) 两处权值加上 ansx,y a n s x , y ，在 (dfnx,endy) ( d f n x , e n d y ) 和 (dfny,endx) ( d f n y , e n d x ) 处权值减去 ansx,y a n s x , y ，那么操作二就是询问一个子矩形的权值和。

考虑对每种颜色分别算贡献，将产生贡献的矩形暴力搞出来，然后离散化，然后暴力给全局打上标记即可。查询可以扫描线加树状数组解决。可以证明产生贡献的矩形个数期望是线性的。

#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
#define Debug(...) fprintf(stderr, __VA_ARGS__)

typedef long long LL;
typedef long double LD;
typedef unsigned int uint;
typedef pair <int, int> pii;
typedef unsigned long long uLL;

template <typename T> inline void Read(T &x) {
  char c = getchar();
  bool f = false;
  for (x = 0; !isdigit(c); c = getchar()) {
    if (c == '-') {
      f = true;
    }
  }
  for (; isdigit(c); c = getchar()) {
    x = x * 10 + c - '0';
  }
  if (f) {
    x = -x;
  }
}

template <typename T> inline bool CheckMax(T &a, const T &b) {
  return a < b ? a = b, true : false;
}

template <typename T> inline bool CheckMin(T &a, const T &b) {
  return a > b ? a = b, true : false;
}

const int N = 200005;
const int M = 2005;

struct Event {
  int x, l, r, f;

  Event(int x = 0, int l = 0, int r = 0, int f = 0):x(x), l(l), r(r), f(f) {}

  bool operator < (const Event &b) const {
    return x < b.x;
  }
};

int n, m, t, cnt, a[N], s[N], p[N], q[N], dic[N], per[N], tag[M][M];
LL ans[N], bit[N], bitx[N], bity[N], bitxy[N];
vector <Event> eve, qry;
vector <int> adj[N];
uint SA, SB, SC;

inline uint Rand() {
  SA ^= SA << 16, SA ^= SA >> 5, SA ^= SA << 1;
  uint t = SA;
  SA = SB, SB = SC, SC ^= t ^ SA;
  return SC;
}

inline void DFS(int x) {
  p[x] = ++t, s[t] = 1;
  for (auto y : adj[x]) {
    DFS(y);
  }
  q[x] = ++t, s[t] = -1;
}

inline void AddDic(int x) {
  dic[++cnt] = p[x] - 1, dic[++cnt] = q[x];
  for (auto y : adj[x]) {
    dic[++cnt] = p[y] - 1, dic[++cnt] = q[y];
  }
}

inline void AddTag(int xl, int xr, int yl, int yr, int f) {
  xl = lower_bound(dic, dic + cnt + 1, xl - 1) - dic + 1;
  yl = lower_bound(dic, dic + cnt + 1, yl - 1) - dic + 1;
  xr = lower_bound(dic, dic + cnt + 1, xr) - dic + 1;
  yr = lower_bound(dic, dic + cnt + 1, yr) - dic + 1;
  tag[xl][yl] += f, tag[xr][yr] += f;
  tag[xr][yl] -= f, tag[xl][yr] -= f;
}

inline void AddRec(int x) {
  AddTag(1, n << 1, p[x], q[x], 1);
  AddTag(p[x], q[x], 1, n << 1, 1);
  for (auto y : adj[x]) {
    AddTag(p[y], q[y], p[y], q[y], -2);
  }
}

inline void Modify(int x, int y, int a) {
  LL b = 1LL * a * s[x - 1], c = 1LL * a * s[y - 1], d = 1LL * a * s[x - 1] * s[y - 1];
  for (int i = y; i <= n << 1; i += i & -i) {
    bit[i] += a, bitx[i] += b, bity[i] += c, bitxy[i] += d;
  }
}

inline LL Query(int x, int y) {
  LL a = 0, b = 0, c = 0, d = 0;
  for (int i = y; i; i -= i & -i) {
    a += bit[i], b += bitx[i], c += bity[i], d += bitxy[i];
  }
  return a * s[x] * s[y] - b * s[y] - c * s[x] + d;
}

int Main() {
  Read(n), Read(m), Read(SA), Read(SB), Read(SC);
  eve.clear(), qry.clear(), t = 0;
  for (int i = 1; i <= n; ++i) {
    adj[i].clear(), per[i] = i;
  }
  for (int i = 2; i <= m; ++i) {
    adj[i - 1].pb(i);
  }
  for (int i = m + 1; i <= n; ++i) {
    adj[Rand() % (i - 1) + 1].pb(i);
  }
  for (int i = 1; i <= n; ++i) {
    a[i] = Rand() % n + 1;
  }
  DFS(1);
  for (int i = 1; i <= n << 1; ++i) {
    bit[i] = bitx[i] = bity[i] = bitxy[i] = 0;
    s[i] += s[i - 1];
  }
  sort(per + 1, per + n + 1, [&](const int &x, const int &y) {
    return a[x] < a[y];
  });
  for (int i = 1, j = 1; i <= n; i = j) {
    for (dic[0] = 0, dic[cnt = 1] = n << 1; j <= n && a[per[j]] == a[per[i]]; ++j);
    for (int k = i; k < j; ++k) {
      AddDic(per[k]);
    }
    sort(dic, dic + cnt + 1), cnt = unique(dic, dic + cnt + 1) - dic - 1;
    for (int k = 1; k <= cnt; ++k) {
      for (int l = 1; l <= cnt; ++l) {
        tag[k][l] = 0;
      }
    }
    for (int k = i; k < j; ++k) {
      AddRec(per[k]);
    }
    for (int k = 1; k <= cnt; ++k) {
      for (int l = 1; l <= cnt; ++l) {
        tag[k][l] += tag[k - 1][l] + tag[k][l - 1] - tag[k - 1][l - 1];
        if (tag[k][l]) {
          eve.pb(Event(dic[k - 1] + 1, dic[l - 1] + 1, dic[l], 1));
          eve.pb(Event(dic[k] + 1, dic[l - 1] + 1, dic[l], -1));
        }
      }
    }
  }
  Read(m);
  for (int i = 1, o, x, y, xl, xr, yl, yr; i <= m; ++i) {
    Read(o), Read(x), Read(y), ans[i] = 0;
    if (o == 1) {
      xl = xr = p[x], yl = yr = p[y];
    } else {
      xl = 1, xr = p[x], yl = 1, yr = p[y];
    }
    qry.pb(Event(xl - 1, yl, yr, -i));
    qry.pb(Event(xr, yl, yr, i));
  }
  sort(eve.begin(), eve.end());
  sort(qry.begin(), qry.end());
  for (int i = 0, j = 0; i < eve.size() || j < qry.size(); ) {
    if (i < eve.size() && (j == qry.size() || eve[i].x <= qry[j].x)) {
      Modify(eve[i].x, eve[i].l, eve[i].f);
      Modify(eve[i].x, eve[i].r + 1, -eve[i].f);
      ++i;
    } else {
      if (qry[j].f > 0) {
        ans[qry[j].f] += Query(qry[j].x, qry[j].r) - Query(qry[j].x, qry[j].l - 1);
      } else {
        ans[-qry[j].f] -= Query(qry[j].x, qry[j].r) - Query(qry[j].x, qry[j].l - 1);
      }
      ++j;
    }
  }
  for (int i = 1; i <= m; ++i) {
    printf("%lld\n", ans[i]);
  }
  return 0;
}

int main() {
#ifdef wxh010910
  freopen("d.in", "r", stdin);
#endif
  int T;
  for (Read(T); T; --T) {
    Main();
  }
#ifdef wxh010910
  Debug("My Time: %.3lfms\n", (double)clock() / CLOCKS_PER_SEC);
#endif
  return 0;
}
旧试题

根据rng_58恒等式，答案为 ∑i∑j∑k[gcd(i,j)=1][gcd(i,k)=1][gcd(j,k)=1]⌊Ai⌋⌊Bj⌋⌊Ck⌋ ∑ i ∑ j ∑ k [ gcd ( i , j ) = 1 ] [ gcd ( i , k ) = 1 ] [ gcd ( j , k ) = 1 ] ⌊ A i ⌋ ⌊ B j ⌋ ⌊ C k ⌋ 。

然后贴一下 蒜头的奖杯 代码就行了。

#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
#define Debug(...) fprintf(stderr, __VA_ARGS__)

typedef long long LL;
typedef long double LD;
typedef unsigned int uint;
typedef pair <int, int> pii;
typedef unsigned long long uLL;

template <typename T> inline void Read(T &x) {
  char c = getchar();
  bool f = false;
  for (x = 0; !isdigit(c); c = getchar()) {
    if (c == '-') {
      f = true;
    }
  }
  for (; isdigit(c); c = getchar()) {
    x = x * 10 + c - '0';
  }
  if (f) {
    x = -x;
  }
}

template <typename T> inline bool CheckMax(T &a, const T &b) {
  return a < b ? a = b, true : false;
}

template <typename T> inline bool CheckMin(T &a, const T &b) {
  return a > b ? a = b, true : false;
}

const int N = 100005;
const int mod = 1e9 + 7;

int n, A, B, C, ans, tot, a[N], b[N], c[N], d[N], e[N], f[N], p[N], q[N], r[N], s[N], t[N], u[N], v[N], w[N], pri[N];
bool chk[N];

inline void F(int *a, int n) {
  for (int i = 1; i <= n; ++i) {
    for (int j = i + i; j <= n; j += i) {
      a[j] = (a[j] - a[i] + mod) % mod;
    }
  }
}

inline void G(int *a, int n) {
  for (int i = 1; pri[i] <= n; ++i) {
    for (int j = n / pri[i]; j; --j) {
      a[j] = (a[j] + a[j * pri[i]]) % mod;
    }
  }
}

inline void H(int *a, int n) {
  for (int i = 1; pri[i] <= n; ++i) {
    for (int j = 1; j * pri[i] <= n; ++j) {
      a[j * pri[i]] = (a[j * pri[i]] + a[j]) % mod;
    }
  }
}

int Main() {
  Read(A), Read(B), Read(C), n = max(max(A, B), C), ans = tot = 0;
  for (int i = 2; i <= n; ++i) {
    if (!chk[i]) {
      pri[++tot] = i;
    }
    for (int j = 1; i * pri[j] <= n; ++j) {
      chk[i * pri[j]] = true;
      if (i % pri[j] == 0) {
        break;
      }
    }
  }
  pri[tot + 1] = n + 1;
  for (int i = 1; i <= n; ++i) {
    a[i] = A / i, b[i] = B / i, c[i] = C / i, d[i] = e[i] = f[i] = i == 1;
  }
  G(c, n), F(e, n), F(f, n);
  for (int i = 1; i <= n; ++i) {
    int m = n / i;
    for (int j = 1; j <= m; ++j) {
      p[j] = e[i * j], q[j] = f[i * j], r[j] = c[i * j], s[j] = a[i * j], t[j] = b[i * j], w[j] = d[i * j];
    }
    F(w, m);
    for (int x = 1; x * x <= m; ++x) {
      for (int j = 1; j <= m; ++j) {
        u[j] = v[j] = 0;
      }
      for (int j = x; j <= m; j += x) {
        u[j] = s[j], v[j] = t[j];
      }
      G(u, m), G(v, m);
      for (int j = 1; j <= m; ++j) {
        u[j] = 1LL * u[j] * w[j] % mod, v[j] = 1LL * v[j] * w[j] % mod;
      }
      H(u, m), H(v, m);
      for (int j = 1; j <= m; ++j) {
        u[j] = 1LL * u[j] * t[j] % mod, v[j] = 1LL * v[j] * s[j] % mod;
      }
      for (int y = x; x * y <= m; ++y) {
        if (__gcd(x, y) == 1) {
          LL g = 0, h = 0;
          for (int j = y; j <= m; j += y) {
            g += u[j], h += v[j];
          }
          g %= mod, h %= mod;
          ans = (1LL * p[x] * q[y] % mod * r[x * y] % mod * g + ans) % mod;
          if (x != y) {
            ans = (1LL * p[y] * q[x] % mod * r[x * y] % mod * h + ans) % mod;
          }
        }
      }
    }
  }
  printf("%d\n", ans);
  return 0;
}

int main() {
#ifdef wxh010910
  freopen("d.in", "r", stdin);
#endif
  int T;
  for (Read(T); T; --T) {
    Main();
  }
#ifdef wxh010910
  Debug("My Time: %.3lfms\n", (double)clock() / CLOCKS_PER_SEC);
#endif
  return 0;
}
荣誉称号

不难发现只有前 2k+1 2 k + 1 个值有用，然后暴力树形DP就行了。

#include <bits/stdc++.h>

using namespace std;

#define X first
#define Y second
#define mp make_pair
#define pb push_back
#define Debug(...) fprintf(stderr, __VA_ARGS__)

typedef long long LL;
typedef long double LD;
typedef unsigned int uint;
typedef pair <int, int> pii;
typedef unsigned long long uLL;

template <typename T> inline void Read(T &x) {
  char c = getchar();
  bool f = false;
  for (x = 0; !isdigit(c); c = getchar()) {
    if (c == '-') {
      f = true;
    }
  }
  for (; isdigit(c); c = getchar()) {
    x = x * 10 + c - '0';
  }
  if (f) {
    x = -x;
  }
}

template <typename T> inline bool CheckMax(T &a, const T &b) {
  return a < b ? a = b, true : false;
}

template <typename T> inline bool CheckMin(T &a, const T &b) {
  return a > b ? a = b, true : false;
}

const int N = 2050;
const int M = 205;
const LL inf = 1LL << 60;

int n, m, k, p, A, B, cur, lim, dis[N];
LL sum[N], f[N][M], g[N][M];
uint SA, SB, SC;

inline uint Rand() {
  SA ^= SA << 16, SA ^= SA >> 5, SA ^= SA << 1;
  uint t = SA;
  SA = SB, SB = SC, SC ^= t ^ SA;
  return SC;
}

inline void Modify(int i, int x, int y) {
  for (; i >> cur > lim; cur += k);
  i >>= cur, x %= m;
  sum[i] += y, g[i][0] += 1LL * (m - x) * y, g[i][x] -= m * y;
}

inline void DFS(int x) {
  if (x << 1 > lim) {
    for (int i = 0; i < m; ++i) {
      f[x][i] = g[x][i];
    }
    return ;
  }
  for (int i = 0; i < m; ++i) {
    f[x][i] = inf;
  }
  DFS(x << 1);
  if ((x << 1 | 1) > lim || dis[x << 1] != dis[x << 1 | 1]) {
    for (int i = 0; i < m; ++i) {
      for (int j = 0; j < m; ++j) {
        CheckMin(f[x][(i + j) % m], g[x][i] + f[x << 1][j]);
      }
    }
    return ;
  }
  DFS(x << 1 | 1);
  for (int i = 0; i < m; ++i) {
    for (int j = 0; j < m; ++j) {
      CheckMin(f[x][(i + j) % m], g[x][i] + f[x << 1][j] + f[x << 1 | 1][j]);
    }
  }
}

int Main() {
  Read(n), Read(k), Read(m), Read(p), Read(SA), Read(SB), Read(SC), Read(A), Read(B), ++k, cur = 0, lim = min((1 << k) - 1, n);
  for (int i = 1; i <= lim; ++i) {
    sum[i] = 0;
    for (int j = 0; j < m; ++j) {
      g[i][j] = 0;
    }
  }
  for (int i = 1, x, y; i <= n; ++i) {
    if (i <= p) {
      Read(x), Read(y);
    } else {
      x = Rand() % A + 1, y = Rand() % B + 1;
    }
    Modify(i, x, y);
  }
  for (int i = 1; i <= lim; ++i) {
    for (int j = 1; j < m; ++j) {
      g[i][j] += g[i][j - 1];
    }
    for (int j = 1; j < m; ++j) {
      g[i][j] += sum[i] * j;
    }
  }
  for (int i = lim; i; --i) {
    dis[i] = i << 1 <= lim ? dis[i << 1] + 1 : 1;
  }
  DFS(1);
  printf("%lld\n", f[1][0]);
  return 0;
}

int main() {
#ifdef wxh010910
  freopen("d.in", "r", stdin);
#endif
  int T;
  for (Read(T); T; --T) {
    Main();
  }
#ifdef wxh010910
  Debug("My Time: %.3lfms\n", (double)clock() / CLOCKS_PER_SEC);
#endif
  return 0;
}
wxh010910
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
P4619 [SDOI2018]旧试题（莫比乌斯反演，建图优化三重枚举，三元环计数，神仙好题，超级清晰易懂）
繁凡さん的博客
06-20 471
整理的算法模板合集： ACM模板 点我看算法全家桶系列！！！ 实际上是一个全新的精炼模板整合计划 P4619 [SDOI2018]旧试题（莫比乌斯反演，三元环计数） Problem 计算： ∑i=1A∑j=1B∑k=1Cd(ijk)mod  (109+7)\sum_{i=1}^{A}\sum_{j=1}^{B}\sum_{k=1}^{C}d(ijk) \mod (10^9+7)i=1∑A​j=1∑B​k=1∑C​d(ijk)mod(109+7) 其中 d(ijk)d(ijk)d(ijk) 表示 i×j
CQOI2018简要题解
04-04
### CQOI2018简要题解 #### D1T1：破解D-H协议 **题意：** 题目要求解决的是一个基于Diffie-Hellman (D-H) 密钥交换协议的问题。该协议允许两个通信方（Alice 和 Bob）在不安全的通道上建立共享密钥。具体步骤如下...
SDOI2018：原题识别
weixin_33905756的博客
05-20 247
题解： https://files.cnblogs.com/files/clrs97/old-solution.pdf   Code： #include&lt;cstdio&gt; #include&lt;cstdlib&gt; #include&lt;algorithm&gt; using namespace std; typedef long long ll; const...
[SDOI2018]旧试题
weixin_34072458的博客
03-08 177
[SDOI2018]旧试题 Bzoj5332: [Sdoi2018]旧试题 反演+三元环计数+卡常 反演部分： 第一步还是“约数个数和”那个结论的扩展：（考虑质数的次数的分配即可证明） 最后一步就是枚举u,v,w考虑满足能凑出u，v，w的x,y,z有哪些 （式子u,v,w,lcm那里有打错了） 还是一脸不可做 三个点形成的三元组(u,v,w)贡献一定。 三元环...
BZOJ 5332 [Sdoi2018]旧试题
dilei7762的博客
01-08 118
题目链接 https://lydsy.com/JudgeOnline/problem.php?id=5332 题解 反演得到 ∑d=1min⁡(A,B)μ(d)∑e=1min⁡(A,C)μ(e)∑f=1min⁡(B,C)μ(f)F(lcm(d,e),A)F(lcm(d,f),B)F(lcm(e,f),C) \sum_{d=1}^{\min(...
loj #2564. 「SDOI2018」原题识别
OZY的博客
01-02 502
链接：https://loj.ac/problem/2564 肝了大半天这题。。 总算是肝出来了 第一问首先可以用莫队莽过去。。 那么问题就是第二问了 第二问的话，先考虑链怎么做？ 显然，链的话就退化为一个序列了 记录一下每个点，上一个点颜色和他相同的位置就可以记录答案了 距离来说，询问的是(x,y) 那么一个点now，设他上一个颜色一样的位置是z（上一个指的是祖先，如果没有就是0） 如果now是...
NOI2018day2后面两个题的题解pdf
07-22
《NOI2018day2后面两个题的题解pdf》是一份针对全国青少年信息学奥林匹克竞赛（NOI）2018年第二天比赛的深度解析文档，由出题人亲自讲解，旨在帮助参赛者理解和掌握比赛中的关键知识点。这份资源包含两个部分：...
WC2018讲课论文 题解
08-02
圆方树simple_immortalCO 一道有趣的字符串题 二维计算几何相关算法与实战应用_徐明宽 Problem-Setting On Topcoder OI 中的超现实数和不平等博弈问题 杜瑜皓 NP-Hard问题求解方法杂谈 ...即时战略 通道 周区划分 题解
BJOI2018 双人猜数游戏 题解
09-14
BJOI2018双人猜数游戏...总之，BJOI2018双人猜数游戏题解涉及到的关键知识点包括C++编程、二分搜索算法、动态规划以及策略优化。理解和掌握这些内容，不仅可以解决这个问题，也能在其他类似的算法问题中发挥重要作用。
重庆八中周赛Round14周赛题解
最新发布
10-03
【重庆八中周赛Round14周赛题解】 本次周赛中包含了三道题目，分别是T1修改序列，T2 Knuth表示法，以及T3魔力塔。以下是每道题目的详细解答： **T1 修改序列 (modify)** 这道题目要求修改一个序列，使得序列的和...
SDOI2018 原题识别（主席树）
WAautomaton的博客
02-14 761
题目链接 题目大意 给定nnn个节点的树，其中包含一条非随机生成的长度为kkk的链，剩下的节点均随机父节点连边。每个节点有一个随机的颜色，维护： 1.给定x,yx,yx,y，求x,yx,yx,y之间不同颜色数。 2.给定x,yx,yx,y，对于所有满足分别在x,yx,yx,y到根的路径上的点a,ba,ba,b，求其询问1的答案之和。 n≤105,m≤2×105n\le 10^5,m\le 2\ti...
SDOI 2018 R2 爆零滚粗退役记
HT008的Blog
05-07 948
先把坑占下…
SDOI 2018 R2 自我总结
was__n的博客
05-19 326
DAY1        由于蒟蒻一轮省选爆炸CE，二轮上了考场之后并不是很慌，先看了看第一道题，感觉很可做，然后几乎整场比赛就栽在第一道题上了.... 第一眼的思路就是现将整个坐标系按照给定直线导轨的角度翻转，然后将挡板（线段）映射到导轨上，然后处理一下导轨之间的遮挡情况，注意到最后选取的区间左右端点一定有一个会卡在挡板的端点上，所以可以O(n)统计答案。       感觉很难码，然后就赶紧开始码...
SDOI2018 二轮培训划水记
weixin_30367169的博客
06-29 90
\(\mathcal{Day -1}\) 嗯，虽然说\(rqy\)看我们这么懵O，并不建议我们去掺和这种神仙打架，但是为了逃文化课学习OI并参观膜拜各路神仙，我毅然决然地选择去参加这次培训…… 这次一定要选一个靠前排的座位……并且要和\(rqy\)靠的近些不然就真的啥都不会了……这次培训居然不让带手机……真是不知其可\(qwq\) 唉……希望比上次好点吧，上次培训的题目我还没整呢！ \(emmmm...
SDOI 2018 R2 T1 物理实验题解
was__n的博客
05-19 366
题目链接：点击打开链接loj      首先将所有挡板翻转至以直线导轨为X轴的坐标系中，记录每个挡板两个端点的坐标，然后我们考虑如何处理挡板之间的遮挡情况，我们可以使用扫描线+set维护，我们将每个端点按x值排序，左端点记1 ，右端点记-1,然后我们可以用set维护在每段区间最下方的直线斜率是多少，遇到左端点在set中插入线段，遇到右端点，删除该线段。      最后我们考虑如何O(n）的统计答案...
bzoj5332&loj2565 [Sdoi2018]旧试题
Elijahqi
05-26 518
http://www.elijahqi.win/archives/3492 题目地址https://loj.ac/problem/2565 https://www.lydsy.com/JudgeOnline/problem.php?id=5332 首先知道结论 原题所求可以转化为 ∑i=1A∑j=1B∑k=1C∑x|i∑y|j∑z|k[gcd(x,y)=1][gcd(x,z)=1][gc...
SDOI2018R2续命记
薇小薇
05-17 1234
序（tan90°） 一步 Day -x（x&amp;amp;amp;amp;amp;amp;gt;=1) 因为知道自己很菜，国级的比赛根本不会摘金摘银，在大家都去CTSC/APIO的时候，窝还是在机房孑然一身，好吧还有zyz R2前差不多一个周就回家蹲着了，被盯着学习= =，看看以前做过的题窝怎么什么都不会啊 Day 0 因为不跟大部队一起，怎么去也是一个问题，最后窝的父母决定带着窝去，据说是如果退役了能有安慰QwQ，于是乎也不跟l...
SDOI2018日常
Coco_T的博客
04-08 1141
Day-1 终于等到这一天，明天出发 省选，OI成就梦想的第一步 恍恍惚惚，就到了今天 曾经也想过放弃，但是终究坚持下来了，心里还是有点小开心的 几个月前，就开始制定计划，一切都是这么的潜移默化 一个月之前，听到的是：没事放平心态，不行就回去学文化课 但是感觉随着时间的推移，大家对我的。。。（不知道是信心还是什么）就越来越强 听到的也从原来的 “ 放平心态 ” 变成了 “ 等你的好...
[SDOI2018]旧试题 题解
GKxx
05-01 609
传送门 简单反演+神仙优化题。 首先我们知道 σ0(xy)=∑i∣x∑j∣y[(i,j)=1]\sigma_0(xy)=\sum\limits_{i|x}\sum\limits_{j|y}[(i,j)=1]σ0​(xy)=i∣x∑​j∣y∑​[(i,j)=1] σ0(xyz)=∑i∣x∑j∣y∑k∣z[(i,j)=1][(i,k)=1][(j,k)=1]\sigma_0(xyz)=\sum\lim...
AtCoder Beginner Contest 134 简要题解
05-23
AtCoder Beginner Contest 134 是一场 AtCoder 的入门级比赛，以下是每道题的简要题解： A - Dodecagon 题目描述：已知一个正十二边形的边长，求它的面积。 解题思路：正十二边形的内角为 $150^\circ$，因此可以...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值