计蒜客 方程的解数

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

\displaystyle k_1x_1^{p_1}+k_2x_2^{p_2}+\ldots+k_nx_n^{p_n}=0k1x1p1+k2x2p2++knxnpn=0

其中:x_1,x_2,\ldots,x_nx1,x2,,xn 是未知数,k_1,k_2,\ldots,k_nk1,k2,,kn 是系数,p_1,p_2,\ldots,p_np1,p2,,pn 是指数。方程中所有数都一定是整数。

假设未知数 1\leq x_i\leq M,i=1 \ldots n1xiM,i=1n。你能帮蒜头君算出这个方程的整数解个数吗?

输入格式

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

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

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

输出格式

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

样例输入
3
100
1 2
-1 2
1 2
样例输出
104

逐次尝试可能解 当存在解时  计数器+1     注意打表  pow多次运算会很慢:


ac:

#include <bits/stdc++.h>

using namespace std;

int n,m,a[10],b[10],value[10],sum=0,ans=0,mypow[160][5];

void fun()

{
	
	memset(mypow,1,sizeof(mypow));
	
	for(int i=1;i<=150;i++)
	
	{
		
		for(int j=1;j<=4;j++)
		
		mypow[i][j]=i*mypow[i][j-1];
	}
}

void dfs(int k)

{
	
	if(k==n+1)
	
	{
		
		sum=0;
		
		for(int i=1;i<=n;i++)
			
		sum+=a[i]*mypow[value[i]][b[i]];
		
		if(sum==0)
		
		ans++;
		
		return ;
	}
	
	for(int i=1;i<=m;i++)
	
	{
		
		value[k]=i;
		
		dfs(k+1);
	}
}

int main()

{
	
	cin>>n>>m;
	
    fun();
    
	for(int i=1;i<=n;i++)
	
	cin>>a[i]>>b[i];
	
	dfs(1);
	
	cout<<ans;
	
	return 0;
}

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值