2019牛客暑期多校训练营（第七场) H Pair（数位DP+记忆化搜索）

最新推荐文章于 2020-11-12 16:23:52 发布

Qo0

最新推荐文章于 2020-11-12 16:23:52 发布

阅读量337

点赞数

分类专栏： 0x5D数位统计DP 2019牛客多校暑期训练营 0x22深度优先搜索

本文链接：https://blog.csdn.net/qq_41835683/article/details/98987176

版权

2019牛客多校暑期训练营同时被 3 个专栏收录

6 篇文章 0 订阅

订阅专栏

0x22深度优先搜索

6 篇文章 0 订阅

订阅专栏

0x5D数位统计DP

2 篇文章 0 订阅

订阅专栏

题目链接：https://ac.nowcoder.com/acm/contest/887/H

题意：给你一个三个整数 A,B,C ， $1\leq x\leq A$ ， $1\leq y\leq B$ ,问至少满足下面任意一个条件的 <x,y> 对数
$(x\ and\ y)>C$

$(x\ xor\ y)<C$

题解：数位DP+记忆化搜索，能写出来的都是神仙呀，Orz神仙们%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%，菜鸡只能赛后补题，代码来自牛逼网友，%%%%%

$int \ dfs(int\ pos,int \ \_and,int \ \_xor,int \ \_up1,int \ \_up2)$

简单解释一下各个参数的意义， pos 是下标， $\_and$ 是其中的一种情况： $&$ $\&$ 运算， $\_xor$ 是另一种情况 xor 运算， $\_up1$ 是枚举的第 pos 位的最大上限， $\_up2$ 是枚举的第 pos 的最大上限，这些说了之后，看代码能更好的理解一些。

简单说一下， A,B,C 都是在 1e9 范围内，那么30位顶死，把 A,B,C 转换成30位的字符串，从高位往低位扫，遇到满足情况则其实就是一种情况，枚举每一位的时候，该位的 x,y 最大上限是由上一位决定的，然后对每一位进行暴力枚举。

最后还要减去枚举x,y时，任意一方出现0的所有情况，0和小于C的正整数。

简单附代码：

# define _CRT_SECURE_NO_WARNINGS
#pragma GCC optimize(2)
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <vector>
#define inf 0x3f3f3f3f
#define Pair pair<int,int>
#define int long long
#define fir first
#define sec second
 
namespace fastIO {
#define BUF_SIZE 100000
#define OUT_SIZE 100000
  //fread->read
  bool IOerror = 0;
  //inline char nc(){char ch=getchar();if(ch==-1)IOerror=1;return ch;}
  inline char nc() {
    static char buf[BUF_SIZE], * p1 = buf + BUF_SIZE, * pend = buf + BUF_SIZE;
    if (p1 == pend) {
      p1 = buf; pend = buf + fread(buf, 1, BUF_SIZE, stdin);
      if (pend == p1) { IOerror = 1; return -1; }
    }
    return *p1++;
  }
  inline bool blank(char ch) { return ch == ' ' || ch == '\n' || ch == '\r' || ch == '\t'; }
  template<class T> inline bool read(T& x) {
    bool sign = 0; char ch = nc(); x = 0;
    for (; blank(ch); ch = nc());
    if (IOerror)return false;
    if (ch == '-')sign = 1, ch = nc();
    for (; ch >= '0' && ch <= '9'; ch = nc())x = x * 10 + ch - '0';
    if (sign)x = -x;
    return true;
  }
  inline bool read(double& x) {
    bool sign = 0; char ch = nc(); x = 0;
    for (; blank(ch); ch = nc());
    if (IOerror)return false;
    if (ch == '-')sign = 1, ch = nc();
    for (; ch >= '0' && ch <= '9'; ch = nc())x = x * 10 + ch - '0';
    if (ch == '.') { double tmp = 1; ch = nc(); for (; ch >= '0' && ch <= '9'; ch = nc())tmp /= 10.0, x += tmp * (ch - '0'); }
    if (sign)x = -x; return true;
  }
  inline bool read(char* s) {
    char ch = nc();
    for (; blank(ch); ch = nc());
    if (IOerror)return false;
    for (; !blank(ch) && !IOerror; ch = nc())* s++ = ch;
    *s = 0;
    return true;
  }
  inline bool read(char& c) {
    for (c = nc(); blank(c); c = nc());
    if (IOerror) { c = -1; return false; }
    return true;
  }
  template<class T, class... U>bool read(T& h, U& ... t) { return read(h) && read(t...); }
#undef OUT_SIZE
#undef BUF_SIZE
};
using namespace fastIO;
using namespace std;
 
 
/*---------------------------------------------------------------------------------------------------------------------------*/
const double eps = 1e-8;
struct Point { double x, y; Point(double x = 0, double y = 0) :x(x), y(y) {} };
typedef Point Vector;
Vector operator + (Vector A, Vector B) { return Vector(A.x + B.x, A.y + B.y); }
Vector operator - (Point A, Point B) { return Vector(A.x - B.x, A.y - B.y); }
Vector operator * (Vector A, double p) { return Vector(A.x * p, A.y * p); }
Vector operator / (Vector A, double p) { return Vector(A.x / p, A.y / p); }
bool operator < (const Point& A, const Point& B) { return A.x < B.x || (A.x == B.x && A.y < B.y); }
int dcmp(double x) { if (fabs(x) < eps)return 0; else return x < 0 ? -1 : 1; }
bool operator == (const Point& A, const Point& B) { return dcmp(A.x - B.x) == 0 && dcmp(A.y - B.y) == 0; }
double Dot(Vector A, Vector B) { return A.x * B.x + A.y * B.y; }
double Length(Vector A) { return sqrt(Dot(A, A)); }
double Dist2(const Point& A, const Point& B) { return (A.x - B.x) * (A.x - B.x) + (A.y - B.y) * (A.y - B.y); }
double Angle(Vector A, Vector B) { return acos(Dot(A, B) / Length(A) / Length(B)); }
double Cross(Vector A, Vector B) { return A.x * B.y - A.y * B.x; }
double Area2(Point A, Point B, Point C) { return Cross(B - A, C - A); }
double DistanceToLine(Point P, Point A, Point B) { Vector v1 = B - A, v2 = P - A; return fabs(Cross(v1, v2)) / Length(v1); }
Vector Normal(Vector A) { double L = Length(A); return Vector(-A.y / L, A.x / L); }
Vector Rotate(Vector A, double rad) { return Vector(A.x * cos(rad) - A.y * sin(rad), A.x * sin(rad) + A.y * cos(rad)); }
Point GetLineIntersection(Point P, Vector v, Point Q, Vector w) { Vector u = P - Q; double t = Cross(w, u) / Cross(v, w); return P + v * t; }
double PolygonArea(vector<Point>p) {//duo area
  double area = 0;
  int n = p.size();
  for (int i = 1; i < n - 1; i++)area += Cross(p[i] - p[0], p[i + 1] - p[0]);
  return area / 2;
}
double torad(double deg) { return deg / 180 * acos(-1); }
vector<Point> ConvexHull(vector<Point>& p) {
  sort(p.begin(), p.end());
  p.erase(unique(p.begin(), p.end()), p.end());
  int n = p.size();
  int m = 0;
  vector<Point> ch(n + 1);
  for (int i = 0; i < n; i++) { while (m > 1 && Cross(ch[m - 1] - ch[m - 2], p[i] - ch[m - 2]) <= 0)m--; ch[m++] = p[i]; }
  int k = m;
  for (int i = n - 2; i >= 0; i--) { while (m > k && Cross(ch[m - 1] - ch[m - 2], p[i] - ch[m - 2]) <= 0)m--; ch[m++] = p[i]; }
  if (n > 1)m--;
  ch.resize(m);
  return ch;
}
double diameter2(vector<Point>& points, Point& A, Point& B) {
  vector<Point>p = ConvexHull(points);
  int n = p.size();
  if (n == 1) return 0;
  if (n == 2) return Dist2(p[0], p[1]);
  //p.push_back(p[0]); // 免得取模
  double ans = 0;
  for (int u = 0, v = 1; u < n; u++) {
    for (;;) {
      //Cross(p[u+1]-p[u],p[v+1]-p[v])<=0时停止旋转
      int diff = Cross(p[u + 1] - p[u], p[v + 1] - p[v]);
      if (diff <= 0) {
        if (Dist2(p[u], p[v]) > ans) {
          ans = Dist2(p[u], p[v]); // u和v是对踵点
          A = p[u]; B = p[v];
        }
        //if(diff==0) ans=max(ans,Dist2(p[u],p[v+1])); // diff==0 时u和v+1也是对踵点
        break;
      }
      v = (v + 1) % n;
    }
  }
  return ans;
}
/*---------------------------------------------------------------------------------------------------------------------------*/
struct Line {
  Point P;
  Vector v;
  double ang;
  Line() {}
  Line(Point P, Vector v) :P(P), v(v) { ang = atan2(v.y, v.x); }
  bool operator < (const Line& L)const {
    return ang < L.ang;
  }
};
Point GetLineIntersection(const Line& a, const Line& b) {
  Vector u = a.P - b.P;
  double t = Cross(b.v, u) / Cross(a.v, b.v);
  return a.P + a.v * t;
}
bool OnLeft(Line L, Point p) { return Cross(L.v, p - L.P) > 0; }
vector<Point>HalfplaneIntersection(vector<Line>L) {
  int n = L.size();
  sort(L.begin(), L.end());
  int first, last;
  vector<Point> p(n);
  vector<Line> q(n);
  vector<Point> ans;
  q[first = last = 0] = L[0];
  for (int i = 1; i < n; i++) {
    while (first < last && !OnLeft(L[i], p[last - 1])) last--;
    while (first < last && !OnLeft(L[i], p[first]))first++;
    q[++last] = L[i];
    if (fabs(Cross(q[last].v, q[last - 1].v)) < eps) {
      last--;
      if (OnLeft(q[last], L[i].P))q[last] = L[i];
    }
    if (first < last)p[last - 1] = GetLineIntersection(q[last - 1], q[last]);
  }
  while (first < last && !OnLeft(q[first], p[last - 1]))last--;//删除无用平面
  if (last - first <= 1)return ans;
  p[last] = GetLineIntersection(q[last], q[first]);//计算首尾两个半平面的交点
  //从deque复制到输出中
  for (int i = first; i <= last; i++)ans.push_back(p[i]);
  return ans;
}
bool OnSegment(const Point& p, const Point& a1, const Point& a2) {
  return dcmp(Cross(a1 - p, a2 - p)) == 0 && dcmp(Dot(a1 - p, a2 - p)) < 0;
}
bool SegmentProperIntersection(const Point& a1, const Point& a2, const Point& b1, const Point& b2) {
  double c1 = Cross(a2 - a1, b1 - a1), c2 = Cross(a2 - a1, b2 - a1),
    c3 = Cross(b2 - b1, a1 - b1), c4 = Cross(b2 - b1, a2 - b1);
  return dcmp(c1) * dcmp(c2) < 0 && dcmp(c3) * dcmp(c4) < 0;
}
int IsPointInPolygon(const Point& p, const vector<Point>& poly) {
  int wn = 0;
  int n = poly.size();
  for (int i = 0; i < n; i++) {
    const Point& p1 = poly[i];
    const Point& p2 = poly[(i + 1) % n];
    if (p1 == p || p2 == p || OnSegment(p, p1, p2)) return -1; // 在边界上
    int k = dcmp(Cross(p2 - p1, p - p1));
    int d1 = dcmp(p1.y - p.y);
    int d2 = dcmp(p2.y - p.y);
    if (k > 0 && d1 <= 0 && d2 > 0) wn++;
    if (k < 0 && d2 <= 0 && d1 > 0) wn--;
  }
  if (wn != 0) return 1; // 内部
  return 0; // 外部
}
 
bool ConvexPolygonDisjoint(const vector<Point> ch1, const vector<Point> ch2) {
  int c1 = ch1.size();
  int c2 = ch2.size();
  cout << "c1 " << c1 << " " << "c2 " << c2 << endl;
  for (int i = 0; i < c1; i++)
    if (IsPointInPolygon(ch1[i], ch2) != 0) return false; // 内部或边界上
  for (int i = 0; i < c2; i++)
    if (IsPointInPolygon(ch2[i], ch1) != 0) return false; // 内部或边界上
  for (int i = 0; i < c1; i++)
    for (int j = 0; j < c2; j++)
      if (SegmentProperIntersection(ch1[i], ch1[(i + 1) % c1], ch2[j], ch2[(j + 1) % c2])) return false;
  return true;
}
/*---------------------------------------------------------------------------------------------------------------------------*/
 
const int N = 200 + 5;
const double pi = acos(-1.0);
typedef long long ll;
//const int mod = 998244353;
 
 
 
int f[31][3][3][2][2];
string sa, sb, sc;
int rp;
int dfs(int pos, int _and, int _xor, int _up1, int _up2) {
  //rp++;
  //printf("%lld %lld %lld %lld %lld\n",pos,_and,_xor,_up1,_up2);
  if (pos == 30) {
    if (_and == 1 || _xor == 2)return 1;
    else return 0;
  }
  if (f[pos][_and][_xor][_up1][_up2] != -1)return f[pos][_and][_xor][_up1][_up2];
  f[pos][_and][_xor][_up1][_up2] = 0;
  int up1 = _up1 ? (sa[pos] - '0') : 1;
  int up2 = _up2 ? (sb[pos] - '0') : 1;
  int n_and, n_xor;
  for (int i = 0; i <= up1; i++) {
    for (int j = 0; j <= up2; j++) {
      int x = i & j, y = i ^ j;
      n_and = _and;
      n_xor = _xor;
      if (_and == 0) {
        if (x > sc[pos] - '0')n_and = 1;
        else if (x < sc[pos] - '0')n_and = 2;
      }
      if (_xor == 0) {
        if (y > sc[pos] - '0')n_xor = 1;
        else if (y < sc[pos] - '0')n_xor = 2;
      }
      f[pos][_and][_xor][_up1][_up2] += dfs(pos + 1, n_and, n_xor, _up1 && (i == up1), _up2 && (j == up2));
    }
  }
  return f[pos][_and][_xor][_up1][_up2];
}
 
signed main() {
  int t; read(t);
  while (t--) {
    memset(f, -1, sizeof f);
    sa.clear(); sb.clear(); sc.clear();
    int a, b, c; read(a, b, c);
 
    /*printf("%lld\n",1<<30);
    int res=0;
    for(int i=1;i<=a;i++){
      for(int j=1;j<=b;j++){
        if((i&j)>c || (i^j)<c)res++;
      }
    }
    printf("%lld ",res);*/
 
    int ans = min(c, a+1) + min(c, b+1)-1;
    int i = 0;
    while (i < 30) { sa += (a & 1) + '0'; a >>= 1; i++; }
    reverse(sa.begin(), sa.end()); i = 0;
    while (i < 30) { sb += (b & 1) + '0'; b >>= 1; i++; }
    reverse(sb.begin(), sb.end()); i = 0;
    while (i < 30) { sc += (c & 1) + '0'; c >>= 1; i++; }
    reverse(sc.begin(), sc.end());
    /*printf("%s\n%s\n%s\n",sa.c_str(),sb.c_str(),sc.c_str());*/
    /*printf("%lld\n",min(c-1,a));
    printf("%lld\n",min(c-1,b));*/
    printf("%lld\n", dfs(0, 0, 0, 1, 1) - ans);
 
  }
  return 0;
}
 
/*
  1
  1 6 5
 
  2 2 2
*/

Qo0

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
2019牛客暑期多校训练营（第七场) H Pair（数位DP+记忆化搜索）

题目链接：https://ac.nowcoder.com/acm/contest/887/H题意：给你一个三个整数，，,问至少满足下面任意一个条件的对数题解：数位DP+记忆化搜索，能写出来的都是神仙呀，Orz神仙们%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%，菜鸡只能赛后补题，代码来自牛逼网友，%%%%%简单解释一下各个参数的意义，是下标，是其中的一种情况：运算...
复制链接

扫一扫

专栏目录