Many existing data protection methods (DPMs) tend to overlook the use cases of forecasting, that can result in noisy forecasts and decreased accuracy. To address this issue, we develop a novel data protection framework for data providers who prioritize data privacy with minimal loss of forecasting accuracy. Recent studies show that k-nearest Time Series (k-nTS) swapping can be applied to achieve usable forecasting accuracy while maintaining an acceptable level of privacy. However, this method is often times inefficient because it requires a few tasks for each of time series data, say, (A_j in mathbb {R}^n) for all (j=1, dots , J), where J is very large: (i) distance computing for all (J-1) vectors from (A_j); (ii) ordering all (J-1) vectors based on the calculated distances from (A_j); iii) randomly selecting one of them to swap the values; (iv) and then, finally, moving to the next time series, to start over the entire process until the last vector (A_J). This entire process is time-consuming, inefficient for practical applications, and lacks any performance guarantee. In this research, we propose a new method called the K-means Time Series (K-mTS) Shuffling, which can be applied to thousands of time series to enhance efficiency and ensure performance guarantees. The efficiency is improved by (i) clustering, which eliminates the need for iterative individual distance calculations, and (ii) simultaneous data swapping (or shuffling) within the same cluster. Furthermore, for performance guarantee, instead of using random selection as in the k-nTS approach, we apply the perfect matching scheme to find optimal vector to swap with for each (A_j) that allows data providers to systematically manage the trade-offs between data privacy and forecasting accuracy. The presented numerical results indicate that our proposed method achieves guaranteed usable forecasting accuracy than that obtained using confidential data protected by traditional methods, while maintaing an acceptable level of privacy.