Dotaz: nahrazeni prvnich x vyrazu jinym vyrazem

[Tomas Hrobarik]

Zdravim vsechny, potreboval bych nahradit v jednom soubore prvnich x vyrazu jinym vyrazem. zkousel jsem

sed 's/vyraz1\{x\}/vyraz2/g'

ale to nefunguje. Nevite mi nekdo prosim poradit. diky

Odpovědi

[Michal Marek (twofish)]

\{x\} má vyšší prioritu než konkatenace, viz info sed, man grep atd.

sed 's/\(vyraz1\)\{x\}/vyraz2/g'

[Tomas Hrobarik]

sed 's/\(vyraz1\)\{x\}/vyraz2/g' nepomohlo, resp. nic nezamenilo. V cem je problem?

[unchallenger]

Zkus napsat, co opravdu nahrazuješ, ne vyraz1 a x, protože bota bude někde jinde; co napsal twofish, to funguje.

[Tomas Hrobarik]

napsal jsem sed 's/\(Na\)\{45\}/ K/g' ala10_nacl.pdb > pokus

[unchallenger]

Hm, možná jsem teprve pochopil, co chceš. Sed umí opak, tj. nahradit výskyty až od nějakého, a nahradit konkrétní výskyt, ale nahradit prvních n AFAIK snadno nejde. Lze to udělat např.

sed '-e s/výraz/náhrada/'{,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,}

kde těch čárek je n-1. Ovšem pouze v případě, kdy náhrady nemohou vytvořit nová místa, kde se výraz matchne.

[Tomas Hrobarik]

Ja netvrdim aby to byl striktne sed. Ale verim, ze existuje nejaky programek, ktery dokaze to co prave potrebuji a to je nahradit v souboru prvnich 45 vyrazu Na za vyraz K. btw.

sed '-e s/výraz/náhrada/'{,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,}

mi taky nepomohlo bez ohlkadu na pocet carek

[Tomas Hrobarik]

pro lepsi nazornost prikladam ukazku souboru ktery chci zmenit

ATOM  15040  Na   Na  5012      36.160  49.470  53.120  1.00  0.00
ATOM  15041  Na   Na  5013      24.520  52.110  42.150  1.00  0.00
ATOM  15042  Na   Na  5014      51.240   1.070  22.880  1.00  0.00
ATOM  15043  Na   Na  5015      51.240   1.070  41.500  1.00  0.00
ATOM  15044  Na   Na  5016       0.710  52.110  11.040  1.00  0.00
ATOM  15045  Na   Na  5017      37.510  53.200  19.790  1.00  0.00
ATOM  15046  Na   Na  5018       1.610  47.900  27.350  1.00  0.00
ATOM  15047  Na   Na  5019      43.140  52.110  23.530  1.00  0.00
ATOM  15048  Na   Na  5020      37.510  53.200  38.410  1.00  0.00
ATOM  15049  Na   Na  5021       1.610  47.900  45.970  1.00  0.00
ATOM  15050  Na   Na  5022      43.140  52.110  42.150  1.00  0.00
ATOM  15051  Cl   Cl  5023      10.130   1.050  17.700  1.00  0.00
ATOM  15052  Cl   Cl  5024       9.030   0.860   1.330  1.00  0.00
ATOM  15053  Cl   Cl  5025       3.900  17.410  15.600  1.00  0.00
ATOM  15054  Cl   Cl  5026       8.590  13.740   0.160  1.00  0.00
ATOM  15055  Cl   Cl  5027      10.010  53.925  11.540  1.00  0.00
ATOM  15056  Cl   Cl  5028      10.130   1.050  36.320  1.00  0.00

[Michal Kubeček]

Od oka by na tohle asi bylo jednodušší použít awk…

[Tomas Hrobarik]

a mohl bys mi poradit nejaky skript. Ja to zas az tak neovladam.

[unchallenger]

{
  if (n < 45)
    n += gsub(/a/, "ble")
  print
}

Teda nic moc, protože se na řádku, kde to překročí 45 výskytů, nahradí všechny. Lze se tomu vyhnout, ale sehlaly mi už tři teorie, co vlastně chceš, takže se nebudu pouštět do větších akcí...

[Michal Marek (twofish)]

Něco na způsob

awk '{while (i < 45 && sub("Na", " K")) i++; print; }'

sub() vrací 1 pokud se provedla substituce. Má to zase ten problém, že náhrada nesmí jít nahradit podruhé.

[unchallenger]

Pakliže víš, co chce, já vím, jak to řešit správně ;-)

[Michal Marek (twofish)]

_Algoritmus_, jak to řešit správně, vím, nevím jak to _úsporně_ zapsat v awku... :-|

[unchallenger]

Hm, OK, nejlepší, na co jsem přišel, používalo gsub() + půlení intervalu na substituci posledních výskytů, aby jich bylo správně 45.

No, než to Pichi napíše v dc, tak Python:

#/usr/bin/python
import re, sys

r = re.compile(r'vyraz')
repl = 'replacement'
m = 13

n = 0
for line in sys.stdin:
    if n < m:
        i = len(r.findall(line))
        line = r.sub(repl, line, max(i + n - m, 0))
        n += i
    sys.stdout.write(line)

[unchallenger]

s/13/45/

[Hynek (Pichi) Vychodil]

dc na řetězce není nic moc. Uvažoval bych o perlu

perl -pe 'next if $a>45; while ($a++ <45 and s/vyraz/replacement/){};'

[Hynek (Pichi) Vychodil]

Nebo prakticky totéž

perl -pe 'next if $a>=45; while ($a++ <45 and s/vyraz/replacement/){};

a nebo

perl -pe 'next unless $a<45; while ($a++ <45 and s/vyraz/replacement/){};

[unchallenger]

To by tady fungovalo. Ale když jsem chtěl, aby se to chovalo korektně, i když replacement může vytvořit nová místa, kde matchne vyraz, skončil jsem u Pythonu.

[Hynek (Pichi) Vychodil]

Jojo. Máš pravdu. Tohle se v perlu moc elegantně nedělá. Nicméně v zadání takový požadavek není. Krom toho jsem si to přes oběd nechal trohu projít lavou a mám tam dost blbou chybu a asi už trpím stařeckou demencí (ve 30) a nejsem schopen opustit myšlenku - next. Takže snad na čtvrtý pokus

perl -pe '$a++ while $a<45 and s/vyraz/replace/'

[Hynek (Pichi) Vychodil]

Elegance se poněkud vytrácí, ale jde to.

perl -pe '$l=0; while($a<45 and substr($_,$l) =~ s/vyraz/replace/){$a++; $l=length($_)-length($'\'')}'

[Tomas Hrobarik]

Diky. Moc jste mi pomohli. Btw. nevite o nejaky dobry strance na webu s tutorialama jak psat skripty?

[tomas]

Hodně skriptů pro práci s PDB soubory je napsáno v Perlu. I původní PDB databáze z BNL byla napsána z 90% v Perlu. Při zpracování PDB souborů se lítají regulárními výrazy jen za druhým a Perl je v tomto prvotřídní.

Co se týče otázky tutoriálu Perlu bych zkusil třeba tady:

http://docs.linux.cz/programming/interpreted/perl/

nebo si vyberte ve studnici vědění

google: perl tutorial site:cz
google: perl pdb basic
http://www.perlmol.org/ ...

tj