Codeforces Round #563 (Div. 2) E. Ehab and the Expected GCD Problem

https://codeforces.com/contest/1174/problem/E

dp

好题

*(if 满足条件)

满足条件 *1

不满足条件 *0

///这代码虽然写着方便，但是常数有点大

 #include <cstdio>
 #include <cstdlib>
 #include <cmath>
 #include <cstring>
 #include <string>
 #include <algorithm>
 #include <iostream>
 using namespace std;
 #define ll long long

 const double eps=1e-;
 const ll inf=1e9;
 const ll mod=1e9+;
 const int maxn=1e6+;

 /*
 dp题
 检查：
 1.认真检查公式
 2.认真检查代码
 3.造一个规模适中的数据
 手动推导
 对应程序的结果
 判断是否相同
 */

 ///+1 20->21 2->3
 ///0->-1 前面加一个辅助数组，这个数组里的所有元素数值为0
 ///x1*y1+x2*y2+x3*y3+...x20*y20 有可能会超long long, 所以所有乘法运算后面都加%mod
 ///1e18*k k<=9
 int f[maxn][21][3],v[21][3];

 int main()
 {
     int n,two,three,i,j,k;
     scanf("%d",&n);
     two=log(n+eps)/log();///+eps
     three=(1.0*n/(<<two)>=1.5);
     for (i=;i<=two;i++)
         for (j=;j<=three;j++)
             v[i][j]=n/((<<i)*(j==?:));

     f[][two][]=;
     if (three==)
         f[][two-][]=;
     for (i=;i<=n;i++)
         for (j=;j<=two;j++)
             for (k=;k<=three;k++)
                 ///乘1ll,最后强制转换long long 转 int
                 f[i][j][k]=( 1ll*f[i-][j][k]*(v[j][k]-(i-))*((v[j][k]-(i-))>=)%mod + 1ll*f[i-][j+][k]*(v[j][k]-v[j+][k])*(j!=two)%mod + 1ll*f[i-][j][k+]*(v[j][k]-v[j][k+])*(k!=three)%mod )%mod;
     printf("%d",f[n][][]);
     return ;
 }

推荐比如说这个xiongdi的代码

https://codeforces.com/contest/1174/submission/56391704

超内存的代码

 #include <cstdio>
 #include <cstdlib>
 #include <cmath>
 #include <cstring>
 #include <string>
 #include <algorithm>
 #include <iostream>
 using namespace std;
 #define ll long long

 const double eps=1e-;
 const ll inf=1e9;
 const ll mod=1e9+;
 const int maxn=1e6+;

 /*
 dp题
 检查：
 1.认真检查公式
 2.认真检查代码
 3.造一个规模适中的数据
 手动推导
 对应程序的结果
 判断是否相同
 */

 ///+1 20->21 2->3
 ///0->-1 前面加一个辅助数组，这个数组里的所有元素数值为0
 ///x1*y1+x2*y2+x3*y3+...x20*y20 有可能会超long long, 所以所有乘法运算后面都加%mod
 ///1e18*k k<=9
 ll f[maxn][21][3],v[21][3];

 int main()
 {
     int n,two,three,i,j,k;
     scanf("%d",&n);
     two=log(n+eps)/log();///+eps
     three=(1.0*n/(<<two)>=1.5);
     for (i=;i<=two;i++)
         for (j=;j<=three;j++)
             v[i][j]=n/((<<i)*(j==?:));

     f[][two][]=;
     if (three==)
         f[][two-][]=;
     for (i=;i<=n;i++)
         for (j=;j<=two;j++)
             for (k=;k<=three;k++)
                 f[i][j][k]=( f[i-][j][k]*(v[j][k]-(i-))*((v[j][k]-(i-))>=)%mod + f[i-][j+][k]*(v[j][k]-v[j+][k])*(j!=two)%mod + f[i-][j][k+]*(v[j][k]-v[j][k+])*(k!=three)%mod )%mod;

     printf("%lld",f[n][][]);
     return ;
 }