hdu 5618 Jam's problem again

  1. #include <bits/stdc++.h>
  2. #define MAXN 100010
  3. using namespace std;
  4. int n,k,T,xx;
  5. int ans[MAXN],c[MAXN],f[MAXN];
  6. struct Node{
  7.     int x,y,z,id;
  8. }a[100010],b[100010];
  9. inline int read(){
  10.     char ch;
  11.     bool f=false;
  12.     int res=0;
  13.     while (((ch=getchar())<'0'||ch>'9')&&ch!='-');
  14.     if (ch=='-')
  15.         f=true;
  16.     else
  17.         res=ch-'0';
  18.     while ((ch=getchar())>='0'&&ch<='9')
  19.         res=(res<<3)+(res<<1)+ch-'0';
  20.     return f?~res+1:res;
  21. }
  22. inline bool cmp(Node x,Node y){
  23.     if (x.x==y.x){
  24.         if (x.y==y.y){
  25.             if (x.z==y.z){
  26.                 return x.id<y.id;
  27.             }
  28.             else return x.z<y.z;
  29.         }
  30.         else return x.y<y.y;
  31.     }
  32.     else return x.x<y.x;
  33. }
  34. inline bool cmp1(Node x,Node y){
  35.     if (x.y==y.y)
  36.         return x.z<y.z;
  37.     else
  38.         return x.y<y.y;
  39. }
  40. inline int lowbit(int x){
  41.     return x&(-x);
  42. }
  43. inline void add(int x,int y){
  44.     while (x<=MAXN){
  45.         c[x]+=y,x+=lowbit(x);
  46.     }
  47. }
  48. inline int sum(int x){
  49.     int summ=0;
  50.     while (x>0){
  51.         summ+=c[x],x-=lowbit(x);
  52.     }
  53.     return summ;
  54. }
  55. void solve(int l,int r){
  56.     if (l>=r)
  57.         return;
  58.     int mid=(l+r)/2;
  59.     for (register int i=l;i<=r;i++){
  60.         b[i-l+1].x=0,b[i-l+1].y=a[i].y,b[i-l+1].z=a[i].z;
  61.         if (i>mid)
  62.             b[i-l+1].id=a[i].id;
  63.         else
  64.             b[i-l+1].id=-1;
  65.     }
  66.     sort(b+1,b+r-l+2,cmp1);
  67.     for (register int i=1;i<=r-l+1;++i)
  68.         if (b[i].id==-1)
  69.             add(b[i].z,1);
  70.         else
  71.             ans[b[i].id]+=sum(b[i].z);
  72.     for(register int i=1;i<=r-l+1;++i)
  73.         if(b[i].id==-1)
  74.             add(b[i].z,-1);
  75.     solve(l,mid);
  76.     solve(mid+1,r);
  77. }
  78. int main(){
  79.     //scanf("%d",&T);
  80.     T=read();
  81.     while (T--){
  82.         memset(c,0,sizeof(c));
  83.         memset(ans,0,sizeof(ans));
  84.         //scanf("%d",&n);
  85.         n=read();
  86.         for (register int i=1;i<=n;++i){
  87.             a[i].x=read(),a[i].y=read(),a[i].z=read(),a[i].id=i;
  88.             //scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].z);
  89.             //a[i].id=i;
  90.         }
  91.         sort(a+1,a+1+n,cmp);
  92.         for (register int i=n-1;i>=1;i--){
  93.             if (a[i].x==a[i+1].x&&a[i].y==a[i+1].y&&a[i].z==a[i+1].z)
  94.                 xx++;
  95.             else xx=0;
  96.             ans[a[i].id]+=xx;
  97.         }
  98.         solve(1,n);
  99.         for (register int i=1;i<=n;++i)
  100.             printf("%d\n",ans[i]);
  101.         /*for (register int i=1;i<=n;++i)
  102.             f[ans[i]]++;
  103.         for (register int i=0;i<n;++i)
  104.             printf("%d\n",f[i]);*/
  105.     }
  106.     return 0;
  107. }

　　

动态逆序对

  1. #include <bits/stdc++.h>
  2. #define MAXN 100010
  3. using namespace std;
  4. long long n,m,sum1,x,gg;
  5. long long ans[MAXN],f[MAXN],c[MAXN];
  6. struct Node{
  7.     long long  x,y,z,id;
  8. }a[MAXN],b[MAXN];
  9. inline long long read(){
  10.     char ch;
  11.     bool f=false;
  12.     long long res=0;
  13.     while (((ch=getchar())<'0'||ch>'9')&&ch!='-');
  14.     if (ch=='-')
  15.         f=true;
  16.     else
  17.         res=ch-'0';
  18.     while ((ch=getchar())>='0'&&ch<='9')
  19.         res=(res<<3)+(res<<1)+ch-'0';
  20.     return f?~res+1:res;
  21. }
  22. inline bool cmp(Node x,Node y){
  23.     if (x.x==y.x)
  24.         return x.y<y.y;
  25.     else
  26.         return x.x<y.x;
  27. }
  28. inline bool cmp1(Node x,Node y){
  29.     return x.z<y.z;
  30. }
  31. inline long long lowbit(long long x){
  32.     return x&(-x);
  33. }
  34. inline void add(long long x,long long y){
  35.     while (x<=n){
  36.         c[x]+=y,x+=lowbit(x);
  37.     }
  38. }
  39. inline long long sum(long long x){
  40.     long long summ=0;
  41.     while (x>0){
  42.         summ+=c[x],x-=lowbit(x);
  43.     }
  44.     return summ;
  45. }
  46. inline void sc(long long l,long long r){
  47.     if (l>=r)
  48.         return;
  49.     long long mid=(l+r)/2;
  50.     for (register long long i=l;i<=r;++i){
  51.         b[i-l+1]=a[i];
  52.         if (i>mid) b[i-l+1].id=1;
  53.     }
  54.     sort(b+1,b+r-l+2,cmp);
  55.     for (register long long i=1;i<=r-l+1;++i){
  56.         if (b[i].id==0)
  57.             add(b[i].y,1);
  58.         else
  59.             ans[b[i].z]+=sum(n)-sum(b[i].y);
  60.     }
  61.     for (register long long i=1;i<=r-l+1;++i)
  62.         if (b[i].id==0)
  63.             add(b[i].y,-1);
  64.     for (register long long i=r-l+1;i>=1;--i){
  65.         if (b[i].id==0)
  66.             add(b[i].y,1);
  67.         else
  68.             ans[b[i].z]+=sum(b[i].y);
  69.     }
  70.     for (register long long i=r-l+1;i>=1;--i){
  71.         if (b[i].id==0)
  72.             add(b[i].y,-1);
  73.     }
  74.     sc(l,mid);
  75.     sc(mid+1,r);
  76. }
  77. int main(){
  78.     n=read(),m=read();gg=n;
  79.     for (register long long i=1;i<=n;++i){
  80.         a[i].x=i,a[i].y=read();
  81.         f[a[i].y]=i;
  82.     }
  83.     for (register long long i=1;i<=m;++i)
  84.         a[f[read()]].z=gg--;
  85.     for (register long long i=1;i<=n;i++)
  86.         if (a[i].z==0)
  87.             a[i].z=gg--;
  88.     sort(a+1,a+n+1,cmp1);
  89.     sc(1,n);
  90.     for (register long long i=1;i<=n;i++)
  91.         sum1+=ans[i];
  92.     for (register long long i=n;i>n-m;--i){
  93.         printf("%lld\n",sum1);
  94.         sum1-=ans[i];
  95.     }
  96.     return 0;
  97. }

hdu 4742 Pinball Game

Description

　　给出n个三维点对(x,y,z)，求三维非严格最长上升子序列长度和最长上升子序列数量

Solution

　　这是一个三维LIS问题，可以用CDQ解决。

　　先按x排序，降一维，然后剩下y、z，在y上进行CDQ分治，按y的大小用前面的更新后面的。

　　z方向离散化后用树状数组维护。

  1. #include <bits/stdc++.h>
  2. #define MAXN 100010
  3. using namespace std;
  4. const int Mo=1<<30;
  5. int ans,T,n;
  6. int dp[MAXN],fa[MAXN],q1[MAXN],q2[MAXN];
  7. int g[MAXN<<1],p[MAXN];
  8. long long num[MAXN],q3[MAXN],ansn;
  9. int t1;
  10. inline int read() {
  11.     char ch;
  12.     bool f=false;
  13.     int res=0;
  14.     while (((ch=getchar())<'0'||ch>'9')&&ch!='-');
  15.     if (ch=='-')
  16.         f=true;
  17.     else
  18.         res=ch-'0';
  19.     while ((ch=getchar())>='0'&&ch<='9')
  20.         res=(res<<3)+(res<<1)+ch-'0';
  21.     return f?~res+1:res;
  22. }
  23. struct Node {
  24.     int x,y,z;
  25. }q[MAXN];
  26. inline bool cmp(Node x,Node y) {
  27.     if (x.x!=y.x)
  28.         return x.x<y.x;
  29.     if (x.y!=y.y)
  30.         return x.y<y.y;
  31.     return x.z<y.z;
  32. }
  33. inline bool cmp1(int x,int y) {
  34.     if (q[x].y!=q[y].y)
  35.         return q[x].y<q[y].y;
  36.     return x<y;
  37. }
  38. inline int lowbit(int x) {
  39.     return x&(-x);
  40. }
  41. inline void add(int x,int v,long long y,int maxx) {
  42.     while (x<=maxx) {
  43.         if (v>q2[x])
  44.             q2[x]=v,q3[x]=y;
  45.         else
  46.             if (v==q2[x])
  47.                 q3[x]+=y;
  48.         x+=lowbit(x);
  49.     }
  50. }
  51. inline int sum(int x,long long &nn) {
  52.     int summ=0;
  53.     while (x>0) {
  54.         if (summ<q2[x])
  55.             summ=q2[x],nn=q3[x];
  56.         else
  57.             if (q2[x]==summ)
  58.                 nn+=q3[x];
  59.         x-=lowbit(x);
  60.     }
  61.     return summ;
  62. }
  63. void sc(int t,int w) {
  64.     /*if (t>w)
  65.         return;*/
  66.     if (t==w) {
  67.         if (dp[t]>ans)
  68.             ans=dp[t],ansn=num[t];
  69.         else
  70.             if (dp[t]==ans) {
  71.                 ansn=ansn+num[t];
  72.             }
  73.         q1[t]=q[t].z;
  74.         return;
  75.     }
  76.     int mid=(t+w)>>1,le=t,ri=mid+1,tt1=w-t+1,maxx=mid-t+1,l,h;
  77.     memmove(g+t1,fa+t,tt1*sizeof(int));
  78.     for (int i=0;i<tt1;++i)
  79.         if (g[t1+i]<=mid)
  80.             fa[le++]=g[t1+i];
  81.         else
  82.             fa[ri++]=g[t1+i];
  83.     t1+=tt1;
  84.     sc(t,mid);
  85.     l=t1;
  86.     t1-=tt1;
  87.     for (int i=1;i<=maxx;++i)
  88.         q2[i]=0;
  89.     for (int i=t1;i<l;++i) {
  90.         if  (g[i]<=mid) {
  91.             h=lower_bound(q1+t,q1+t+maxx,q[g[i]].z)-(q1+t)+1;
  92.             add(h,dp[g[i]],num[g[i]],maxx);
  93.             continue;
  94.         }
  95.         h=lower_bound(q1+t,q1+t+maxx,q[g[i]].z)-(q1+t);
  96.         if (h>=maxx||q1[t+h]>q[g[i]].z)
  97.             h--;
  98.         if (h>=0) {
  99.             long long ttt=0;
  100.             int tt=sum(h+1,ttt)+1;
  101.             if (tt>dp[g[i]])
  102.                 dp[g[i]]=tt,num[g[i]]=ttt;
  103.             else
  104.                 if (tt==dp[g[i]]&&tt!=1)
  105.                     num[g[i]]+=ttt;
  106.         }
  107.     }
  108.     sc(mid+1,w);
  109.     merge(q1+t,q1+t+maxx,q1+t+maxx,q1+t+tt1,p);
  110.     memmove(q1+t,p,tt1*sizeof(int));
  111. }
  112. int main() {
  113.     T=read();
  114.     while (T--) {
  115.         n=read();
  116.         ansn=t1=ans=0;
  117.         for (int i=1;i<=n;++i)
  118.             q[i].x=read(),q[i].y=read(),q[i].z=read();
  119.         sort(q+1,q+n+1,cmp);
  120.         for (int i=1;i<=n;++i)
  121.             fa[i]=i;
  122.         for (int i=1;i<=n;++i){
  123.             dp[i]=1;
  124.             num[i]=1;
  125.         }
  126.         sort(fa+1,fa+n+1,cmp1);
  127.         sc(1,n);
  128.         printf("%d %lld\n",ans,ansn%Mo);
  129.     }
  130.     return 0;
  131. }

hdu 4456 crowd

Description

　　详见

Solution

　　本来想继续写cdq的（毕竟cdq专题2333）但网上一查题解发现好像写二维树状数组的人更多。

　　一想反正我也不会，我这么菜，那就学这个吧。（而且看起来码量小多了）

　　然后...我调了一上午QWQ

　　这题就是将矩形旋转45度，每个询问就相当于询问一个矩阵，用容斥口胡下然后用二维树状数组维护

  1. #ifndef _GLIBCXX_NO_ASSERT
  2. #include <cassert>
  3. #endif
  4. #include <cctype>
  5. #include <cerrno>
  6. #include <cfloat>
  7. #include <ciso646>
  8. #include <climits>
  9. #include <clocale>
  10. #include <cmath>
  11. #include <csetjmp>
  12. #include <csignal>
  13. #include <cstdarg>
  14. #include <cstddef>
  15. #include <cstdio>
  16. #include <cstdlib>
  17. #include <cstring>
  18. #include <ctime>
  19.  
  20. #if __cplusplus >= 201103L
  21. #include <ccomplex>
  22. #include <cfenv>
  23. #include <cinttypes>
  24. #include <cstdalign>
  25. #include <cstdbool>
  26. #include <cstdint>
  27. #include <ctgmath>
  28. #include <cwchar>
  29. #include <cwctype>
  30. #endif
  31.  
  32. // C++
  33. #include <algorithm>
  34. #include <bitset>
  35. #include <complex>
  36. #include <deque>
  37. #include <exception>
  38. #include <fstream>
  39. #include <functional>
  40. #include <iomanip>
  41. #include <ios>
  42. #include <iosfwd>
  43. #include <iostream>
  44. #include <istream>
  45. #include <iterator>
  46. #include <limits>
  47. #include <list>
  48. #include <locale>
  49. #include <map>
  50. #include <memory>
  51. #include <new>
  52. #include <numeric>
  53. #include <ostream>
  54. #include <queue>
  55. #include <set>
  56. #include <sstream>
  57. #include <stack>
  58. #include <stdexcept>
  59. #include <streambuf>
  60. #include <string>
  61. #include <typeinfo>
  62. #include <utility>
  63. #include <valarray>
  64. #include <vector>
  65.  
  66. #if __cplusplus >= 201103L
  67. #include <array>
  68. #include <atomic>
  69. #include <chrono>
  70. #include <condition_variable>
  71. #include <forward_list>
  72. #include <future>
  73. #include <initializer_list>
  74. #include <mutex>
  75. #include <random>
  76. #include <ratio>
  77. #include <regex>
  78. #include <scoped_allocator>
  79. #include <system_error>
  80. #include <thread>
  81. #include <tuple>
  82. #include <typeindex>
  83. #include <type_traits>
  84. #include <unordered_map>
  85. #include <unordered_set>
  86. #endif
  87. #define MAXN 4000010
  88. using namespace std;
  89. int t[MAXN],q[MAXN],p[MAXN],X[MAXN],Y[MAXN],Z[MAXN];
  90. int k,n,m,ma;
  91. inline int read() {
  92.     char ch;
  93.     bool f=false;
  94.     int res=0;
  95.     while (((ch=getchar())<'0'||ch>'9')&&ch!='-');
  96.     if (ch=='-')
  97.         f=true;
  98.     else
  99.         res=ch-'0';
  100.     while ((ch=getchar())>='0'&&ch<='9')
  101.         res=(res<<3)+(res<<1)+ch-'0';
  102.     return f?~res+1:res;
  103. }
  104. inline int find (int x) {
  105.     return lower_bound(q+1,q+k,x)-q;
  106. }
  107. inline int lowbit(int x) {
  108.     return x&(-x);
  109. }
  110. inline void Hash(int x,int y) {
  111.     int yy=y;
  112.     while (x<=ma) {
  113.         y=yy;
  114.         while (y<=ma) {
  115.             q[k++]=x*ma+y;
  116.             y+=lowbit(y);
  117.         }
  118.         x+=lowbit(x);
  119.     }
  120. }
  121. inline void add(int x,int y,int z) {
  122.     int yy=y;
  123.     while (x<=ma) {
  124.         y=yy;
  125.         while (y<=ma) {
  126.             //int zz=lower_bound(q+1,q+k,(x*ma+y))-q;
  127.             //t[zz]+=z;
  128.             t[find(x*ma+y)]+=z;
  129.             //printf("%d\n",y);
  130.             y+=lowbit(y);
  131.         }
  132.         x+=lowbit(x);
  133.     }
  134. }
  135. inline int sum(int x,int y) {
  136.     int summ=0,yy=y;
  137.     while (x) {
  138.         y=yy;
  139.         while (y) {
  140.             int z=find(x*ma+y);
  141.             if (q[z]==x*ma+y)
  142.                 summ+=t[z];
  143.             //printf("%d\n",z);
  144.             y-=lowbit(y);
  145.         }
  146.         x-=lowbit(x);
  147.     }
  148.     //printf("%d\n",summ);
  149.     return summ;
  150. }
  151. int main() {
  152.     while (true) {
  153.         n=read();
  154.         if (n==0)
  155.             return 0;
  156.         k=1,ma=2*n,m=read();
  157.         memset(t,0,sizeof(t));
  158.         for (int i=1;i<=m;++i){
  159.             p[i]=read(),X[i]=read(),Y[i]=read(),Z[i]=read();
  160.             if (p[i]==1)
  161.                 Hash(X[i]-Y[i]+n,X[i]+Y[i]);
  162.         }
  163.         sort(q+1,q+k);
  164.         k=unique(q+1,q+k)-q;
  165.         for (int i=1;i<=m;++i){
  166.             if (p[i]==1)
  167.                 add(X[i]-Y[i]+n,X[i]+Y[i],Z[i]);
  168.             else {
  169.                 int a=max(1,X[i]-Y[i]+n-Z[i]);
  170.                 int b=max(1,X[i]+Y[i]-Z[i]);
  171.                 int c=min(ma,X[i]-Y[i]+n+Z[i]);
  172.                 int d=min(ma,X[i]+Y[i]+Z[i]);
  173.                 printf("%d\n",sum(c,d)-sum(c,b-1)-sum(a-1,d)+sum(a-1,b-1));
  174.             }
  175.         }
  176.     }
  177.     return 0;
  178. }

 [HEOI2016/TJOI2016]序列

  1. #include<bits/stdc++.h>
  2. #define N 100010
  3. using namespace std;
  4. int a[N],maxx[N],minn[N],fa[N],c[N],f[N];
  5. int num,ans,x,y,n,m;
  6. inline int read() {
  7.     char ch;
  8.     bool f=false;
  9.     int res=0;
  10.     while (((ch=getchar())<'0'||ch>'9')&&ch!='-');
  11.     if (ch=='-')
  12.         f=true;
  13.     else
  14.         res=ch-'0';
  15.     while ((ch=getchar())>='0'&&ch<='9')
  16.         res=(res<<3)+(res<<1)+ch-'0';
  17.     return f?~res+1:res;
  18. }
  19. inline int lowbit(int x){
  20.     return x&(-x);
  21. }
  22. inline void add(int x,int y){
  23.     while (x<=num){
  24.         if (y)
  25.             c[x]=max(c[x],y);
  26.         else c[x]=0;
  27.         x+=lowbit(x);
  28.     }
  29. }
  30. inline int sum(int x){
  31.     int summ=0;
  32.     while (x>0){
  33.         summ=max(summ,c[x]),x-=lowbit(x);
  34.     }
  35.     return summ;
  36. }
  37. inline bool cmp(int x,int y){
  38.     return a[x]<a[y];
  39. }
  40. inline bool cmp1(int x,int y){
  41.     return minn[x]<minn[y];
  42. }
  43. void sc(int l,int r){
  44.     if (l==r){
  45.         if (f[l]<1)
  46.             f[l]=1;
  47.         return;
  48.     }
  49.     int mid=(l+r)>>1;
  50.     sc(l,mid);
  51.     for (int i=l;i<=r;++i)
  52.         fa[i]=i;
  53.     sort(fa+l,fa+1+mid,cmp);
  54.     sort(fa+mid+1,fa+1+r,cmp1);
  55.     for (int i=mid+1,j=l;i<=r;++i){
  56.         while (j<=mid&&a[fa[j]]<=minn[fa[i]])
  57.             add(maxx[fa[j++]],f[fa[j]]);
  58.         f[fa[i]]=max(f[fa[i]],sum(a[fa[i]])+1);
  59.         if (i==r)
  60.             for (int k=l;k<j;++k)
  61.                 add(maxx[fa[k]],0);
  62.     }
  63.     sc(mid+1,r);
  64. }
  65. int main(){
  66.     n=read(),m=read();
  67.     for (int i=1;i<=n;++i)
  68.         a[i]=read();
  69.     for (int i=1;i<=n;++i)
  70.         maxx[i]=minn[i]=a[i];
  71.     while (m--){
  72.         x=read(),y=read();
  73.         maxx[x]=max(maxx[x],y);
  74.         minn[x]=min(minn[x],y);
  75.     }
  76.     for (int i=1;i<=n;++i)
  77.         num=max(num,a[i]);
  78.     sc(1,n);
  79.     for (int i=1;i<=n;++i)
  80.         ans=max(ans,f[i]);
  81.     printf("%d\n",ans);
  82.     return 0;
  83. }

[NOI2007]货币兑换

题意不想描述，写的很清晰了。。。

Solution：

　　好像隐隐约约中听大佬说不止要用cdq，还要用斜率优化、splay维护....感觉好可怕

　　然后就搁置了两三天...（颓了两三天）

　　设f[i]表示前i天的最大收益。

　　第i天可以换成B券最大数目f[i]*(1/(Rate[i]*A[i]+B[i]))

　　第i天可以换成的A券最大数目f[i]*(Rate[i]/(Rate[i]*A[i]+B[i]))

　　第i天将第j天的券全卖掉A[i]*X(j)+B[i]*Y(j)

　　所以f[i]=max{f[i-1],A[i]*X(j)+B[i]*Y(j)}

　　则我们需要求 max p=A[i]*X(j)+B[i]*Y(j)

　　即我们要最大化直线方程Y(j)=-A[i]/B[i]*X(j)+p/B[j]的截距

  1. #ifndef _GLIBCXX_NO_ASSERT
    #include <cassert>
  2. #endif
  3. #include <cctype>
  4. #include <cerrno>
  5. #include <cfloat>
  6. #include <ciso646>
  7. #include <climits>
  8. #include <clocale>
  9. #include <cmath>
  10. #include <csetjmp>
  11. #include <csignal>
  12. #include <cstdarg>
  13. #include <cstddef>
  14. #include <cstdio>
  15. #include <cstdlib>
  16. #include <cstring>
  17. #include <ctime>
  18.  
  19. #if __cplusplus >= 201103L
  20. #include <ccomplex>
  21. #include <cfenv>
  22. #include <cinttypes>
  23. #include <cstdalign>
  24. #include <cstdbool>
  25. #include <cstdint>
  26. #include <ctgmath>
  27. #include <cwchar>
  28. #include <cwctype>
  29. #endif
  30.  
  31. // C++
  32. #include <algorithm>
  33. #include <bitset>
  34. #include <complex>
  35. #include <deque>
  36. #include <exception>
  37. #include <fstream>
  38. #include <functional>
  39. #include <iomanip>
  40. #include <ios>
  41. #include <iosfwd>
  42. #include <iostream>
  43. #include <istream>
  44. #include <iterator>
  45. #include <limits>
  46. #include <list>
  47. #include <locale>
  48. #include <map>
  49. #include <memory>
  50. #include <new>
  51. #include <numeric>
  52. #include <ostream>
  53. #include <queue>
  54. #include <set>
  55. #include <sstream>
  56. #include <stack>
  57. #include <stdexcept>
  58. #include <streambuf>
  59. #include <string>
  60. #include <typeinfo>
  61. #include <utility>
  62. #include <valarray>
  63. #include <vector>
  64.  
  65. #if __cplusplus >= 201103L
  66. #include <array>
  67. #include <atomic>
  68. #include <chrono>
  69. #include <condition_variable>
  70. #include <forward_list>
  71. #include <future>
  72. #include <initializer_list>
  73. #include <mutex>
  74. #include <random>
  75. #include <ratio>
  76. #include <regex>
  77. #include <scoped_allocator>
  78. #include <system_error>
  79. #include <thread>
  80. #include <tuple>
  81. #include <typeindex>
  82. #include <type_traits>
  83. #include <unordered_map>
  84. #include <unordered_set>
  85. #endif
  86. using namespace std;
  87. const int MAXN=100010;
  88. const double inf=1e20,Inf=1e-8;
  89. struct Node {
  90.     double x,y,a,b,k,z;
  91.     int num;
  92. }p[MAXN],tt[MAXN];
  93. int a[MAXN],n,cnt;
  94. double f[MAXN];
  95. inline bool cmp(Node x,Node y) {
  96.     return x.k>y.k;
  97. }
  98. inline int read() {
  99.     char ch;
  100.     bool f=false;
  101.     int res=0;
  102.     while (((ch=getchar())<'0'||ch>'9')&&ch!='-');
  103.     if (ch=='-')
  104.         f=true;
  105.     else
  106.         res=ch-'0';
  107.     while ((ch=getchar())>='0'&&ch<='9')
  108.         res=(res<<3)+(res<<1)+ch-'0';
  109.     return f?~res+1:res;
  110. }
  111. inline double pd(int a,int b){
  112.     if (!b)
  113.         return -inf;
  114.     else
  115.         if (fabs(p[a].x-p[b].x)<Inf)
  116.             return inf;
  117.     else
  118.         return (p[b].y-p[a].y)/(p[b].x-p[a].x);
  119. }
  120. void sc(int t,int w){
  121.     int mid=(t+w)>>1;
  122.     if (t==w){
  123.         if (f[t-1]>f[t])
  124.             f[t]=f[t-1];
  125.         p[t].y=f[t]/(p[t].a*p[t].z+p[t].b);
  126.         p[t].x=p[t].y*p[t].z;
  127.         return;
  128.     }
  129.     int j=1,cnt1=t,cnt2=mid+1;
  130.     for (int i=t;i<=w;++i)
  131.         if (p[i].num<=mid)
  132.             tt[cnt1++]=p[i];
  133.         else
  134.             tt[cnt2++]=p[i];
  135.     for (int i=t;i<=w;++i)
  136.         p[i]=tt[i];
  137.     sc(t,mid);
  138.     cnt=0;
  139.     for (int i=t;i<=mid;++i){
  140.         while (cnt>1&&pd(a[cnt-1],a[cnt])<pd(a[cnt-1],i)+Inf)
  141.             cnt--;
  142.         a[++cnt]=i;
  143.     }
  144.     for(int i=mid+1;i<=w;i++) {
  145.         while(j<cnt&&pd(a[j],a[j+1])+Inf>p[i].k)
  146.             j++;
  147.         f[p[i].num]=max(f[p[i].num],p[a[j]].x*p[i].a+p[a[j]].y*p[i].b);
  148.     }
  149.     sc(mid+1,w);
  150.     int t1=t,ww1=mid+1;
  151.     for(int i=t;i<=w;i++) {
  152.         if(((p[t1].x<p[ww1].x||(fabs(p[t1].x-p[ww1].x)<Inf&&p[t1].y<p[ww1].y))||ww1>w)&&t1<=mid)
  153.             tt[i]=p[t1++];
  154.         else
  155.             tt[i]=p[ww1++];
  156.     }
  157.     for(int i=t;i<=w;i++) p[i]=tt[i];
  158. }
  159. int main(){
  160.     n=read();
  161.     scanf("%lf",&f[0]);
  162.     for (int i=1;i<=n;++i){
  163.         scanf("%lf%lf%lf",&p[i].a,&p[i].b,&p[i].z);
  164.         p[i].k=-p[i].a/p[i].b;
  165.         p[i].num=i;
  166.     }
  167.     sort(p+1,p+n+1,cmp);
  168.     sc(1,n);
  169.     printf("%.3lf\n",f[n]);
  170.     return 0;
  171. }

gcd论文！

cdq分治（hdu 5618 Jam's problem again[陌上花开]、CQOI 2011 动态逆序对、hdu 4742 Pinball Game、hdu 4456 Crowd、[HEOI2016/TJOI2016]序列、[NOI2007]货币兑换 ）的更多相关文章

  1. HDU 5618 Jam's problem again（三维偏序，CDQ分治，树状数组，线段树）

    Jam's problem again Time Limit: 5000/2500 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Othe ...

  2. HDU 5618 Jam's problem again CDQ分治 BC ROUND 70

    题意:给你1e5个点(x,y,z),对于每一个点询问有多少个点(x1,y1,z1)满足x1<=x&&y1<=y&&z1<=z 分析:(官方题解奉上)很 ...

  3. HDU 5618 Jam's problem again (cdq分治+BIT 或 树状数组套Treap)

    题意:给n个点,求每一个点的满足 x y z 都小于等于它的其他点的个数. 析:三维的,第一维直接排序就好按下标来,第二维按值来,第三维用数状数组维即可. 代码如下: cdq 分治: #pragma ...

  4. HDU 5618 Jam's problem again

    题意: 三维坐标,对于1个点,找出有多少个点,3个坐标都比该点小! Sample Input 1 4 10 4 7 10 6 6 8 2 5 7 3 10   Sample Output 1 1 0 ...

  5. bzoj 3295 动态逆序对 （三维偏序，CDQ+树状数组）

    链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3295 思路: 可以将这道题看成倒着插入,这样就可以转化成求逆序对数,用CDQ分治降维,正反用 ...

  6. 【Luogu1393】动态逆序对（CDQ分治）

    [Luogu1393]动态逆序对(CDQ分治) 题面 题目描述 对于给定的一段正整数序列,我们定义它的逆序对的个数为序列中ai>aj且i < j的有序对(i,j)的个数.你需要计算出一个序 ...

  7. [BZOJ3295][Cqoi2011]动态逆序对 CDQ分治&树套树

    3295: [Cqoi2011]动态逆序对 Time Limit: 10 Sec  Memory Limit: 128 MB Description 对于序列A,它的逆序对数定义为满足i<j,且 ...

  8. BZOJ 3295 动态逆序对 | CDQ分治

    BZOJ 3295 动态逆序对 这道题和三维偏序很类似.某个元素加入后产生的贡献 = time更小.pos更小.val更大的元素个数 + time更小.pos更大.val更小的元素个数. 分别用类似C ...

  9. 【算法】CDQ分治 -- 三维偏序 & 动态逆序对

    初次接触CDQ分治,感觉真的挺厉害的.整体思路即分而治之,再用之前处理出来的答案统计之后的答案. 大概流程是(对于区间 l ~ r): 1.处理 l ~mid, mid + 1 ~ r 的答案: 2. ...

随机推荐

  1. Android开发如何轻松实现基于Tesseract的Android OCR应用程序

    介绍 此应用程序使用Tesseract 3的Tesseract OCR引擎,该引擎通过识别字符模式( https://github.com/tesseract-ocr/tesseract )来工作. ...

  2. 南京邮电大学java程序设计作业在线编程第七次作业

    王利国的"Java语言程序设计第7次作业(2018)"详细 主页 我的作业列表 作业结果详细 总分:100 选择题得分:60  1. 下列叙述中,错误的是( ). A.Java中, ...

  3. UE3客户端加入DS过程

    拉起DS进程 客户端将比赛地图及相关参数发送给ZoneSvr请求开赛,收到消息后,ZoneSvr会分配一个ip和端口号,并与客户端发过来的地图及其他参数,来构建一个命令行来拉起一个DS进程, DS启动 ...

  4. Linux 环境 Java JDK 安装&基本配置

    索引: 目录索引 参看代码 GitHub: jdk.txt 一.Linux (DeepinOS) 环境 1.官网下载 jdk-8u112-linux-x64.tar.gz 2.创建目录 mkdir - ...

  5. windows下安装MongoDB扩展和配置

    windows下安装MongoDB扩展和配置 1.下载mongoDB扩展,根据当前php版本进行下载 地址如下:http://pecl.php.net/package/mongo 我本地php版本是 ...

  6. windows 为qt5.7.1 安装openssl

    本人使用qt5.7.1+msvc2015写一个https的客户端程序,但是用到解析https协议时,报出如下错误 qt.network.ssl: QSslSocket: cannot call unr ...

  7. LeetCode算法题-Range Addition II（Java实现）

    这是悦乐书的第271次更新,第285篇原创 01 看题和准备 今天介绍的是LeetCode算法题中Easy级别的第138题(顺位题号是598).给定一个m行n列的新二维数组M,其初始值为0.提供一个二 ...

  8. java类加载及类初始化

    1.前言 java是跨平台语言,主要是因为它的java虚拟机的存在,java有事编译语言,所以需要将编写的java文件编译成jvm可运用的class字节码文件.在java中一切皆对象.对于Java虚拟 ...

  9. Spring类型转换（Converter）

    Spring的类型转换 以前在面试中就有被问到关于spring数据绑定方面的问题,当时对它一直只是朦朦胧胧的概念,最近稍微闲下来有时间看了一下其中数据转换相关的内容,把相应的内容做个记录. 下面先说明 ...

  10. uml类图关系

    原文地址http://www.jfox.info/uml-lei-tu-guan-xi-fan-hua-ji-cheng-shi-xian-yi-lai-guan-lian-ju-he-zu-he 在 ...