方程的解

最新推荐文章于 2023-06-17 16:37:10 发布

FrEddy_Ffy

最新推荐文章于 2023-06-17 16:37:10 发布

阅读量343

点赞数

分类专栏：题解

本文链接：https://blog.csdn.net/fengfeiyangxx/article/details/68073639

版权

题解专栏收录该内容

5 篇文章 0 订阅

订阅专栏

题目描述

　佳佳碰到了一个难题，请你来帮忙解决。　　对于不定方程a1+a2+……+ak-1+ak=g(x)，其中k>=2且k∈N*，x是正整数，g(x)=x^x mod 1000（即xx除以1000的余数），x，k是给定的数。我们要求的是这个不定方程的正整数解组数。　　举例来说，当k=3,x=2时，分别为(a1,a2,a3)=(2,1,1),(1,2,1),(1,1,2).

输入描述

输入文件有且只有一行，为用空格隔开的两个正整数，依次为k，x。

输出描述

　输出文件有且只有一行，为方程的正整数解组数。

范围

对于40%的数据，ans<=10^16; 对于100%的数据，k<=100,x<=2^31-1,k<=g(x)。

输入

3 2

输出

3

题意

很明显，不必多讲。

分析

设 x’=（（X^X）mod 1000）。相当于求（（X^X）mod 1000）个球，放进（k-1）个抽屉里的方案数，答案是C（x’-1，k-1）。
X^X可以用快速幂求，“一边幂一边mod”，输出较大，要用到高精度。

标程

#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
using namespace std;
long long get(long long x,long long z)
{
    if(x==1)
        return z;
    if(x==0)
        return 1;
    long long w=get(x/2,z);
    if(x%2==0)
        return (w*w)%1000;
    else
        return (w*w*z)%1000;
}
int ans[600];
void cheng(int x)
{
    for(int i=1;i<=ans[0];i++)
        ans[i]=ans[i]*x;
    for(int i=1;i<ans[0];i++)
    {
        ans[i+1]+=ans[i]/10;
        ans[i]=ans[i]%10;
    }
    while(ans[ans[0]]>=10)
    {
        ans[ans[0]+1]+=ans[ans[0]]/10;
        ans[ans[0]]=ans[ans[0]]%10;
        ans[0]++;
    }
}
void chu(int x)
{
    for(int i=ans[0];i>=1;i--)
    {
        ans[i-1]+=ans[i]%x*10;
        ans[i]/=x;
    }
    while(ans[ans[0]]==0&&ans[0]>1)
        ans[0]--;
}
void C(int x,int y)
{
    for(int i=1;i<=y;i++)
        cheng(x-i+1);
    for(int i=1;i<=y;i++)
        chu(i);
    for(int i=ans[0];i>=1;i--)
        printf("%d",ans[i]);
}
int main()
{
    long long k,x;
    scanf("%lld%lld",&k,&x);
    long long z=get(x,x%1000);
    ans[1]=1;
    ans[0]=1;
    C(z-1,k-1);
    return 0;
}