The_10958_problem_brute.py

ops = '+-*/^'
def calc(nums):
    if(nums==[10958]):
        print('hooray!')
        raise TypeError("found 10958")  
    for space in range(1,len(nums)):
        for op in ops:
            if op == '+':
                calc(nums[:space-1] + [nums[space-1]+nums[space]] + nums[space+1:])
            if op == '-':
                calc(nums[:space-1] + [nums[space-1]-nums[space]] + nums[space+1:])
            if op == '*':
                calc(nums[:space-1] + [nums[space-1]*nums[space]] + nums[space+1:])
            if op == '/':
                if nums[space]!=0:
                    try:
                        calc(nums[:space-1] + [nums[space-1]/nums[space]] + nums[space+1:])
                    except:
                        a=1 #do nothing, dead end
            if op == '^':
                if(nums[space]<1000):
                    try:
                        calc(nums[:space-1] + [nums[space-1] ** nums[space]] + nums[space+1:])
                    except:
                        a=1 #do nothing, dead end
                    
def splitnum(depth, root, num):
    #reversing the place where the space is put means that those with fewwer spaces are processed first
    for p in reversed(range(1,len(num))):
        lnum = num[0:p]
        rnum = num[p:]
        #convert to array of numbers
        txtnum = root+' '+lnum+' '+rnum
        nums = [int(x) for x in txtnum.split()]
        print(depth, nums)
        calc(nums)
        if(depth<6):#ship those with too many spaces, because they take too long to check all ops on
            splitnum(depth+1,root+' '+lnum, rnum)
            
orig = '123456789'
splitnum(1,'',orig)

Exhaustive search (with some cases skipped) No solution ¯_(ツ)_/¯