写成数论风格


19

使用符号编写数学陈述:

  • There exists at least one non-negative integer(写为E存在量词)
  • All non-negative integers(写为A,通用量词)
  • + (加成)
  • * (乘法)
  • = (平等)
  • ><(比较运算符)
  • &(和),|(或),!(非)
  • ()(用于分组)
  • 变量名

等价于以下陈述

存在一个有理数a,使得π+ e * a是有理数。

(当然,π=3.1415 ...是等于圆周除以圆直径的数学常数,Ë=2.7182 ...欧拉数

您必须证明您的陈述确实等同于上述陈述。

显然,执行此操作的“最短”方法是证明该语​​句为真或为假,然后用琐碎的真或假声明进行回答,因为所有的真声明彼此都等效,所有的假陈述也彼此等同。

但是,给定语句的真值在数学上是一个未解决的问题:我们甚至不知道π+Ë是否不合理!因此,除非进行开创性的数学研究,否则面临的挑战是找到一个“简单”的等效语句,证明其等效性,并尽可能简短地对其进行描述。

计分

E A + * = > < & |并将!每个分数加1。(并且)不要添加任何分数。每个变量名称将分数加1。

例如,E x (A ba x+ba>x*(x+ba))得分13(E x A ba x + ba > x * x + ba

最低分获胜。


注意:

免责声明:本说明不是由OP编写的。

  • 这是不是一个挑战。答案不需要包含代码。
  • 这类似于但不像挑战,因为您需要编写一条语句并证明它与另一条语句等效。
  • 如果您可以证明上面的陈述是对/错,则可以提交平凡的Ax x=x陈述(例如,对于所有x,x = x )或平凡的陈述(例如,对于所有x,x> x Ax x>x)。
  • 允许您使用其他符号(类似于证明高尔夫球中的引理),但是分数将被计为与不使用它们相同。
    例如,如果定义a => b(!a) | b,则每次=>在证明中使用时,分数将增加2。
  • 由于常量未在允许的符号中列出,因此您不得使用它们。
    例如:语句1 > 0可以写成

    
    Forall zero: ( zero + zero = zero ) =>
    Forall one: ( Forall x: x * one = x ) =>
    one > zero
    

    得分为23。(请记住,=>每次使用费用为2)。

提示

  • 要使用自然常数,您可以这样做E0, 0+0=0 & E1, At 1*t=t &(因此您不需要=>哪个更广泛);对于大于1的数字,只需加上一些1

5
我喜欢这里的概念,但是声明确实很难写,无论评分如何,我都会对任何解决方案印象深刻。我建议使用更简单的方法,以便更多的人参与。
xnor

1
您需要一个与给定的等价的数学陈述。从什么意义上说它们应该等效?如果我是正确的,给定的陈述是错误的。因此,我很难理解它与其他陈述的等效性。例如,是否等价于存在一个有理数a,使得i + e * a是有理数(其中i是虚数单位)?
路易斯·门多

1
目前的说明只是说You are allowed to submit a trivially-true (e.g., for all x, x = x Ax x=x) or a trivially-false statement (e.g., for all x, x > x Ax x>x) if you can prove the statement above is true/false.。该声明现在尚未得到证明或反对,因此我真的不介意问题是否因为解决了问题而变得无聊
l4m2

1
所写的问题似乎在很大程度上掩盖了人们的观点,并且避免解释到底发生了什么,所以我在笔记中写了一点解释(挑战的非平凡性取决于给定声明当前未知的真值) 。
林恩

I'd be impressed by any solution no matter the score.得分只是针对那些可以解决这个问题的人。
l4m2

Answers:


27

671

E a (a+a>a*a & (E b (E c (E d (A e (A f (f<a | (E g (E h (E i ((A j ((!(j=(f+f+h)*(f+f+h)+h | j=(f+f+a+i)*(f+f+a+i)+i) | j+a<e & (E k ((A l (!(l>a & (E m k=l*m)) | (E m l=e*m))) & (E l (E m (m<k & g=(e*l+(j+a))*k+m)))))) & (A k (!(E l (l=(j+k)*(j+k)+k+a & l<e & (E m ((A n (!(n>a & (E o m=n*o)) | (E o n=e*o))) & (E n (E o (o<m & g=(e*n+l)*m+o))))))) | j<a+a & k=a | (E l (E m ((E n (n=(l+m)*(l+m)+m+a & n<e & (E o ((A p (!(p>a & (E q o=p*q)) | (E q p=e*q))) & (E p (E q (q<o & g=(e*p+n)*o+q))))))) & j=l+a+a & k=j*j*m))))))) & (E j (E k (E l ((E m (m=(k+l)*(k+l)+l & (E n (n=(f+m)*(f+m)+m+a & n<e & (E o ((A p (!(p>a & (E q o=p*q)) | (E q p=e*q))) & (E p (E q (q<o & j=(e*p+n)*o+q))))))))) & (A m (A n (A o (!(E p (p=(n+o)*(n+o)+o & (E q (q=(m+p)*(m+p)+p+a & q<e & (E r ((A s (!(s>a & (E t r=s*t)) | (E t s=e*t))) & (E s (E t (t<r & j=(e*s+q)*r+t))))))))) | m<a & n=a & o=f | (E p (E q (E r (!(E s (s=(q+r)*(q+r)+r & (E t (t=(p+s)*(p+s)+s+a & t<e & (E u ((A v (!(v>a & (E w u=v*w)) | (E w v=e*w))) & (E v (E w (w<u & j=(e*v+t)*u+w))))))))) | m=p+a & n=(f+a)*q & o=f*r)))))))) & (E m (m=b*(h*f)*l & (E n (n=b*(h*f+h)*l & (E o (o=c*(k*f)*i & (E p (p=c*(k*f+k)*i & (E q (q=d*i*l & (m+o<q & n+p>q | m<p+q & n>o+q | o<n+q & p>m+q))))))))))))))))))))))))))

怎么运行的

首先,乘以a和(π+ e·a)的公称分母来重写条件为:存在a,b,c∈ℕ(并非全为零),其中a·π+ b·e = c或a·π-b·e = c或−a·π+ b·e = c。处理标志问题需要三个案例。

然后,我们需要重写它以通过有理逼近来谈论π和e:对于所有有理逼近π₀<π<π₁和e₀<e <e₁,我们有a·π₀+ b·e₀<c <a·π₁+ b·e₁或a·π₀-b·e₁<c <a·π₁+ b·e₀或-a·π₁+ b·e₀<c <-a·π₀+ b·e₁。(请注意,我们现在免费获得“并非全零”条件。)

现在是困难的部分。我们如何获得这些有理逼近?我们想使用像

2/1·2/3·4/3·4/5⋯(2·k)/(2·k +1)<π/ 2 <2/1·2/3·4/3·4/5⋯ (2·k)/(2·k +1)·(2·k + 2)/(2·k +1),

((k + 1)/ k)k <e <((k + 1)/ k)k +1

但是没有明显的方法来编写这些产品的迭代定义。因此,我们建立了一些我在Quora帖子中首先介绍的机制。限定:

除(d,a):=∃b,a = d·b,

powerOfPrime(a,p):=∀b,((b> 1并除以(b,a))⇒除以(p,b)),

如果a = 1或p = 1或p是素数且a是它的幂,则满足。然后

isDigit(a,s,p):= a <p和∃b,(powerOfPrime(b,p)和∃qr,(r <b和s =(p·q + a)·b + r))

当a = 0或a是基数p的数字时满足。这使我们可以使用一些基数p的数字表示任何有限集。现在,我们可以通过大致书写来转换迭代计算,其中存在一组中间状态,以使最终状态在该组中,并且该组中的每个状态都为初始状态,或者与该状态中的某个其他状态相距一个步骤。组。

详细信息在下面的代码中。

Haskell中生成代码

{-# LANGUAGE ImplicitParams, TypeFamilies, Rank2Types #-}

-- Define an embedded domain-specific language for propositions.
infixr 2 :|

infixr 3 :&

infix 4 :=

infix 4 :>

infix 4 :<

infixl 6 :+

infixl 7 :*

data Nat v
  = Var v
  | Nat v :+ Nat v
  | Nat v :* Nat v

instance Num (Nat v) where
  (+) = (:+)
  (*) = (:*)
  abs = id
  signum = error "signum Nat"
  fromInteger = error "fromInteger Nat"
  negate = error "negate Nat"

data Prop v
  = Ex (v -> Prop v)
  | Al (v -> Prop v)
  | Nat v := Nat v
  | Nat v :> Nat v
  | Nat v :< Nat v
  | Prop v :& Prop v
  | Prop v :| Prop v
  | Not (Prop v)

-- Display propositions in the given format.
allVars :: [String]
allVars = do
  s <- "" : allVars
  c <- ['a' .. 'z']
  pure (s ++ [c])

showNat :: Int -> Nat String -> ShowS
showNat _ (Var v) = showString v
showNat prec (a :+ b) =
  showParen (prec > 6) $ showNat 6 a . showString "+" . showNat 7 b
showNat prec (a :* b) =
  showParen (prec > 7) $ showNat 7 a . showString "*" . showNat 8 b

showProp :: Int -> Prop String -> [String] -> ShowS
showProp prec (Ex p) (v:free) =
  showParen (prec > 1) $ showString ("E " ++ v ++ " ") . showProp 4 (p v) free
showProp prec (Al p) (v:free) =
  showParen (prec > 1) $ showString ("A " ++ v ++ " ") . showProp 4 (p v) free
showProp prec (a := b) _ =
  showParen (prec > 4) $ showNat 5 a . showString "=" . showNat 5 b
showProp prec (a :> b) _ =
  showParen (prec > 4) $ showNat 5 a . showString ">" . showNat 5 b
showProp prec (a :< b) _ =
  showParen (prec > 4) $ showNat 5 a . showString "<" . showNat 5 b
showProp prec (p :& q) free =
  showParen (prec > 3) $
  showProp 4 p free . showString " & " . showProp 3 q free
showProp prec (p :| q) free =
  showParen (prec > 2) $
  showProp 3 p free . showString " | " . showProp 2 q free
showProp _ (Not p) free = showString "!" . showProp 9 p free

-- Compute the score.
scoreNat :: Nat v -> Int
scoreNat (Var _) = 1
scoreNat (a :+ b) = scoreNat a + 1 + scoreNat b
scoreNat (a :* b) = scoreNat a + 1 + scoreNat b

scoreProp :: Prop () -> Int
scoreProp (Ex p) = 2 + scoreProp (p ())
scoreProp (Al p) = 2 + scoreProp (p ())
scoreProp (p := q) = scoreNat p + 1 + scoreNat q
scoreProp (p :> q) = scoreNat p + 1 + scoreNat q
scoreProp (p :< q) = scoreNat p + 1 + scoreNat q
scoreProp (p :& q) = scoreProp p + 1 + scoreProp q
scoreProp (p :| q) = scoreProp p + 1 + scoreProp q
scoreProp (Not p) = 1 + scoreProp p

-- Convenience wrappers for n-ary exists and forall.
class OpenProp p where
  type OpenPropV p
  ex, al :: p -> Prop (OpenPropV p)

instance OpenProp (Prop v) where
  type OpenPropV (Prop v) = v
  ex = id
  al = id

instance (OpenProp p, a ~ Nat (OpenPropV p)) => OpenProp (a -> p) where
  type OpenPropV (a -> p) = OpenPropV p
  ex p = Ex (ex . p . Var)
  al p = Al (al . p . Var)

-- Utility for common subexpression elimination.
cse :: Int -> Nat v -> (Nat v -> Prop v) -> Prop v
cse uses x cont
  | (scoreNat x - 1) * (uses - 1) > 6 = ex (\x' -> x' := x :& cont x')
  | otherwise = cont x

-- p implies q.
infixl 1 ==>

p ==> q = Not p :| q

-- Define one as the unique n with n+n>n*n.
withOne ::
     ((?one :: Nat v) =>
        Prop v)
  -> Prop v
withOne p =
  ex
    (\one ->
       let ?one = one
       in one + one :> one * one :& p)

-- a is a multiple of d.
divides d a = ex (\b -> a := d * b)

-- a is a power of p (assuming p is prime).
powerOfPrime a p = al (\b -> b :> ?one :& divides b a ==> divides p b)

-- a is 0 or a digit of the base-p number s (assuming p is prime).
isDigit a s p =
  cse 2 a $ \a ->
    a :< p :&
    ex
      (\b -> powerOfPrime b p :& ex (\q r -> r :< b :& s := (p * q + a) * b + r))

-- An injection from ℕ² to ℕ, for representing tuples.
pair a b = (a + b) ^ 2 + b

-- πn₀/πd < π/4 < πn₁/πd, with both fractions approaching π/4 as k
-- increases:
-- πn₀ = 2²·4²·6²⋯(2·k)²·k
-- πn₁ = 2²·4²·6²⋯(2·k)²·(k + 1)
-- πd = 1²⋅3²·5²⋯(2·k + 1)²
πBound p k cont =
  ex
    (\s x πd ->
       al
         (\i ->
            (i := pair (k + k) x :| i := pair (k + k + ?one) πd ==>
             isDigit (i + ?one) s p) :&
            al
              (\a ->
                 isDigit (pair i a + ?one) s p ==>
                 ((i :< ?one + ?one :& a := ?one) :|
                  ex
                    (\i' a' ->
                       isDigit (pair i' a' + ?one) s p :&
                       i := i' + ?one + ?one :& a := i ^ 2 * a')))) :&
       let πn = x * k
           πn = πn + x
       in cont πn πn πd)

-- en₀/ed < e < en₁/ed, with both fractions approaching e as k
-- increases:
-- en₀ = (k + 1)^k * k
-- en₁ = (k + 1)^(k + 1)
-- ed = k^(k + 1)
eBound p k cont =
  ex
    (\s x ed ->
       cse 3 (pair x ed) (\y -> isDigit (pair k y + ?one) s p) :&
       al
         (\i a b ->
            cse 3 (pair a b) (\y -> isDigit (pair i y + ?one) s p) ==>
            (i :< ?one :& a := ?one :& b := k) :|
            ex
              (\i' a' b' ->
                 cse 3 (pair a' b') (\y -> isDigit (pair i' y + ?one) s p) ==>
                 i := i' + ?one :& a := (k + ?one) * a' :& b := k * b')) :&
       let en = x * k
           en = en + x
       in cont en en ed)

-- There exist a, b, c ∈ ℕ (not all zero) with a·π/4 + b·e = c or
-- a·π/4 = b·e + c or b·e = a·π/4 + c.
prop :: Prop v
prop =
  withOne $
  ex
    (\a b c ->
       al
         (\p k ->
            k :< ?one :|
            Bound p k $ n πn πd ->
               eBound p k $ \en en ed ->
                 cse 3 (a * πn * ed) $ \x ->
                   cse 3 (a * πn * ed) $ \x ->
                     cse 3 (b * en * πd) $ \y ->
                       cse 3 (b * en * πd) $ \y ->
                         cse 6 (c * πd * ed) $ \z ->
                           (x + y :< z :& x + y :> z) :|
                           (x :< y + z :& x :> y + z) :|
                           (y :< x + z :& y :> x + z))))

main :: IO ()
main = do
  print (scoreProp prop)
  putStrLn (showProp 0 prop allVars "")

在线尝试!


“如果a = 1或p是素数且a是它的幂,则可以满足”-您还可以拥有p =1。尽管p> 1隐含表示isDigit,但您使用它的唯一位置。
与Orjan约翰森

@ØrjanJohansen谢谢,我修正了那张纸条。(只要有某种方法可以表示每个有限集,实际上在哪个集合中出现powerOfPrimeisDigit最终表示都没关系。)
Anders Kaseorg

2
a我认为,如果得分为7或更高,那么就值得添加ex (\a' -> a' := a :& ... )包装器isDigit
与Orjan约翰森

@ØrjanJohansen当然,可以节省68。谢谢!
安德斯·卡塞格

我相信您需要使用k>0,因为eBound在这种k==0情况下给出了一个零分母(和一个零分子),所以所有替代方法都失败了。
与Orjan约翰森

3

270

E1                                                                              { Exist 1, defined when Any k introduced }
Ec1 Ec2 Ec3 Ec4 Ec5 Ak k*1=k & c3>1 & ( En0 An n<n0 |                           { for large enough n, |(c1-c4)e+c3(4-pi)/8+(c2-c5)|<1/k }
Ex Ep Ew Emult At (Eb ((b>1 & Eh b*h=t) &! Eh h*p=b)) |                         { x read in base-p, then each digit in base-w. t as a digit }
Ee1 Ee2 Ehigher Elower e2<p & lower<t & ((higher*p+e1)*p+e2)*t+lower=x &        { last digit e1, this digit e2 }
    { Can infer that e2=w+1 | e1<=e2 & u1<=u2 & i1<=i2 & s1<=s2 & t1<=t2, so some conditions omitted }
Ei1 Es1 Et1 Eu1 (((u1*w)+i1)*w+t1)*w+s1=e1 &                                    { (u,i,t,s) }
Ei2 Es2 Et2 Eu2 i2<w & s2<w & t2<w & (((u2*w)+i2)*w+t2)*w+s2=e2 &               { e2=1+w is initial state u=i=0, s=t=1 }
(e2=w+1 | e1=e2 | i2=i1+1+1 & s2=s1*(n+1) & t2=t1*n &                           { i=2n, s=(n+1)^n, mult=t=n^n, s/mult=e }
Eg1 Eg2 g1+1=(i2+i2)*(i2+i2) & g1*u1+mult=g1*u2+g2 & g2<g1) &                   { u/mult=sum[j=4,8,...,4n]1/(j*j-1)=(4-pi)/8. mult=g1*(u2-u1)+g2 }
(t>1 | i2=n+n & t2=mult & Ediff Ediff2                                          { check at an end t=1 }
c1*s2+c2*mult+c3*u2+diff=c4*s2+c5*mult+diff2 & k*(diff+diff2)<mult))            { |diff-diff2|<=diff+diff2<mult/k, so ...<1/k }

a|b&ca|(b&c)因为我想删除这些括号使它看起来更好,反正他们是免费的。

使用JavaScript "(expr)".replace(/\{.*?\}/g,'').match(/[a-z0-9]+|[^a-z0-9\s\(\)]/g)对令牌进行计数。


为什么要服用mult = t?另外,由于x只能包含有限的多个数字,因此您需要e1 = e2 = 0考虑足够大的t。此外,对于类似的歧义结构,您还需要更多括号或其他歧义_ & _ | _
安德斯·卡塞格

@AndersKaseorg我将每个项目相乘multmult=t2最后看不到任何问题。e1=e2=0应该是固定的,但不确定,因此我目前不更改接受度。
l4m2

如果a & b | c是,(a & b) | c那么您t*1=t肯定在错误的地方。同样,您也没有排除平凡的解决方案c1 = c4 & c2 = c5 & c3 = 0 & diff = diff2
安德斯·卡塞格

@AndersKaseorg我的diff≠diff2工作理由吗?
l4m2

无论如何!(c2=c5),我们可以使用,因为我们已经知道这e是不合理的,所以即使这个无效的分数也不会增加
l4m2
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.