2021牛客暑期多校训练营5 G-Greater Integer, Better LCM（暴力+子集dp）

G-Greater Integer, Better LCM

看到校大佬的代码就瞬间会做了。

当时和队友想的是，首先转化题意即找两个数$x\ge a,y\ge b[\text{lcm}(a,b)=c]$，首先不考虑$x\ge a,y\ge b$的限制，显然对于$c=p_1^{q_1}{p}_2^{q_2}\dots p_n^{q_n}$，考虑其中一项$p_i^{q_i}$，对于该项显然要么$[p_i^0{]}_a,[p_i^{q_i}{]}_b$要么$[p_i^{q_i}{]}_a,[p_i^0{]}_b$，这样只需要暴搜就行了，时间复杂度为$O(2^n)$，但是加上了这个限制就不知道怎么搞了。

题目中有一个非常重要的地方就是$\sum q_i\le 18$，这意味这我们不需要考虑是$[p_i^0{]}_a,[p_i^{q_i}{]}_b$还是$[p_i^{q_i}{]}_a,[p_i^0{]}_b$，只需要暴力枚举每一个$q_i$的次幂最后只针对$x\ge a,y\ge b$考虑即可。 详细看代码

最后一部分就是枚举超集dp转化一下。

#include<bits/stdc++.h>
using namespace std;
using ll=long long;
template <class T=int> T rd()
{
    T res=0;T fg=1;
    char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') fg=-1;ch=getchar();}
    while( isdigit(ch)) res=(res<<1)+(res<<3)+(ch^48),ch=getchar();
    return res*fg;
}
void print(__int128 v)
{
    if(v==0) return puts("0"),void();
    string ans;
    while(v)
    {
        ans.push_back(v%10+'0');
        v/=10;
    }
    reverse(ans.begin(),ans.end());
    printf("%s\n",ans.c_str());
}
ll p[25],q[25];
int n;
__int128 a,b;
__int128 fa[(1<<18)+5],fb[(1<<18)+5];
void dfs(int u,int s,__int128 cur)
{
    if(cur>=a) fa[s]=min(fa[s],cur);
    if(cur>=b) fb[s]=min(fb[s],cur);
    
    if(u>=n) return;
    
    dfs(u+1,s,cur);
    for(int i=1;i<q[u];i++) 
    {
        cur*=p[u];
        dfs(u+1,s,cur);
    }
    cur*=p[u];
    dfs(u+1,s^(1<<u),cur);
}
int main()
{
    
    n=rd();
    for(int i=0;i<n;i++) p[i]=rd<ll>(),q[i]=rd<ll>();
    
    a=rd<__int128>(),b=rd<__int128>();
    
    memset(fa,0x3f,sizeof fa);
    memset(fb,0x3f,sizeof fb);
    
    dfs(0,0,1);
    
    __int128 ans; memset(&ans,0x3f,sizeof ans);
    
    for(int i=(1<<n)-1;i>=0;i--)
        for(int j=0;j<n;j++)
            if((i>>j)&1)
            {
                fa[i^(1<<j)]=min(fa[i^(1<<j)],fa[i]);
                fb[i^(1<<j)]=min(fb[i^(1<<j)],fb[i]);
            }
    for(int i=0;i<(1<<n);i++) ans=min(ans,fa[i]+fb[((1<<n)-1)^i]-a-b);
    print(ans);
    return 0;
}