A Yet Another Remainder The 2022 ICPC Asia Regionals Online Contest (II)

最新推荐文章于 2024-10-05 11:02:16 发布

timidcatt

最新推荐文章于 2024-10-05 11:02:16 发布

阅读量205

点赞数

分类专栏：数论文章标签：算法 c++

本文链接：https://blog.csdn.net/ashbringer233/article/details/132782910

版权

数论专栏收录该内容

55 篇文章 1 订阅

订阅专栏

PTA | 程序设计类实验辅助教学平台

题目大意：有一个n位长的隐藏数x，从高位到低位依次标号为1到n，sum[i][j]表示从第i为开始每j位上的数的和，有q次询问，每次给出一个100以内除了5以外的质数p，问这个数%p等于多少

1<=t<=100；1<=n<=1e6；1<=i<=min(100,n)；1<=j<=i；1<=q<=10

思路：对于x的第i位，可以将其写成x[i]* $10^{n-i+1}$ ，x%p就相当于x[1]* $10^{n-1}$ %p+x[2]* $10^{n-2}$ %p+...+x[n]* $10^{0}$ %p。

如果 $10^{n-i+1}$ %p=x，那么我们可以将每一项x[i]* $10^{n-i+1}$ 写成x[i]*( $10^{n-i+1}$ -x)，这样的每一项%p都等于0，那么原式子不等于0的部分就是 $10^{n-1}$ %p*x[1]+ $10^{n-2}$ %p*x[2]+...+ $10^{0}$ %p*x[n]。

然后又发现 $10^{n-i+1}$ %p是有循环节的，且循环节的长度是p-1，那么我们可以将上式拆成p-1个式子例如第一个式子是 $10^{n-1}$ %p*(x[1]+x[p]+x[2*p-1]+...)，第二个式子就是 $10^{n-2}$ %p*(x[2]+x[p+1]+x[2*p]+...)。

每个式子里后面括号中的数都是题目中给出的即sum[p-1][第i位]，所以可以在O(qtp)的时间复杂度内求出

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int n;
ll sum[105][105];
void init()
{
   
}
void solve()
{
    cin>>n;
    for(int i=1;i<=min(100,n);i++)
    {
        for(int j=1;j<=i;j++)
        {
            cin>>sum[i][j];
        }
        
    }
    int q;
    cin>>q;
    for(int i=1;i<=q;i++)
    {
        ll mod;
        cin>>mod;
        vector<ll>m;
        ll x=1;
        for(int i=1;i<=mod-1;i++)
        {//求出一个循环节
            m.push_back(x%mod);
            x=x*10%mod;
        }
        ll nn=n;
        ll ans=0;
        ll y=min(mod-1,nn);//如果p比n大，直接取sum的最后一行即可
        for(int i=1;i<=mod-1;i++)
        {
            ll temp=nn%(mod-1);//第几位
            if(temp==0)
            {
                temp=mod-2;
            }
            else
            {
                temp--;
            }
            ans=(ans+m[temp]*sum[y][i]%mod)%mod;
            nn--;
        }
        cout<<ans<<endl;
    }
}
int main()
{
    cin.tie(0);
    cout.tie(0);
    ios::sync_with_stdio(false);
    int t;
    cin>>t;
    while(t--)
    {
        solve();
    }
    return 0;
}