方程组的解数

蒜头君在求解一个 nnn 元的高次方程:

k1x1p1+k2x2p2+…+knxnpn=0\displaystyle k_1x_1{p_1}+k_2x_2{p_2}+\ldots+k_nx_n^{p_n}=0k1​x1p1​​+k2​x2p2​​+…+kn​xnpn​​=0

其中:x1,x2,…,xnx_1,x_2,\ldots,x_nx1​,x2​,…,xn​ 是未知数,k1,k2,…,knk_1,k_2,\ldots,k_nk1​,k2​,…,kn​ 是系数,p1,p2,…,pnp_1,p_2,\ldots,p_np1​,p2​,…,pn​ 是指数。方程中所有数都一定是整数。

假设未知数 1≤xi≤M,i=1…n1\leq x_i\leq M,i=1 \ldots n1≤xi​≤M,i=1…n。你能帮蒜头君算出这个方程的整数解个数吗?
输入格式

第一行输入一个整数 n(1≤n≤4)n(1 \leq n \leq 4)n(1≤n≤4)。

第二行输入一个整数 M(1≤M≤150)M(1 \leq M \leq 150)M(1≤M≤150)。

第 333 行到第 n+2n+2n+2 行,每行输入两个整数,分别表示 ki(∣ki∣≤20)k_i(|k_i| \leq 20)ki​(∣ki​∣≤20) 和 pi(1≤pi≤4)p_i(1 \leq p_i \leq 4)pi​(1≤pi​≤4)。两个整数之间用一个空格隔开。
输出格式

输出一行,输出一个整数,表示方程的整数解的个数。

输出时每行末尾的多余空格,不影响答案正确性
样例输入

3
100
1 2
-1 2
1 2

样例输出

104

#include <stdio.h>
#include <stdlib.h>
#include<string.h>//使用memset函数所需的头文件
int n,m;
int a[5];   //系数
int b[5];   //指数
int res = 0;//统计解的个数
int mypow[153][5];
void MyPow()  //递归计算次方,暴力打表(没想到比pow快?)
{
    int i,j;
    memset(mypow,1,sizeof(mypow));//赋值比二重循环快,使用二重循环就通过不了
    for(i=1;i<=150;i++)
        for( j=1;j<=4;j++)
            mypow[i][j] = i*mypow[i][j-1];//powmy[i][j]表示i^j
}
void dfs(int k,int val)
{
    if(k>n) return ;   //void 返回值不能写零,直接return,返回到调用函数处
    if(k==n)
    {
        if(val==0)
            res++;
        else
            return ;
    }
    int i;
	for(i=1;i<=m;i++)
    {
        dfs(k+1,val+a[k]*mypow[i][b[k]]);
    }
}

int main()
{
    scanf("%d%d",&n,&m);
    MyPow();
    int i;
    for( i=0;i<n;i++)
    {
       scanf("%d%d",&a[i],&b[i]);
    }
    dfs(0,0);
     printf("%d",res);
     return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值