We calculate both analytically and numerically the ac susceptibility chi(omega) and the nonlinear electromagnetic response of thin superconductor strips and disks of constant thickness in a perpendicular time-dependent magnetic field B-a(t) = B(0)cos omega t, taking account of the strong nonlinearity of the voltage-current characteristics below the irreversibility line. We consider integral equations of nonlinear nonlocal flux diffusion for a wide class of thermally activated creep models. It is shown that thin superconductors, despite being fully in the critical state, exhibit a universal Meissner-like electromagnetic response in the dissipative flux-creep regime. The expression for the linear ac susceptibility during flux creep appears to be similar to the susceptibility of Ohmic conductors, but with the relaxation time constant replaced by the time t elapsed after flux creep has started. This result is independent of any material parameter or temperature or de field. For omega t much greater than 1, we obtain chi(omega)approximate to-1+pIn(qi omega t)/(i omega t), where p and q are constants. Above a critical ac amplitude B-0=B-l, the local response of the electric field becomes nonlinear, and there are two distinctive nonlinear regimes at B-0>B-l, where B-l similar to s(d/a)B-1/2(p), B-p is a characteristic field of full flux penetration, s(T,B)=\dInj/dInt\ is the dimensionless flux-creep rate and d and a are the sample thickness and width, respectively. For B-l<B-0<B-h(omega) the response of the electric field is strongly nonlinear but nonhysteretic, since the ac field B-a(t) does not cause a periodic inversion of the critical state. As a result, the magnetic moment exhibits a Meissner-Like nondissipative response, in stark contrast to the Bean model. For B-0>B-h(omega) the ac field causes hysteresis dissipation due to a periodic remagnetization of the critical state that gives rise to the hysteretic magnetic response of the Bean model at B-0 much greater than B-h. Here B-h(omega) weakly depends on omega and is of order (d/a)B-1/2(p) for a very wide frequency range, well below the irreversibility field, where s(T,B)much less than 1. Magnetization and ac losses at B-0 much greater than B-h are calculated accounting for the nonlinearity of E(J) at J<J(c) and a crossover between flux flow and flux creep at J similar or equal to J(c). All these regimes were confirmed by our computer simulations of nonlinear flux diffusion in strips and disks.