首页 > 其他分享> 文章详细

「APIO2019」奇怪装置

2022-05-02 19:31:54 阅读：220 来源： 互联网

标签：node 装置 frac ll bmod APIO2019 ans cases 奇怪

Problem

传送门：左转右转

Solution

至高无上的xxx250说这道题有循环节，且观察样例可大胆猜测周期一定是\(AB\)的因数。

题目要求求这一坨玩意：

\[\begin{cases} x=((t+\lfloor\frac{t}{B}\rfloor)\ \bmod A\\ y=(t\ \bmod B) \end{cases}\]

设循环节为\(k\)，因为\(t=0\)时，\(\begin{cases} x=0 \\ y=0 \end{cases}\)，所以，\(\begin{cases} x = ((k + \frac{k}{B})\ \bmod A) = 0\\ y = (k \ \bmod B) = 0 \end{cases}\)

即\(\begin{cases} 0 = (k + \frac{k}{B}) \ \bmod A \\ 0 = (k \ \bmod B) \end{cases}\)

因为\(0 = (k \ \bmod B)\) 所以式子中的向下取整可以去掉。\(0 = (k + \frac{k}{B}) \ \bmod A\)于是可化为\(0 = \frac{k(B+1)}{B} \ \bmod A\)

所以\(0 \equiv \frac{k(B+1)}{B} \ \bmod A\)

式子可化为\(0 \equiv \frac{k}{B} \ (\bmod \ \frac{A}{gcd(A,B+1)})\)
移项可得\(0 \equiv k \ (\bmod \ \frac{AB}{gcd(A,B+1)})\)

所以最小循环节为\(k(\bmod \ \frac{AB}{gcd(A,B+1)}) = 0\)，所以\(k = \frac{AB}{gcd(A,B+1)}\)

然后题目所求可以转化为区间覆盖问题。

分类讨论亿一下，区间\([l,r]\)

若\(k \leq r - l + 1\)，则答案为\(k\)

否则，设\(L = l \bmod k , R = r \bmod k\)。将点对映射到数轴上，将其转化为一个线段覆盖问题。（需要分别考虑\(L \leq R\)和\(L > R\)的情况）

\(code:\)

//最小循环节 + 线段覆盖
#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

const int MAX_N = 1000000 + 5;

int n,m;
ll A,B,k,ans;
struct node{
    ll l,r;
}s[MAX_N];
inline bool cmp(node a,node b){
    return a.l<b.l;
}
int main(){
    scanf("%d%lld%lld",&n,&A,&B);
    k=A/__gcd(A,B+1)*B;
    if(k<0) k=2e18;
    for(int i=1;i<=n;++i){
        ll l,r;
        scanf("%lld%lld",&l,&r);
        if(r-l+1>=k){
            ans=k;
            continue;
        }
        else{
            l%=k,r%=k;
            if(l>r) s[++m]=(node){0,r},s[++m]=(node){l,k-1};
            else s[++m]=(node){l,r};
        }
    }
    if(ans){
        printf("%lld\n",ans);
        return 0;
    }
    sort(s+1,s+1+m,cmp);
    ans=0;
    ll L=s[1].l,R=s[1].r;
    s[++m]=(node){k,k};
    for(int i=2;i<=m;++i){
        if(R<s[i].l){
            ans+=R-L+1;
            L=s[i].l,R=s[i].r;
        }
        else{
            R=max(s[i].r,R);
        }
    }
    printf("%lld\n",ans);
    return 0;
}

标签：node,装置,frac,ll,bmod,APIO2019,ans,cases,奇怪
来源： https://www.cnblogs.com/TheAutumnGlory/p/16216481.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

「APIO2019」奇怪装置

Problem

Solution