A Fun Perfect Square Checker Using Integer Arithmetic Only... ;o)

 
Thread Tools Search this Thread
Operating Systems OS X (Apple) A Fun Perfect Square Checker Using Integer Arithmetic Only... ;o)
# 22  
Old 08-24-2015
Quote:
Originally Posted by wisecracker
Hi Don...

Phew, that took some reading...

As 0^2 is 0, I never even considered it to be a perfect square so I stand corrected.
My error; however my code can very easily be changed to cater for it.

As you say it is slow, but it was coded from an observation not from any mathematical set of algorithms, that is, every _line_no_ is the perfect square of the odd number series...
Yes, it was a long post. I hope a few people found the narrative informative.

My code rejected 0 since your code didn't process it correctly. If you'd like to modify my script to accept 0, here is a diff showing the changes needed:
Code:
48,49c48,49
< 	    [[ num -lt 1 ]]
< 	then	printf "Out of range: 0 < x <= $LONG_MAX: x=$num\n"
---
> 	    [[ num -lt 0 ]]
> 	then	printf "Out of range: 0 <= x <= $LONG_MAX: x=$num\n"
58c58
< 		low=$((10 ** ((${#num} - 1) / 2)))
---
> 		low=$((10 ** ((${#num} - 1) / 2) - 1))

Your code used the observation that the sum of all of the odd numbers from 0 to x is a perfect square. My code used the observation that the square root of any positive x digit number has a square root between 10^((x-1)/2) and 10^(((x-1)/2)+1)-1 inclusive (except for the number 0 [since 10^0 is 1; not 0]).

Your code starts with the lowest perfect square and calculates successive perfect squares until it has a perfect square that is greater than or equal to the number being checked. My code performs a binary search of the possible square roots that could match the target number until it determines that there is a square root in range that is the square root of the target or determines that one does not exist.

For a 63-bit integer value, your code may need to loop through millions of possible perfect squares. A binary search like my code uses needs to loop through no more than 32 possible square roots. My code needs less than that since it restricts the high and low possible square roots based on the number of digits in the target number.
# 23  
Old 08-24-2015
OK, here is another take at finding the square root: it does not use any operation other than subtraction/addition and it is no approximation method.

The math is simple: we know that

Code:
           2    2          2
  ( a + b )  = a  + 2ab + b

or, in the special case of b = 1:

Code:
           2    2          
  ( n + 1 )  = n  + 2n + 1

which we can more conveniently write as:

Code:
           2    2
  ( n + 1 )  = n  + n + ( n + 1 ))

This means that
Code:
   n
  ---             2
  \   2i + 1   = n
  /__ 
  i=0

In less mathematical terms: we have two counters, x and y. We start with x=0, y=1 and we add them both to a sum, incrementing the counters at each pass. Here is a little table:

Code:
  Pass  x y   sum
  1     0 1     1
  2     1 2     4
  3     2 3     9
  4     3 4    16

As you can see the "y"-column always holds the root to the square in the "sum"-column.


To find the root i simply reversed that and subtracted (instead of added) from the value originally given down to 0. The code is very simple:

Code:
#! /bin/bash

typeset -i iInput=$1
typeset -i iCnt1=0
typeset -i iCnt2=1

while [ $iInput -gt 0 ] ; do
     (( iInput = iInput - iCnt1 ))
     if [ $iInput -lt $iCnt2 ] ; then
          echo "Not a perfect square."
          exit 0
     fi
     (( iInput = iInput - iCnt2 ))
     (( iCnt1++ ))
     (( iCnt2++ ))
done

echo Solution: $iCnt1

exit 0

bakunin
# 24  
Old 08-24-2015
Hi,
For fun in bash:
Code:
#!/bin/bash

XX=$1
YY=$((${#XX}&1 ? 1 : 2))
ZZ=${XX:0:$YY}
AA=$(((${#XX}-YY)/2))
VV=${XX:$YY:$((AA*2))}
BB=0
DD=1
EE=1
while (($EE))
do
	ZZ=$((ZZ-DD))
	BB=$((ZZ>=0 ? BB+1 : BB))
	DD=$((ZZ>=0 ? DD+2 : DD-2))
	EE=$((ZZ>=0 ? 1 : 0))
	ZZ=$((ZZ>=0 ? ZZ : ZZ+DD+2))
done
rac=$rac$BB
II=0
while (($AA))
do
	BB=0
	EE=1
	DD=$((((DD+1)*10)+1))
	ZZ=$((ZZ*100))
	UU=${VV:$II:1}
	UU=$((UU*10))
	UU=$((UU+${VV:II+1:1}))
	ZZ=$((ZZ+UU))
	while (($EE))
	do
		ZZ=$((ZZ-DD))
		BB=$((ZZ>=0 ? BB+1 : BB))
		DD=$((ZZ>=0 ? DD+2 : DD-2))
		EE=$((ZZ>=0 ? 1 : 0))
		ZZ=$((ZZ>=0 ? ZZ : ZZ+DD+2))
	done
	rac=$rac$BB	
	II=$((II+2))
	AA=$((AA-1))
done
(( ! ZZ )) && echo $XX is perfect square of $rac || echo $XX is not a perfect square

Examples:
Code:
$ ./p_square.sh 391362529049165025
391362529049165025 is perfect square of 625589745
$ ./p_square.sh 98743978937000250000
98743978937000250000 is perfect square of 9937000500
$ ./p_square.sh 98743978937000250011
98743978937000250011 is not a perfect square

Login or Register to Ask a Question

Previous Thread | Next Thread

8 More Discussions You Might Find Interesting

1. Shell Programming and Scripting

Not able to find the perfect code...Geting confused in between

I have to find last delimiter in each line of a file and store the value after the last '/' in a variable in ksh script...Pls Pls help me:(The file is as shown below: /opt/apps/cobqa/apps/abadv/bind/advc0007.bnd /opt/apps/cobqa/apps/abbrio/bind/naac6115.bnd... (5 Replies)
Discussion started by: bhavanabahety
5 Replies

2. Shell Programming and Scripting

egrep line with perfect mach

Hi Input File A L006 AL01 0 (OCK) L006 A006 0 (OCK) L011 AR11 1 (NLOCK) Input File B L006 AL01 0 (OCK) L006 A006 0 (OCK) Need Egrep Command for perfect Match Thanks (4 Replies)
Discussion started by: asavaliya
4 Replies

3. UNIX for Dummies Questions & Answers

Can you perfect my sed ?

I want to print only the lines that meet the criteria : "worde:" and "wordo;" I got this far: sed -n '/\(*\)\1e:\1o;/p;' But it doesn't quite work. Can someone please perfect it and tell me exactly how its a fixed version/what was wrong with mine? Thanks heaps, (1 Reply)
Discussion started by: maximus73
1 Replies

4. Shell Programming and Scripting

how to compare string integer with an integer?

hi, how to I do this? i="4.000" if ; then echo "smaller" fi how do I convert the "4.000" to 4? Thanks! (4 Replies)
Discussion started by: h0ujun
4 Replies

5. Shell Programming and Scripting

Delete text between square brackets and also delete those square brackets using sed or awk

Hi All, I have a text file which looks like this: computer programming systems engineering I want to get rid of these square brackets and also the text that is inside these brackets. So that my final text file looks like this: computer programming systems engineering I am using... (3 Replies)
Discussion started by: shoaibjameel123
3 Replies

6. AIX

I want the perfect user-interface

I've got an aix-box somewhere on the network and a PC on my desk. Nothing fancy so far. The PC is made dual-boot: - windowsXP with putty & winSCP or - slackware 13 with xfce4 installed. The aix-box runs DB2 v8.2 and I've installed db2top to monitor the database. db2top is a character... (0 Replies)
Discussion started by: dr_te_z
0 Replies

7. UNIX for Dummies Questions & Answers

A perfect number shell program

Here's my work of testing whether a number input is perfect or not.. echo Enter a number read no i=1 ans=0 while do if then ans='expr $ans + $i' fi i='expr $i + 1' done if then echo $no is perfect else echo $no is NOT perfect fi (12 Replies)
Discussion started by: Cyansnow
12 Replies

8. Shell Programming and Scripting

extraction of perfect text from file.

Hi All, I have a file of the following format. <?xml version='1.0' encoding='utf-8'?> <tomcat-users> <role rolename="tomcat"/> <role rolename="role1"/> <role rolename="manager"/> <role rolename="admin"/> <user username="tomcat" password="tomcat" roles="tomcat"/> <user... (5 Replies)
Discussion started by: nua7
5 Replies
Login or Register to Ask a Question