AtCoder Regular Contest 111

A Simple Math 2

$$\frac{10^n}{m}\mathrm{m}\mathrm{o}\mathrm{d}m=\frac{10^n\mathrm{m}\mathrm{o}\mathrm{d}{m}^2+k\ast m^2}{m}\mathrm{m}\mathrm{o}\mathrm{d}m$$
Meanwhile,$$\frac{k\ast m^2}{m}\mathrm{m}\mathrm{o}\mathrm{d}m=0$$
Then the answer is $$\frac{10^n\mathrm{m}\mathrm{o}\mathrm{d}{m}^2}{m}\mathrm{m}\mathrm{o}\mathrm{d}m$$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll n,m;
ll qpow(ll a,ll n,ll mod)
{
    ll ans=1;
    for(;n;n>>=1,a=a*a%mod)
        if(n&1) ans=ans*a%mod;
    return ans;
}
int main()
{
    cin>>n>>m;
    int r=qpow(10,n,m);
    int ans=qpow(10,n,m*m);
    ans=ans/m%m;
    printf("%d\n",ans);
}

B Reversible Cards

2020牛客多校I。
https://ac.nowcoder.com/acm/contest/5673/I

C Too Heavy

$i$现在背的包是$p_j$的包，可以得到一个环$i\to p_i\to p_{p_i}\to ...\to i$。如果环长度只有$1$，略过。
只要环中有一个$b_{p_i}\ge a_i$的，答案就是$-1$。
否则，令环的序列为$i_1,i_2,i_3,..,i_k$，有$$max(a_{i_1},..,a_{i_k})>max(b_{i_1},..,b_{i_k})$$
只需要每次和$a$最大的那个位置交换即可。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=2e5+5;
int n,a[N],b[N],p[N];
bool vis[N];
vector<int>v;
vector<pair<int,int>>ans;
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);
    for(int i=1;i<=n;i++) scanf("%d",&b[i]);
    for(int i=1;i<=n;i++) scanf("%d",&p[i]);
    for(int i=1;i<=n;i++)
        if(!vis[i])
    {
        vis[i]=true;
        if(p[i]==i) continue;
        if(b[p[i]]>=a[i])
        {
            printf("-1\n");
            return 0;
        }
        v.clear();
        v.push_back(i);
        vis[i]=true;
        for(int j=p[i];j!=i;j=p[j])
        {
            if(b[p[j]]>=a[j])
            {
                printf("-1\n");
                return 0;
            }
            v.push_back(j);
            vis[j]=true;
        }
        int pos=0;
        for(int i=1;i<v.size();i++)
            if(a[v[i]]>a[v[pos]]) pos=i;
        for(int i=pos+1;i<v.size();i++)
            ans.push_back({v[i],v[pos]});
        for(int i=0;i<pos;i++)
            ans.push_back({v[i],v[pos]});
    }
    printf("%d\n",ans.size());
    for(pair<int,int>x:ans)
        printf("%d %d\n",x.first,x.second);
}

D Orientation

一条边$u,v$，如果$c_v\ne c_u$，连边显然。否则，它们在同一个强连通分量里，只需要模仿$tarjan$进行一次$dfs$即可，因为数据保证有解。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=105;
int n,m,a[N*N],b[N*N],c[N],mp[N][N];
string ans[N*N];
bool vis[N];
void dfs(int u)
{
    vis[u]=true;
    for(int i=1;i<=n;i++)
    {
        if(mp[u][i])
        {
            mp[i][u]=false;
            if(!vis[i]) dfs(i);
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&a[i],&b[i]);
    }
    for(int i=1;i<=n;i++) scanf("%d",&c[i]);
    for(int i=1;i<=m;i++)
        if(c[a[i]]>c[b[i]])
        ans[i]="->";
        else if(c[b[i]]>c[a[i]])
        ans[i]="<-";
        else mp[a[i]][b[i]]=mp[b[i]][a[i]]=true;
    for(int i=1;i<=n;i++)
        if(!vis[i]) dfs(i);
    for(int i=1;i<=m;i++)
        if(mp[a[i]][b[i]]) ans[i]="->";
    else if(mp[b[i]][a[i]]) ans[i]="<-";
    for(int i=1;i<=m;i++)
        cout<<ans[i]<<'\n';
}

E Simple Math 3

显然问题等价于求：
$\sum _{i=1}^k[\lfloor \frac{(A-1+B\ast i)}{D}\rfloor =\lfloor \frac{(A+C\ast i)}{D}\rfloor ]$，当$k>\lfloor \frac{D-2}{C-B}\rfloor$两边一定不相等，因为它们直接的距离相差大于等于$D$。
否则，注意到$\lfloor \frac{(A-1+B\ast i)}{D}\rfloor \le \lfloor \frac{(A+C\ast i)}{D}\rfloor$，且当$k\le \lfloor \frac{D-2}{C-B}\rfloor$，他们的差最多为$1$。也就是差为$1$的个数等于不相等的个数。
注意到这个，答案就是$\sum _{i=1}^k\lfloor \frac{(A+C\ast i)}{D}\rfloor -\lfloor \frac{(A-1+B\ast i)}{D}\rfloor$。
这是类欧几里得算法模板。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll A,B,C,D;
const int mod=998244353;
ll n,a,b,c;
struct node
{
    ll f,g,h;
    node(ll f=0,ll g=0,ll h=0):f(f),g(g),h(h){}
};
ll qpow(ll a,ll n)
{
    ll ans=1;
    for(;n;n>>=1,a=a*a%mod)
        if(n&1) ans=ans*a%mod;
    return ans;
}
const ll inv6=qpow(6,mod-2),inv2=mod-mod/2;
ll hs(ll a,ll b){return a*2000000000+b;}
ll s2(ll x){return x*(x+1)/2%mod;}
ll s3(ll x){return x*(x+1)%mod*(2*x+1)%mod*inv6%mod;}
node F(ll a,ll b,ll c,ll n)
{
    if(c==0||a==0&&b<c) return node(0,0,0);
    if(n==0) return node(b/c,0,b/c*(b/c)%mod);
    node ans;
    if(a>=c||b>=c)
    {
        node res=F(a%c,b%c,c,n);
        ans.f=(a/c*s2(n)%mod+b/c*(n+1)%mod+res.f)%mod;
        ans.g=(res.g+a/c*s3(n)%mod+b/c*s2(n)%mod)%mod;
        ans.h=(res.h+a/c*2%mod*res.g+b/c*2%mod*res.f)%mod;
        ans.h=(ans.h+(a/c)*(a/c)%mod*s3(n)+(b/c)*(b/c)%mod*(n+1))%mod;
        ans.h=(ans.h+(a/c)*(b/c)%mod*2%mod*s2(n))%mod;
    }
    else
    {
        ll m=(a*n+b)/c;
        node res=F(c,c-b-1,a,m-1);
        ans.f=(n*m%mod-res.f+mod)%mod;
        ans.g=(s2(n)*m%mod-res.h*inv2%mod-res.f*inv2%mod+mod+mod)%mod;
        ans.h=n*m%mod*(m+1)%mod-2*res.g%mod-2*res.f%mod-ans.f;
        ans.h=(ans.h%mod+mod)%mod;
    }
    return ans;
}
int main()
{
    int T;scanf("%d",&T);
    while(T--)
    {
        scanf("%lld%lld%lld%lld",&A,&B,&C,&D);
        int up=(D-2)/(C-B);
        node ans1=F(B,A-1,D,up),ans2=F(C,A,D,up);
        ll ans=(ans2.f-ans1.f+mod)%mod;
        if(A/D!=(A-1)/D) ans--;
        ans=up-ans;
        printf("%lld\n",ans);
    }
}

F Do you like query problems?

$ans=期望值\ast 方案数$，因此只计算其期望$E$即可。
对于一个固定的位置$i$，其被操作的概率为$P_i=\frac{i\ast (n-i+1)\ast 2}{n\ast (n+1)}$，只需解决以下问题即可：
给定两个数$x,su{m}_i$初始都为$0$，随机进行$q$次操作，每次操作：

1. 有$1-P_i$概率什么都不做
2. 有$P_i\ast \frac{1}{2\ast m+1}$的概率使得$su{m}_i=su{m}_i+x$
3. 有$P_i\ast \frac{1}{2\ast m+1}$的概率使得 x = m i n ( x , v ) ( v ∈ { 0 , 1 , . . . , m − 1 } ) x=min(x,v)(v\in\{0,1,...,m-1\}) x=min(x,v)(v∈{0,1,...,m−1})
4. 有$P_i\ast \frac{1}{2\ast m+1}$的概率使得 x = m a x ( x , v ) ( v ∈ { 0 , 1 , . . . , m − 1 } ) x=max(x,v)(v\in\{0,1,...,m-1\}) x=max(x,v)(v∈{0,1,...,m−1})

求$su{m}_i$的期望值是多少。
设$p_{j,v}$为第$j$次操作后$x=v$的概率。$p_{j,\ge v}$为第$j$次操作后$x\ge v$的概率。则$$\begin{gathered} su{m}_i=\sum _{j=0}^{q-1}\sum _{v=0}^{m-1}{p}_{j,v}\ast v\ast P_i\ast \frac{1}{2\ast m+1} \\ =\sum _{j=0}^{q-1}\sum _{v=1}^{m-1}{p}_{j,\ge v}\ast P_i\ast \frac{1}{2\ast m+1} \\ {p}_{j+1,\ge v}=p_{j,\ge v}\ast (1-P_i)+P_i\ast p_{j,\ge v}\ast \frac{m+1+m-v}{2\ast m+1}+ \\ {P}_i\ast (1-p_{j,\ge v})\ast \frac{m-v}{2\ast m+1}= \\ {p}_{j,\ge v}\ast (1-\frac{m}{2\ast m+1}\ast P_i)+P_i\ast \frac{m-v}{2\ast m+1} \end{gathered}$$
令$a_i=1-\frac{m}{2\ast m+1}\ast P_i,b_v=\frac{m-v}{2\ast m+1}$，可以得到$$\begin{gathered} p_{0,\ge v}=0(v\ge 1) \\ {p}_{j+1,\ge v}=p_{j,\ge v}\ast a_i+P_i\ast b_v \end{gathered}$$
观察该式子可以发现$p_{j,\ge v}=P_i\ast b_v\ast \sum _{k=0}^{j-1}{a}_i^j(j\ge 1)$。
则可以得到$$\begin{gathered} E=\sum _{i=1}^{n}su{m}_i=\sum _{i=1}^n\sum _{v=1}^{m-1}\sum _{j=1}^{q-1}\ast P_i\ast \frac{1}{2\ast m+1}{P}_i\ast b_v\ast \sum _{k=0}^{j-1}{a}_i^j \\ =\frac{1}{2\ast m+1}\ast (\sum _{v=1}^{m-1}{b}_v)\ast \sum _{i=1}^n{P}_i^2(\sum _{j=1}^{q-1}\sum _{k=0}^{j-1}{a}_i^j) \end{gathered}$$
对$({\sum }_{j=1}^{q-1}{\sum }_{k=0}^{j-1}{a}_i^j)$可以设计矩阵快速幂$O(logn)$计算，或者经过数学推导可以证明：$$(\sum _{j=1}^{q-1}\sum _{k=0}^{j-1}{a}_i^j)=\frac{q-1-\frac{a_i(1-a_i^{q-1})}{1-a_i}}{1-a_i}$$
从而整个过程只需要枚举$i$，时间复杂度为$O(nlogn)$。
最后答案为$E\ast (\frac{n\ast (n+1)}{2}\ast (2\ast m+1){)}^q$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=4e5+5,mod=998244353;
int n,m,q;
ll inv[N];
ll qpow(ll a,ll n)
{
    ll ans=1;
    for(;n;n>>=1,a=a*a%mod)
        if(n&1) ans=ans*a%mod;
    return ans;
}
struct node
{
    int n,m;
    node(int n=0,int m=0):n(n),m(m){memset(a,0,sizeof(a));}
    ll a[4][4];
    node operator*(const node&o)const
    {
        node ans(n,o.m);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=o.m;j++)
            for(int k=1;k<=m;k++)
            ans.a[i][j]=(ans.a[i][j]+a[i][k]*o.a[k][j])%mod;
        return ans;
    }
};
node qpow(node a,ll n)
{
    node ans(a.n,a.n);
    for(int i=1;i<=a.n;i++) ans.a[i][i]=1;
    for(;n;n>>=1,a=a*a)
        if(n&1) ans=ans*a;
    return ans;
}
int main()
{
    inv[1]=1;for(int i=2;i<N;i++) inv[i]=(mod-mod/i)*inv[mod%i]%mod;
    cin>>n>>m>>q;
    ll ans=qpow(1ll*n*(n+1)/2%mod*(2*m%mod+1)%mod,q)%mod;
    ll sum=0;
    for(int v=1;v<m;v++)
        sum=(sum+(m-v)%mod*inv[2*m+1])%mod;
    ans=ans*sum%mod*inv[2*m+1]%mod;
    ll E=0;
    for(int i=1;i<=n;i++)
    {
        ll p=1ll*i*(n-i+1)%mod*inv[n]%mod*inv[n+1]%mod*2%mod;
        ll a=1-m*inv[2*m+1]%mod*p%mod;
        a=(a%mod+mod)%mod;
        //node A(3,3);
        //A.a[1][1]=A.a[2][1]=A.a[2][2]=A.a[3][2]=1;
        //A.a[3][3]=a;A=qpow(A,q-1);
        //sum=(A.a[2][1]+a*A.a[3][1]%mod)%mod;
        ll s=qpow((1+mod-a)%mod,mod-2);
        sum=q-1-a*(1-qpow(a,q-1))%mod*s%mod;
        sum=(sum%mod+mod)%mod;
        sum=sum*s%mod;
        E=(E+sum*p%mod*p%mod)%mod;
    }
    ans=ans*E%mod;
    printf("%lld\n",ans);
}