Luogu P5444 [APIO2019]奇怪装置

题目
这种题目看上去就是有循环节的对吧。
在考场上,一个可行的方式是打表。
现在我们手推一下这个循环节。
记函数\(f(t)=(((t+\lfloor\frac tB\rfloor)\%A),(t\%B))\),那么\(f(t_1)=f(t_2)\)的充要条件为:
\[ t_1+\lfloor\frac {t_1}B\rfloor\equiv t_2+\lfloor\frac {t_2}B\rfloor(mod\ A) \]
\[ t_1\equiv t_2(mod\ B) \]
看到第二个很熟悉的式子,我们不妨设\(t_2=t_1+xB(x>0)\),代入第一个式子得到
\[ t_1+\lfloor\frac {t_1}B\rfloor\equiv t_1+xB+\lfloor\frac {t_1+xB}B\rfloor(mod\ A) \]
后面那坨东西我们把它化简一下,拆开和左边的抵掉就变成了
\[ xB+x\equiv0(mod\ A) \]

\[ A|x(B+1) \]
\(B+1\)移到左边去
\[ \frac{A(B+1)}{(A,B+1)}|x \]
也就是说最小循环节为\(\frac{A(B+1)}{(A,B+1)}\)
剩下的做一个简单的区间覆盖即可。
区间覆盖的话左端点升序排序,记录最大右端点即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
namespace IO
{
    char ibuf[(1<<21)+1],*iS,*iT;
    char Get(){return (iS==iT? (iT=(iS=ibuf)+fread(ibuf,1,(1<<21)+1,stdin),(iS==iT? EOF:*iS++)):*iS++);}
    ll read(){ll x=0;int c=Get();while(!isdigit(c))c=Get();while(isdigit(c))x=x*10+c-48,c=Get();return x;}
}
using namespace IO;
const int N=1000007;
struct node{ll l,r;}a[N<<1];
int operator<(node a,node b){return a.l<b.l;}
ll max(ll a,ll b){return a>b? a:b;}
int main()
{
    int n=read(),m=0,i;ll A=read(),B=read(),M=A/__gcd(A,B+1)*B,l,r,ans=0,mx=0;
    for(i=1;i<=n;++i)
    {
    l=read(),r=read();
    if(r-l+1>=M) return !printf("%lld",M);
    l%=M,r%=M;
    if(l>r) a[++m]=(node){0,r},a[++m]=(node){l,M-1}; else a[++m]=(node){l,r};
    }
    sort(a+1,a+m+1);
    for(i=1;i<=m;++i)
    {
    if(a[i].l>mx) ans+=a[i].l-mx;
    mx=max(mx,a[i].r+1);
    }
    if(mx<M) ans+=M-mx;
    return !printf("%lld",M-ans);
}

Luogu P5444 [APIO2019]奇怪装置

上一篇:net core WebApi——缓存神器Redis


下一篇:wpf source path