poj 1840(移动几项到等式另外一边减少循环深度)

Eqs
Time Limit: 5000MS Memory Limit: 65536K
Total Submissions: 11220 Accepted: 5480

Description

Consider equations having the following form:
a1x1 3+ a2x2 3+ a3x3 3+ a4x4 3+ a5x5 3=0
The coefficients are given integers from the interval [-50,50].
It is consider a solution a system (x1, x2, x3, x4, x5) that verifies the equation, xi∈[-50,50], xi != 0, any i∈{1,2,3,4,5}.

Determine how many solutions satisfy the given equation.

Input

The only line of input contains the 5 coefficients a1, a2, a3, a4, a5, separated by blanks.

Output

The output will contain on the first line the number of the solutions for the given equation.

Sample Input

37 29 41 43 47

Sample Output

654

Source



上来写了一个四层for循环,结果直接tle了,优化了几发还是tle。

将前三项移动到等式右面,就可以降低到O(n^2 + n^3)的复杂度:

-(a1x13+ a2x23+ a3x33 )=  a4x43+ a5x53

把等式左面可能的取值放到一个数组中做hash,然后遍历右面的取值即可。

提交记录:

1,Wrong Answer 没有看到题目中的xi!=0 这个条件。

2.TLE 四层循环直接挂掉。。

3.Accepted 修改为三等循环+两层循环


代码:

#include 
   
   
    
    
#include 
    
    
     
     
#include 
     
     
      
      
#include 
      
      
       
       
#include 
       
       
         #include 
         #include 
         
           #include 
          
            #include 
           
             #include 
            
              #define oo 0x3f3f3f3f #define MAX_N 30010 #define MAX_M 510 using namespace std; int a1, a2, a3, a4, a5; int x1, x2, x3, x4, x5; int data[50*50*50*10] = {0}; int data1[60] = {0}; short int hash[50*50*50*50*4 + 100] = {0}; map 
             
               m; int main() { int i; for (i = 0; i <= 50; i++) { data[i*i*i] = i; data1[i] = i*i*i; } vector 
              
                ::iterator it; int mod = 50*50*50*50*2; while (EOF != scanf("%d%d%d%d%d", &a1, &a2, &a3, &a4, &a5)) { int result = 0, tmp; memset(hash, 0, sizeof(hash)); for (x1 = -50; x1 <= 50; x1++) { int x13; if (x1 < 0) x13 = - data1[-x1] * a1; else x13 = data1[x1] * a1; if (x1 != 0) for (x2 = -50; x2 <= 50; x2++) { int x23; if (x2 < 0) x23 = - data1[-x2] * a2; else x23 = data1[x2] * a2; if (x2 != 0) for (x3 = -50; x3 <= 50; x3++) { int x33; if (x3 < 0) x33 = - data1[-x3] * a3; else x33 = data1[x3] * a3; if (x3 != 0) { tmp = x13 + x23 + x33; tmp = -tmp; if (tmp < mod && tmp+mod >= 0) hash[tmp+mod]++; //if (m.find(tmp) != m.end()) m[tmp]++; //else m[tmp] = 1; } } } } for (x4 = -50; x4 <= 50; x4++) { if (x4 != 0) { int x43; if (x4 < 0) x43 = -data1[-x4] * a4; else x43 = data1[x4] * a4; for (x5 = -50; x5 <= 50; x5++) { if (!x5) continue;//排除了那种特殊情况 int x53; if (x5 < 0) x53 = -data1[-x5] * a5; else x53 = data1[x5] * a5; tmp = x43 + x53; result += hash[tmp+mod]; //if (m.find(tmp) != m.end()) result += m[tmp]; } } } cout << result << endl; } return 0; } 
               
              
             
            
           
          
       
      
      
     
     
    
    
   
   

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值