|
| |
|
|
A081229
|
|
a(n) is the sum of the common digits of n and the n-th prime in base 10, or -1 if there are no common digits.
|
|
2
| |
|
|
-1, -1, -1, -1, -1, -1, 7, -1, -1, -1, 1, -1, 1, 4, -1, -1, -1, 1, -1, -1, -1, -1, 3, -1, -1, -1, -1, -1, 9, 3, 1, 3, 3, 3, -1, -1, 7, 3, -1, -1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 2, -1, 5, 5, 6, -1, -1, -1, -1, -1, 2, 3, -1, -1, -1, -1, -1, -1, -1, -1, -1, 10, 7, 7, -1, -1, 7, -1, 0, 1, 2, 3, 4, -1, -1, -1, -1, -1, -1, -1, 9, -1, 13, 9, -1
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,7
|
|
|
COMMENTS
| Repeated digits are to be ignored; e.g. 2001 and 4050 have only one digit in common, 0.
|
|
|
EXAMPLE
| a(73)=10 since the 73rd prime is 367 (3+7=10).
|
|
|
MAPLE
| digcomp := proc(A, B) local a, b, crit, f; description "returns the digits that 'A' and 'B' have in common; if a third argument is given and it is 'sum' then the sum of the common digits is returned, else if it is 'num' then the number of common digits is returned."; if nargs>2 then crit := args[3] else crit := NULL fi; a := convert(`if`(A=0, [0], (convert(A, base, 10))), set); b := convert(`if`(B=0, [0], (convert(B, base, 10))), set); f := a intersect b; if crit=sum then return `if`(nops(f)>0, `+`(op(f)), -1) elif crit=num then return nops(f) else return f; fi; end proc; F := ([seq(digcomp(i, ithprime(i), sum), i=1..75)]);
|
|
|
PROG
| (PARI) digitset(n) = local(v, d); v=[]; while(n>0, d=divrem(n, 10); n=d[1]; v=concat(d[2], v)); Set(v)
{for(n=1, 94, s=setintersect(digitset(n), digitset(prime(n))); v=eval(s); print1(if(v==[], -1, sum(j=1, #v, v[j])), ", "))} [Klaus Brockhaus, Dec 19 2006]
|
|
|
CROSSREFS
| Cf. A081227, A081228.
Sequence in context: A070672 A055061 A074465 * A109010 A117825 A010143
Adjacent sequences: A081226 A081227 A081228 * A081230 A081231 A081232
|
|
|
KEYWORD
| base,sign
|
|
|
AUTHOR
| Francois Jooste (pin(AT)myway.com), Mar 11 2003
|
|
|
EXTENSIONS
| Entries checked by Klaus Brockhaus, Dec 19 2006
|
| |
|
|
